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Introduction to Formal Realizability 
Theory—I 


By BROCKWAY McMILLAN 
(Manuscript received October 15, 1951) 


This paper offers a general approach to the realizability theory of net- 
works with many accessible terminals. The methods developed are applied 
to give a complete characterization of all finite passive networks. 


I. SUMMARY 


1.0 A principal result of this paper is to characterize those matrices 
Z(p), functions of the frequency parameter p, which can be realized as 
open-circuit impedance matrices of finite passive networks. This char- 
acterization is provided by the following theorem: 


1.1 Theorem:* Let Z(p) be an n X n matrix whose elements are Z,;(p), 
1< r,s < n, where 
| (i) Each Z,s(p) is a rational function 

(ii) Zi s(p) = Z,;(p) (the bar denotes complex conjugate) 

(ili) Zs(p) = Zsr(p) 


(iv) For each set of real constants ki , --- ,k, , the function 


el) = DL, Zp) 


has a non-negative real part whenever Re(p) > 0. 
Then there exists a finite passive network, a 2n-pole, which has the 
impedance matrix Z(p). 


* Presented to the American Mathematical Society, April 17, 1948. Abstract 
260, Bulletin of the A.M.S. No. 54, July, 1948. 
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Conversely, if a finite passive 2n-pole has an impedance matrix Z(p), 
that matrix has the properties (i), (ii), (iii), (iv). 

A formally identical dual theorem holds for open-circuit admittance 
matrices Y(p). 


1.2 A general realizability theorem, applicable to and characterizing 
completely all finite passive networks, whether having impedance ma- 
trices or not, is also proved. 


1.3 An effort is made to lay a foundation adequate for the realizability 
theory of both active and passive multi-terminal devices. To this end, 
a large part of the paper is devoted to the scrutiny of fundamental 
properties of networks. 


II. INTRODUCTION AND FOREWORD 


2.0 Network theory provides direct means for associating with an 
electrical network a mathematical description which characterizes the 
behavior of that network. Typically, this results in shifting engineering 
attention from a detailed, possibly quite intricate, electrical structure 
to a mathematical entity which succinctly describes the relevant be- 
havior of that structure. An essential feature of this shift in focus is 
emphasized by the word ‘relevant’: only those terminals of the net- 
work which are directly relevant to the problem at hand are considered 
in the mathematical description. Design work can then be done in 
terms of constructs relating explicitly to these accessible terminals, the 
effect of the internal structure being felt only by implication. 

The physical origins of these mathematical constructs, and the im- 
plications of the internal structure upon them, cannot however be en- 
tirely forgotten, for they have mathematical consequences which are 
not always immediately evident. Until he knows these limitations— 
imposed upon him by the physical nature or the necessary structural 
form of the networks he is designing—a design engineer cannot make 
free use of the mathematical tools that network theory has provided. 

We give the name “realizability theory” to that part of network 
‘theory which aims at the isolation and understanding of those broad 
limitations upon network performance, i.e., upon the mathematical 
constructs which describe that performance—which are imposed by 
limitations on the network structure. One may also include in the 
province of realizability theory some of the converse questions: the 
study of those structural features common to all networks whose per- 
formance is limited in some specified way. 

Realizability theory would have little content were it not that “per- 


FORMAL REALIZABILITY THEORY—I 219 


formance” here must be construed to mean performance as viewed from 
the accessible terminals only. Were all branch currents and node poten- 
tials in a network available to observation, a mathematical statement 
of performance would be equivalent to stating the full system of dif- 
ferential equations governing these quantities, i.e., equivalent to giving 
the detailed network diagram. 


2.1 With a few important exceptions, the converse kind of problem in 
realizability theory does not lead to a strict implication from functional 
limitations to structural features, because the field of equivalent struc- 
tures for a specified performance is very broad. Typically, it is only by 
imposing some general a priori limitations on structure that further 
conclusions can be firmly drawn from a functional limitation. In study- 
ing this kind of problem one is rapidly led from those basic issues which 
are clearly part of realizability theory toward general, difficult, and 
usually unsolved problems of network synthesis. One cannot, and should 
not, draw a sharp boundary here, but Nature so far has provided a 
fairly definite one for us, in that most of these problems have proved 
too difficult of solution. 


2.2 The direct realizability problems, the passage from structural prop- 
erties to functional properties, have been somewhat more tractable. 
Here, again, there-is no clear dividing line between general realizability 
theory and the sort of design theory in which, for example, one specifies 
a particular filter structure depending on a limited number of param- 
eters and examines the performance of the structure as a function of 
these parameters. There is an extensive literature at or near this latter 
level of generality, most of it relating to filters or filter-like structures 
(e.g., interstage couplers in amplifiers). 

At a more basic level, the limitations on a network’s structure which 
are commonly met in practice are of the following kinds: 

a. Limitations on the kind of elements appearing, e.g., to passive 
networks, networks without coupled’ coils, networks whose elements 
have specified parasitics, etc; . 

b. Limitations on the general form of the network diagram, e.g., to 
ladder or lattice structures, without limitation to a specified number of 
elements or parameters. 

Here the problems are varied and difficult. We survey briefly the 
present status of some of them. 


2.3 Networks with two accessible terminals, two-poles, are basic in 
network technology. Fortunately, also, two-poles are unique among 
networks in that there is always a simple way to describe their perform- 
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ance. Except for the trivial limiting case of an open circuit, every two- 
pole has a well-defined impedance, Z(p), a function of the complex 
frequency parameter p, which describes its performance in a way which 
is by now well understood. Dually, except for the limiting case of a 
short circuit, every two-pole has a well-defined admittance function 
Y(p). Even the limiting cases are tractable: every open circuit has the 
admittance function Y(p) = O and every short circuit the impedance 
function Z(p) = 0. 

In other words, by exercising his option to speak in terms either of 
impedance or of admittance, one can always describe the performance 
of a two-pole by using a single function of frequency. 

The descriptive simplicity and practical importance of two-poles led 
early to a fairly complete realizability theory for them. In 1924 R. M. 
Foster’ gave a function-theoretic characterization of the impedance 
functions of finite passive two-poles containing only reactances. The 
corresponding problem for two-poles which are not at all limited as to 
structure, beyond being finite and passive, was solved by O. Brune’ in 
1931. The effects of various structural limitations have since been 
studied by several writers (cf. Darlington,’ Bott and Duffin’’). 


2.4 Technology, and the promptings of conscience, have meanwhile 
urged the study of devices with more than two accessible terminals. 
Here, however, Nature has been less kind, in that no uniquely simple 
method is available for describing the performance of such devices as 
viewed from their terminals. 

Indeed, basic network theory has been remiss here, in not even mak- 
ing available a mode of description which is generally applicable— 
whether simple or not. 

W. Cauer’ showed that, when one admits ideal transformers among 
his network components, it is sufficient to study networks which are 
natural and direct generalizations of two-poles, namely, 2n-poles,* for 
arbitrary values of n. The corresponding natural generalization of the 
impedance function Z(p) of a two-pole is the impedance matrix of a 
2n-pole: just as one multiplies a scalar current by a scalar impedance to 
get a scalar voltage, one multiplies a vector current by an impedance 
matrix to get a vector voltage. 

2.41 Not all descriptive difficulties are resolved, however, by consider- 
ing 2n-poles and their impedance or admittance matrices. For the 
moment, a simple example will suffice to show this: the 2 X 2-pole which 
consists simply of one pair of short-circuited terminals and one pair of 


* Defined in Cauer,® and also later here. 
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open-circuited terminals is a finite passive 2n-pole (n = 2) which has 
neither an impedance matrix nor an admittance matrix. 


2.42 When one eliminates this kind of descriptive difficulty by fixing 
his attention only upon 2n-poles for which an impedance matrix (or, 
dually, an admittance matrix) is available, the general realizability 
problem for finite passive devices is solved. A partial solution, for the 
case n = 2, was published by C. M. Gewertz*® in 1933. The solution 
(Theorem 1.1) of the problem for a general value of n has been an- 
nounced recently by three authors, independently: Y. Oono,’° in 1946,* 
the present author, in 1948, and M. Bayard, in 1949. The problem 
for reactive 2n-poles is much simpler and was solved by Cauer,® in 
1931. 


2.5 Intermediate between the fairly specific problems of filter theory on 
the one hand and the general realizability theory of multi-terminal 
devices on the other, lies the study of four-poles as transducers. There 
is a considerable literature on the realization of transfer functions or 
transfer impedances under various structural limitations. The basic 
and extensive work of Bode™ on active systems belongs also in this 
category since it is avowedly limited to single-loop structures. 


2.6 Beyond the important result that, by sufficiently elaborate circum- 
ventions, one may avoid the use of transformers in the synthesis of any 
two-pole, (Bott and Duffin’’) little in general is known about networks 
which do not have transformers. 


2.7 The present paper is an essay in the realizability theory of devices 
with many accessible terminals. Ideal transformers are admitted as 
network elements; indeed, their use is essential. This fact is indicated 
by the adjective “formal” appearing in the title. 

The availability of ideal transformers makes it possible to exploit 
the simplification noted by Cauer and to consider only networks which 
are 2n-poles in his sense. The aim of the paper, therefore, is to set a 
foundation for realizability theory which is completely general within 
this framework. 


2.71 The first problem is that of description. We observed above an 
example—entirely trivial—of a passive four-pole which had neither an 
impedance nor an admittance matrix. Unfortunately, opportunities 


* Date of Japanese publication. This reference, and manuscript of Oono!® 11, 
were sent by Professor Oono in a personal communication to R. L. Dietzold, who 
showed them to me in December, 1948, while an early draft of the present paper 
was in preparation. The recent (1950) American republication of Oono!® unfor- 
tunately omits reference to the original. 

+ Cf. footnote to 1.1. 
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for this kind of degeneracy become manifold in multi-terminal devices, 
and some degree of degeneracy is the rule rather than the exception. 
Consider an entirely practical example: that of an amplifier chassis from 
which the tubes have been removed.* Here the degeneracy is essential 
and intrinsic; 1t would be highly artificial to regard it as the mere acci- 
dent of a limiting case. True, given any particular degenerate network, 
there is usually evident a method for representing or describing its be- 
havior. What is needed, however, is a mode of representation which is 
applicable generally to any 2n-pole without a priori knowledge of its 
structure or peculiar degeneracies. 


2.72 The mode of representation adopted in this paper, embodied in 
the notions of general 2n-pole (Section 4) and linear correspondence 
(Section 6), is an obvious one, and so completely general that it solves 
no problems other than the elemental one for which it was introduced. 
It provides a definite mathematical construct whose properties one can 
discuss with mathematical precision. This is all that we ask of it. 

Realizability theory begins and ends with the study of these proper- 
ties. It would be more accurate to say that the notion of general 2n-pole 
describes a particular, but still very large, class of mathematical en- 
tities; realizability theory consists in the study of certain subclasses of 
the whole class of these entities, the particular subclasses being distin- 
guished by special, and to us interesting, properties. 


2.73 Despite its avowed aim at generality, the paper is oriented toward 
the realizability theory of finite passive networks. It ultimately provides 
a proof of 1.1 and indeed a complete characterization of finite passive 
2n-poles, however degenerate. This characterization is accomplished in 
a sequence of postulates, each one delineating a property of general 
2n-poles, 1.e., a subclass consisting of all 2n-poles having this property. 
The class of 2n-poles having all of these properties is:then identified 
with the class of 2n-poles obtained from finite passive networks. 


2.74 If we have succeeded here in our hope to set an adequate founda- 
tion for the realizability theory of devices with many terminals, it will 
be because of the nature and organization of the postulates themselves. 
They describe what at present seem to be individually significant prop- 
erties of 2n-poles, of progressively greater specificity, which in the 
aggregate characterize finite passive devices. By eliminating them in 
various combinations one obtains larger classes of objects. Further re- 

* It is exactly this example, and the practical need of an adequate theory for 


it, which led the author first to study the realizability theory of passive multi- 
terminal devices. 
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search alone will tell whether or not one obtains in this way the kinds 
of device which are significant. For example, one would like general 
realizability theorems for structures containing vacuum tubes with 
frequency-independent transconductances, vacuum tubes with non- 
vanishing transit times, unilateral devices with specified parasitics, ete. 


2.75 Actually, the postulates as we have given them are certainly not 
adequate for such an ambitious program. Exigencies of the presentation 
have dictated a number of condensations and compromises. It is hoped 
that the basic ideas are still evident even if not isolated individually in 
separate and entirely independent postulates. In any event, it is the 
author’s firm belief that the presentation as given is at least illustrative 
of the kind of approach, and the level of mathematical detail, which 
will be needed if one is ever to provide a truly adequate realizability 
theory: a theory which will cover, for example, the broad range of active 
linear systems which present-day technology allows us to consider. 


2.8 Apart from the network theoretic concepts, which must be evalu- 
ated by their effectiveness in solving problems—an assessment which is 
by no means yet complete—this paper is strongly marked by an idio- 
syncracy of its author: a consistent and insistent use of geometric ideas 
and terminology. This is based on the personal experience that linear 
algebra achieves logical unity and a freedom from encumbering notation 
when viewed in this way. A general reference covering most of the linear 
algebra (geometry) required here is P. R. Halmos’ elegant monograph’. 


2.9 For a proof solely of 1.1, which has already been three times proved 
in the literature," “°'' this paper provides an apparatus which is too 
cumbersome. There is even a sense in which 1.1 alone provides a charac- 
terization of all finite passive devices, for it seems to be generally ac- 
cepted that, by the use of ideal transformers, any finite passive network 
can be represented as a network which has an impedance matrix to 
which is adjoined suitable ideal transformers. Therefore we cannot claim 
that, in using this cumbrous apparatus to characterize all finite passive 
2n-poles (including the degenerate ones), we have offered anything not 
already provided by a simpler proof of 1.1. 

Three things may be said in rebuttal. First, we have already empha- 
sized that the apparatus here exhibited was designed for more problems 
than that to which it is here applied. It is presented in the belief that 
it will prove of further use. 

Second, even in the study of passive networks, it has seemed to the 
author helpful to look at the manifold things which are not passive net- 
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works. One gets then a clearer view of the unique Position occupied by 
passive devices among all linear systems. 

Third, there is a kind of semantic issue here: the assertion that any 
finite passive network (sic) can be put in such a form that 1.1 applies 
seems to this author to give a kind of circular characterization of such 
devices. A characterization which did not itself involve the concept of 
a network seems more satisfying. Logically, there is no circle here, but 
this is a fact requiring proof. A careful reading of this paper will show 
that it provides a proof. This particular subtlety does not of itself justify 
the lengths to which we have gone. It is, however, no longer a subtlety 
if one wishes to consider devices which do not have a representation in 
terms of something non-degenerate to which ideal transformers have 


been added. 


2.91 The present Part I of the paper is so organized that at the end of 
Section 8 the reader is in possession of all of its principal results and its 
basic ideas. The remaining Sections, 9 through 20, may then be regarded 
as an Appendix containing the details of proofs. Indeed, Part II will be 
largely devoted to further details of proof, though there will be there 
one important idea not mentioned, save casually, in Part I—the idea 
of degree for a matrix. 

In Sections 4 through 11, (ene paragraphs have been distinguished 
from explanatory or peunikie ones by starring the paragraph numeral. 

Part II of the paper contains the bulk of the proof of 1.1. This proof 
is modelled after that of Brune’ for the realizability of two-poles. One 
familiar with the Brune process will probably find Part II readable 
without extensive reference to Part I. 

Let the reader be warned that the Brune process is not a practical 
one for realizing networks because of its critical dependence upon a 
difficult minimization and balancing operation. The same criticism 
applies to the generalized Brune process of Part II. 

The Brune process is of theoretical importance because it does realize 
a network with the minimum number of reactive elements. These facts 
will be brought to light in Part IT. - 

The proofs of Oono” and Bayard’ are different from ours. That of 
Oono" again follows the Brune model. 


III, INTRODUCTION TO PART I 


3.0 We keep before us first the problem of finding a mathematical de- 
scription applicable to and characterizing the behavior of all finite pas- 
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sive networks. Second, we seek to make mathematically precise those 
ideas which appear to form the basis of general realizability theory. 
Sections 4 through 7 introduce the immediate mathematical machinery 
for this. Section 8 states the fundamental realizability theorem and 
outlines its proof. At this point the reader has had an introduction to 
the results of the paper. The remainder of the paper is then devoted to 
the technical details of proof. Beginning with Section 12, the device of 
starring the technical passages will be dropped. 


3.1 Cauer’ distinguished precisely the class of networks called 2n-poles 
from the class of all multi-terminal networks. He also showed that, by 
the use of ideal transformers, any multi-terminal network is equivalent 
to a network which is a 2n-pole (for some n) in his sense. We shall in 
Section 4 define a class of objects to be called general 2n-poles. This 
class includes all electrical networks which are 2n-poles in Cauer’s sense. 
Its definition abstracts the significant properties isolated by Cauer. 

For the study, alone, of finite passive networks, this definition is 
unnecessary, since one can in fact so put the arguments as to deal only 
with 2n-poles which are finite passive networks, and therefore to deal 
only with concepts already defined in Cauer’. The somewhat physical 
notion of a general 2n-pole is a convenient backdrop against which to 
display the important physical properties of finite passive networks, 
and, indeed, of networks in general. Having it available, we use it 
throughout the realizability arguments. 


IV. DEFINITION OF GENERAL 2n-POLE 


4.0* Network theory establishes a correspondence between oriented 
linear graphs and systems of differential equations. With each node of 
the graph is associated a potential #, = E,(¢) and with each oriented 
branch acurrent J, = J,(¢). These potentials and currents are constrained, 
first by Krichoff’s laws, and second by differential equations which de- 
pend upon the nature of the branches but not upon the topology of the 
graph. 

~4.01* A finite passive network is one whose graph has the following 
properties: 


(i) There are finitely many nodes, 1, 2,-:-: , N. 
(ii) There are finitely many branches, 1, 2, --: , B. 


* Technical paragraph as explained in Section 2.91. 
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(in) Let the b-th branch begin at node n, and end at nm, . Let vw» = 
i, — E,, . Then for each 6, one of 


v5 = eli, R,> 0 (a) 
d 
=O, C,>0 (b) 
T , 
vo = > Loy: oe (c) 
holds, where the matrix J, is real, symmetric, and semi-defi- 


nite. 


Cauer has shown’ how an ideal transformer can be defined as the 
limiting case of a finite passive network. It is indeed no more nor less 
ideal than an open circuit (2p = © or Cy = 0) or a short circuit (Ry = 0 
or Ch = ©). 


4.02 We seldom deal with networks in the detail which is implicit in 
(iii) above. We are usually interested in the external characteristics, so 
to speak, of such networks as viewed from a relatively small number of 
terminals (nodes). These multi-terminal devices, however, we continue 
to incorporate into larger network diagrams. It is usually clear how 
Kirchoff’s laws are to be applied in these cases, and what the differential 
equations of the resulting system are. We are obliged, however, to make 
these matters precise before we can deal intelligently with the most 
general physical properties of networks. 


4.1 We have seen the two kinds of constaint that a multi-terminal de- 
vice imposes on the branch currents and node voltages in a network in 
which it is incorporated: the topological ones contained in Kirchoff’s 
laws and the dynamical ones described by differential equations. Cor- 
respondingly, there are two aspects to the concept of general 2n-pole. 


4.11* In its relation to Kirchoff’s laws, a general 2n-pole is indicated as 
an object with n pairs of terminals (7, , 7), 1 <r < n. Each terminal 
can be made a node in an arbitrary finite diagram constructed out of 
network elements and other general 2m-poles, with arbitrary values of 
m. This diagram is not an oriented linear graph, so we have no basis 
for the use of Kirchoff’s laws. From it, however, we construct an ori- 
ented linear graph, called the zdeal graph of the diagram, by the follow- 
ing rule: 

The nodes of tiie ideal graph are those of the original diagram. Every 


* Technical paragraph as explained in Section 2.91. 


i) 
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oriented branch of the original diagram is repeated in the ideal graph, 
similarly situated and oriented. Between those nodes which, in the 
original diagram, were the (7, , T) of a 2n-pole N, is drawn a branch 
6, , called the r-th ideal branch of N, oriented from T, to T,. This is 
done for each such terminal pair. 

Kirchoff’s laws now apply to this ideal graph. 


4.12* Consider a particular 2n-pole N. Let ZH, be the potential of 7, , 
E, that of T. . Define 


v(t) = E, — E,. 


Then z,(¢) is the voltage across the ideal branch 8, so oriented that 
v,(t) > 0 when 7’, is positive relative to 7; . Let k,(t) represent the cur- 
rent entering 7’, . Then k,(t) = J,(t), the current in 8, , so k,(t) is also 

the current leaving 7, . This is the force of the notion of ideal branch 
~ and the fact which distinguishes a network which is a 2n-pole from an 
arbitrary network with 2n terminals. 


4.13 For example, the network at (a) of Fig. 1 is not a2 X 2 pole because 
its currents are not constrained to meet the ideal branch requirement. 
The addition of ideal transformers in either of the ways shown in (b) 
or (c) of the figure converts it to a 2 X 2 pole. Of course in a circuit in 
which the currents are constrained externally—as they would be, for 


T T2 
T/ T/ 
1 (a) 2 
gl tts elle 
3 pe 2 3 
T’ Tz 
(b) ‘ 


T 


Ty 





(Cc) 
Fig. 1—Conversions of a four pole to a 2 X 2 pole. 


* Technical paragraph as explained in Section 2.91. 


228 THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1952 


example, when the 2 X 2 pole is driven by separate generators in the two 
external meshes—these transformers can be eliminated. The definition 
of 2n-pole requires however that in every context the ideal branch con- 
cept is valid. 


4.2* The second aspect of the concept of general 2n-pole is that it im- 
poses some kind of constraint—other than that implied by 4.11 and 
Kirchoff’s laws—upon the voltages across and currents in its ideal 
branches. Define the symbols 


v = v(t) = [ni(Z), w(t), +--+ , on] 
and 


as the n-tuples, respectively, of voltages across (T, , T.) and currents 
into T,,1 <r <n. These are added and multiplied by scalars by the 
usual rules of vector algebra. If vy and & represent s¢multaneous values 
of voltage and current in the 2n-pole N—i.e., values satisfying all the 
constraints—then we say that N admits the pair |v, k]. 

We say that N admits v if there is a & such that N admits the pair 
[v, kj. This & is said to correspond to v. Dually, N admits & if there is a 
v (corresponding to k) such that N admits [v, &]. 

The constraints imposed by a general 2n-pole N on voltages and cur- 
rents are completely described by the totality of pairs [v, k] which N 
admits. We shall define a general 2n-pole, therefore, as 


(i) a collection of n oriented ideal branches, as in 4.11, and 


(ii) a list of pairs [v, &] of voltages and currents admitted in these 
branches. 


Hereafter we shall usually drop the adjective “general.” 


4.21 The definition we have just given is, in a way, a postulational 
form of an n-dimensional Thevenin’s theorem. It postulates that a 
2n-pole is a thing} which, as far as the outside world is concerned, can 
be represented by a collection of two-poles 6, , 1 <7 < n, among which 
there exists a complicated agreement as to what currents and voltages 
will be admitted. 


4.22 The passive networks (b) and (ec) of Fig. 1 define 2 X 2 poles, be- 
cause they satisfy 2.01 and clearly permit a complete specification of 
the admissible pairs [v, &]. Any equivalent network would also specify 


* Technical paragraph as explained in Section 2.91. 
} In fact, at this level of generality, any thing. 
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the same 2 X 2 pole, because—by its very equivalence—it would admit 
the same pairs. The closest association we can make between a 2n-pole 
and a network, then, is to identify the 2n-pole with an equivalence 
class of networks. 


4.23 The completely symmetric role played by voltages and currents 
in this definition of general 2n-pole will make it possible to take early 
advantage of the well-known duality principle of network theory. We 
shall do so freely. 


4.3* We shall call a 2n-pole physically realizable if its admissible pairs 
[v, k| are the solutions of a system of differential equations obtained 
from a finite passive network, admitting the limiting elements: ideal 
transformers, open circuits, and short circuits. 


V. PHYSICAL PROPERTIES OF NETWORKS 


5.0 There are clearly a great many properties of finite passive networks 
which are not yet possessed by the general 2n-poles now introduced. It 
is instructive to examine these properties physically. 


5.1 In the first place, the dynamical constraints (a), (b), and (c) of 4.01 
are expressed by linear, time invariant, differential equations. Accord- 
ingly, the 2n-poles of network theory are: 


5.11 Linear, in that the class of admissible pairs [v, &] 1s a linear space; 


5.12 Time invariant, admitting with each [v(t), k(] also all 
[u(t + r), k(t + 7)] for aribtrary r. 


5.2 In the second place, a physical network N cannot predict the future, 
i.e., it cannot respond before it is excited. This can be formalized in 
terms of the pairs [v, k] admitted by N, but to do so would require some 
digression. The reasons will be seen under 5.7 below. 


5.8 We have already mentioned the constraints imposed on voltages 
and currents in a network by the topology of the network, through the 
medium of Kirchoff’s laws. These constraints have three important 
properties: 

5.381 They are workless, since they are imposed by resistanceless 


connections, leakless nodes, and, in the formal theory, by ideal 
transformers. 


5.32 Though it seems scarcely necessary to say it, they are the only 
workless constraints. All other constraints are dynamical and have 
powers or energies associated with them. 
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5.338 They are frequency independent, that is, holonomic in the sense 
of dynamics. 


5.4 The workless and the dynamical constraints in a physical network 
are all defined by relations with real coefficients. The space of admissible 
pairs is then a real linear space. . 


5.5 The positivities specified in 4.01 are characteristic of passive sys- 
' tems. They correspond to the fact that the power dissipation and the 
stored energies are all positive. 


5.6 By definition, finite passive networks contain finitely many lumped 
elements. Correspondingly, their resonances and anti-resonances are 
finite in number. 


5.7 We are accumstomed to dealing with networks which have, in addi- 
tion to the properties listed above, a kind of non-degeneracy, in that 
the list of admissible pairs [v, &] satisfies: 


5.71 At least one of v or k can be specified arbitrarily—any real function 
is admitted; 


5.72 When the free number of [v, &] is specified, the other is uniquely 
determined. 

For these non-degenerate networks, the property 5.2 above is easily 
formalized: if, say, k is determined by v, then | 


v(t) =v(t) for t< bt 
implies . 
BQ) =h@ for t<b, 


where [v’, &'] are admissible pairs, i = 1,2. The general statement of 5.2 
involves this condition and some discussion of the v’s for which N ad- 
mits [v,.0], and the dual notions. 


5.8 The reason for speaking in terms of pairs [v, &], instead of in terms 
of ‘‘cause” and “effect,” or “impulse” and “response,” is hinted at by 
5.7 above. For the tacit implications of the cause and effect language 
completely obscure the fact that 5.71 and 5.72 are properties which are 
not automatically possessed by electrical networks. In fact, the simple 
four-pole of 2.41—a pair of unconnected terminals 7, , 7, , and a pair 
of shorted terminals T; , 7;—has neither property, yet it is a perfectly 
good linear time invariant four pole. Its admissible pairs are 


Ia ? 0), (0, ke)I, 


where v,; and ky are arbitrary real functions of the time 
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VI. LINEAR CORRESPONDENCES 


6.0 In developing the formal properties of 2n-poles which are equivalent 
to the physical ones just listed, it would be instructive to adjoin re- 
quirements piecemeal, much in the order given in Section 5. Space does 
not permit us full enjoyment of this luxury, but the reader will find a 
rough parallel between Section 5 and the developments of this Section 
and Section 7. | 


6.1 It is well known that linear time invariant systems are best studied 
by the tools of Fourier or Laplace analysis. We make this fact the basis 
of our first step in characterizing physically realizable 2n-poles simply 
by phrasing our whole discussion in the frequency language. The con- 
tent of the following paragraph will be obvious enough, but it does de- 
fine terms to be used later. 


6.11* Let v and k, without underscores, represent n-tuples of complex 
numbers: 
v = [v1 ,%2,°°° , Val, (1) 
b= [hy . ke rence k,]. (2) 
These are to be manipulated by the rules of vector algebra. Let p be a 
complex number. We shall say that a 2n-pole N admits the pair [v, k] 


at frequency p, if in the sense of 4.2 N admits the pair [v, k] (with under- 
scores) where v has components 


v,(t) = Re(v,e”'), 1<r<n, (3) 
and & has components 
k(t) = Re(k,e?), l<r<n. (4) 


Also analogously to 4.2, we say that N admits v at frequency p if 
there is a k such that N admits [v, k] at frequency p, and that this k 
corresponds to v (at frequency p). Similarly, N admits k at frequency p 
if there is a (corresponding) v such that N admits [v, k] (at p). 


6.12* Let V denote the aggregate of all n-tuples (1), and K the aggre- 
gate of all n-tuples (2). These are then complex linear spaces. 


6.2* As our first step toward characterizing realizable 2n-poles, let us 
consider a linear correspondence L between V and K described by the 
postulates: 
Pl. There is a set IT, of complex numbers and for each peDy, a list 
L(p) of pairs [v, kK], veV, keK. 
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P2. If [v', k']eL(p) and [v’, k’JeL(p), then 
[aw + ay’, ak’ + ack’ |eL(p) 
for any complex numbers 4 , a . 


6.21* Given such a linear correspondence L, we can always describe a 
2n-pole Nz, by: 

Nz admits [v, k] at frequency p if and only if [v, k]eL(p). 

That is, we can always interpret the pairs [v, k] generated by (8) and 
(4) from the |v, kleL(p), for each pel, , as the voltages across and 
currents in a set of n ideal branches. We call N; the 2n-pole associated 
with LD. 


6.22* We call I, the frequency domain of L (or of Nz). 


6.23 From here on, the words “2n-pole”’ can with some strain be re- 
garded as suggestive but unnecessary. We in fact deal with linear corre- 
spondences—having properties as yet unspecified—and shall show how 
physical networks can be constructed which admit the pairs [v, k]eZ(p). 
Actually we use freely the concept of general 2n-pole and thereby avoid 
some elaborate circumlocutions. 

6.24* We identify two correspondences J, and ZL», as being the same if 
(i) their frequency domains differ only by a finite set, and (ii) for each 
p where both are defined the lists Z;(p) and L2(p) are the same. 

6.3 The simplest linear correspondences are those generated by ma- 
trices. For example, let Z(p) be ann X n matrix with, say, elements 
Zrs(p) which are rational functions of p, 1 < r,s < n. Let Ty, consist of 
all the values of p at which Z(p) is defined. For pel, , define L(p) as 
the class of all pairs 


[v, ki] | (5) 


obtained by letting k range over K, where for each k, v is defined by the 
matrix equation 


v = Z(p)k. (6) 


This kind of matrix equation will be used throughout to symbolize the 
nm component equations 


vr = >, Zrs(p) Ke, l<r<n. (7) 
s=1 


The list of pairs L(p) defined by (5) clearly satisfies Pl and P2. It 
can therefore be used to define a 2n-pole N;, . It is easy to see that N, 
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in fact is non-degenerate in a sense similar to that of 5.7, for the current 
amplitudes k can be specified arbitrarily, and the resulting voltage 
amplitudes v are then fixed by k and p, by (6). 

Z(p) is called the impedance matrix of the 2n-pole N,. It is also 
sometimes called the open-circuit impedance matrix, because each 
Z,s(p) is, by (7), the voltage amplitude across (7, , T’.) when the current 
amplitudes at all terminals save (7, , T’,) are zero—i.e., when all pairs 
save the s-th are on open circuit. 


6.31 Dually, the pairs 
[v, Y(p)o] 


defined by an admittance matrix Y(p) as v ranges over V define a linear 
time invariant 2n-pole which is non-degenerate. 


VII. WORK AND ENERGY 


7.0* A linear correspondence satisfying Pl and P2 is something which 
abstracts the properties of linearity and time invariance. Most of the 
remaining properties of physical networks involve the mention of work 
or energy. These concepts enter our picture by way of the scalar product 
(v, k) between a voltage n-tuple (1) and a current n-tuple (2), of 6.11. 
This scalar product is defined by 


(v,k) = >> vk. (1) 
r=1 
7.01 If p = ww, one easily calculates from (3) and (4) of 6.11 that 
1 T n 
2 Re(v, k) = lim — [ - sshd | dt. 
Too 27 —T Lr=1 


That is, when p = 2m, the real part of 2(v, k) measures the average total 
power dissipated by the system of currents k,(¢) against the driving 
voltages u,(¢). 

When p is not a pure imaginary, the interpretation of the scalar 
product (v, k) is not so clearly physical as this. The reader will ulti- 
mately observe that our significant statements about such products can 
all be reduced to statements applicable when p = iw, 1.e., when the 
power interpretation is valid. 


7.1* An important concept in what follows is that of the annihilator of 
a linear manifold (Halmos’, par. 16). Let Vi € V be a linear manifold. 
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Then its annihilator (V;)° is the set of all & such that 
veV; implies (v, k) = 0. 


(V;)° is a linear manifold in K. 


Dually, given K; C K, (K,)° is the linear manifold of all veV such 
that 


keK, implies (v, k) = 0. 


The annihilator concept is the analog in our general geometric frame- 
work of the idea of orthogonality. It clearly suggests a connection with 
workless constraints. 


7.2* The complex conjugate of an n-tuple v (or k) is defined in the 
obvious way: if 


v= [r,s , nl 
then 
j= [a1 , Peed Vals 


This conjugation operation clearly has the properties 


———-  _. (2) 
ag + by = a& + bF 


where a and 0 are scalars and & and 7 are (consistently) elements of V 
or K. Furthermore, at once from (1) of 7.0, 


(v, k) = (, k). (3) 


7.21* A linear manifold will be called real if it contains, with any 
n-tuple also the conjugate of that n-tuple. 


7.22* A real manifold is spanned by real n-tuples. This will be pronee 
in the Appendix, Section 20. 


7.23* The annihilator of a real manifold is real. For let Ky be real and 
k’, «++, k" be real n-tuples which span K, . Then if ve (K,)° every 


(v, k’) = 0, 


and conversely. But then also 


(@, k*) = 0, *) = 0 = 0, 





sO be(K,)°. 
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7.3* Given a linear correspondence L, we make several definitions: 
V.(p) is the set of all veV such that there is a k with [v, k]eL(p). 
K_(p) is the set of all keK such that there is a v with [v, k]eL(p). 
Vio(p) is the set of veVz(p) such that 

[v, O]eL(p). 
Kyo(p) is the set of keK,(p) such that 
[0, k]eL(p). 


7.31* The postulate P2 implies that for each Pela Vi(p), Kx(p). 
Vio(p) and Ky, (p) are all linear manifolds. 


7.32 V.(p), for example, is the set of veV such that N, admits v at 
frequency p. 


7.4* We now postulate 
P38. There exist fixed linear manifolds Vz, Cc V, K, © K such that 


(A) For every pel’, , Vi(p) = Vi = (Ki(p))° 
(I) For every pel’, , Ki(p) = Ki = (Vio(p))”. 


7.41* We may henceforth write Vio , Kio , for Vio(p), Kro(p), knowing 
that, under P3 


Vio = (K:)’, 
Kio aoe Vi): 


7.42 Linear correspondences satisfying P3 abstract the properties men- 
tioned in 5.8. The equalities Vi(p) = Vz, Ki(p) = Ky, guarantee the 
frequency-independence of the workless constraints. The equalities 
Vilp) = (Kiw(p))’, Ki(p) = (Vio(p))° in a sense guarantee that the 
only constraints imposed upon admissible currents and voltages (as 
opposed to constraints relating currents and voltages) are those which 
arise from open or short circuits, l.e., are workless. 


7.48 An illustrative consequence of P3, for example, is that if L satisfies 
P3 and if N, is such that all of the current amplitudes can be specified 
arbitrarily, then indeed the voltages are determined by the currents. 
This will appear as a consequence of 8.1. It is a very general theorem 
about networks of a kind that this author, at least, has not heretofore 
encountered. 
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7.5* Continuing toward realizability, we introduce 
P4. If pel, , then pel. If [v, kJeL(p), then [d, k]eL(p). 


This postulate embodies most of the reality properties of networks. 
It has as an immediate consequence the 


7.51* Lemma: If L satisfies P1, P2, P3, and P4, then all of 
Vi, Vi, Ki, Kio 


are real. 
Proof: By P4, veVi(p) = Vx implies teVi(p) = Vz. Hence Vz is 
real. Then Kyo = (V,)° is real, and dually. — 


7.6* The three remaining postulates on Z refer to scalar products. 
They are concerned with the energy questions related to passivity, 
rather than with the workless constraint questions. 

P5. If [u, jJeL(p) and [v, k]eL(p), and if 


(A) wu and »v are real, or if 
(I) j and k are real, 
then. 
(u, k) = (2, j). 


7.61 This is the property which provides the reciprocity law. In its 
presence, the relations in P3 may be weakened to 


Vi(p) = Vi 2D (Ku(p))’, 

Ki(p) = Ki D (Vio(p))”. 
This fact will appear as a consequence of the lemma of Section 12. 
7.7* Lemma: A consequence of P2 and P3(A) is that if 

[v, kleL(p), r= 1,2, 
then for any ueVz , 

(u, ky) = (u, ke). 
For by P2 we have that a 
[vy — v, ky — ke] = [0, ki — keleL(p), 

hence ky — kyeKy» . Then however, by P3(A), weV, implies we(K1)", so 
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that 
0 = (u, ky = ke) = (u, ki) _— (u, ke). 

Q.E.D. A dual result follows from P3(I). 
7.71* The result of 7.7 above means that the scalar product (v, k) is 
fixed by v alone when we know that [v, k]eL(p). This means that, given 
veV, , there is a unique function F’,(p) defined for pel, by 

F,(p) = (2, k) 
where [v, k]eL(p). Dually, 

Ji(p) = (», ) 
is defined for each fixed keK, . 
7.72* (P6.) The complement of IT, is finite and 


(1) For each veVz , F,(p) is rationai 
(A) For each keK, , J;(p) is rational. 
7.73* (P7.) (A) Re(p) > 0 implies Re(F.(p)) > 0 
(I) Re(p) > 0 implies Re(Jz(p)) = 0. 


VIII. THE FUNDAMENTAL REALIZABILITY THEOREM 


\ 
8.0* We can now state our fundamental realizability theorem: If a 
linear correspondence JL satisfies Pl, --- , P7, the associated 2n-pole 
N, is physically realizable. Conversely, given a physically realizable 
2n-pole N, the associated linear correspondence satisfies Pl, --- , P7. 


8.01 Actually, the postulates Pl, --- , P7 are not unique nor even en- 
tirely independent. Many changes may be rung on them. We indicated 
one above. At the expense of apparent asymmetry, the (A) or (I) por- 
tions, in various combinations, can be deleted or weakened. We shall 
not pursue this subject further at this point, but must come back to it 
in Section 12. 


8.02 We close this Section by outlining the proof of 8.0. The details are 
then contained in the remainder of the paper. 


8.03 The proof that P1 through P7 are necessary for physical realiza- 
bility will be a direct one: it will be shown that, considered individually, 
each network branch and each ideal transformer satisfies the postulates. 
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By an application of Kron’s method (described by Synge”), it will then 
be shown that the imposition of Kirchoff’s laws preserves the postulates. 
This work is most efficiently performed after the full machinery of the 
sufficiency proofs is available, and will be done in Section 19. 


8.04 The sufficiency of P1 through P7 can be deduced—and we will do 
so—from the lemmas to be quoted below. Apart from Section 19 on 
necessity, the remainder of the paper is devoted to the proofs of these 
lemmas. 


8.1* Lemma: If L is a linear correspondence satisfying Pl, P2, P3, and 
P4, then there exists a fixed real nonsingular matrix W such that 


8.11 The list Z w(p) of all pairst 
(Wy, W’k], 
where [v, k]eL(p), describes a linear correspondence Dw satisfying 


Pl, P2, P38, and P4. 


8.12 The 2n-pole Nw (= Nz,,) associated with Ly consists of 


(i) Some number r of open-circuited terminal pairs (71, Ty), °-°, 


(T, ’ T.), 
(ii) Some number s of short-circuited terminal pairs (Tn—o41, Tr—e41)s 
“s (Tr T, )s 
T" 
(iii) A set oe =n — r — 8 terminal pairs (Ty41, Tri), +++, 
(Tas Dain) 


8.13 Either m = 0, or the terminal pairs in (iii) are those of a 2m-pole N, 
which has a nonsingular impedance matrix Zi(p). 
' This lemma, and the following, will be proved in 13.2. 


8.2* Lemma: If L satisfies P5, P6, and P7, then Zi(p) is a positive 
realt matrix, that is, Z:(p) satisfies (i), --- , (iv) of 1.1. 


8.3* Lemma: If a 2m-pole N, has a positive real impedance matrix, then 
Ni is physically. realizable. 

This is the sufficiency half of the matrix realizability theorem 1.1. 
Part II will be devoted to its proof. 


8.4* Lemma: If Nw is physically realizable, then N can be constructed 
from it by the use of ideal transformers. 

This is Cauer’s Transformation Theorem’ about which we shall say 
more in Section 9. 
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8.5* The sufficiency half of 8.0 is now clear. By 8.2 and 8.3, N; is physi- 
cally realizable. Clearly then Nw is, simply by the adjunction of the 
necessary open and short circuits. Finally N is by Cauer’s theorem, 8.4. 


8.6* We can see now how to prove the necessity of positive reality for 
the realizability of a positive real matrix Z(p). This is the necessity half 
of the matrix theorem 1.1. Let Z(p) be the matrix of a realizable N. 
Then N has an associated linear correspondence L satisfying P1, --- , P7, 
by the necessity half of 8.0. The pairs of D are the pairs 


[Z(p)k, ke] 


generated as k ranges over all n-tuples. By definition, then, the pairs of 
L w are 
. [WZ (p)k, W’k]. 


As k ranges over all n-tuples, the nonsingularity of W implies that 
W’k does also. Let U = W’. Then the pairs above are the same as 


[UZ(p)U’k, k] 


as k ranges over all n-tuples. Hence Lw has the impedance matrix 
UZ(p)U’, where U = W™ is real and nonsingular. Because Lw has an 
impedance matrix, r = 0 in 8.12. 

Now by 8.1 and 8.2, Zi(p) is positive real and the matrix UZ(p)U’ 
of Ly is just Z:(p) bordered by s rows and columns of zeros. It is then 
easy to see that UZ(p)U’ is positive real, and finally also that Z(p) is. 
These last two facts will be proved formally in Section 16. 


IX. CAUER’S TRANSFORMATION THEOREM 

: 5. P 
9.0 Cauer’s transformation theorem” is the cornerstone of formal reali- 
zability theory. In one form, the theorem reads: 


9.1* Let Z(p) be the impedance matrix of a physically realizable 2n-pole 
N. Let U be a real, constant, nonsingular matrix. Then 


UZ(p)U’ (1) 
is again the impedance matrix of a physically realizable 2n-pole, Nv . 
Nv can be constructed from N by the use of ideal transformers. 


9,2* A superficial generalization of this theorem can be obtained at once 
from Cauer’s proof. It asserts that if N is physically realizable and is 
described by the linear correspondence L, then there is a physically 
realizable 2n-pole Nw, obtainable from N by the use of ideal trans- 
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formers, which is described by the linear correspondence L w whose pairs 
at each : are the pairs 


[W'v, W'k], (2) 


where [v, k]eL(p). 
We refer to Cauer’ for the proof. It is straightforward. 


9.21 We shall use the second form (9.2) of Cauer’s theorem in our 
realization, process. Notice that it is in a sense a “physical” theorem, 
about the way one physical network is related to another. It is used in 
this way: we shall always solve a realizability problem by finding some 
network N which is easily realized, and then a W such that Ny, which 
is now realizable, provides a solution to the given problem. 


9,22* We shall call the 2n-pole N w a Cauer equivalent of N. 


9.3 Although Cauer’s theorem will be applied, in a sense, only a posteriori, 
its effect is fundamental. For it implies that formal physical realizability 
is a property of matrices which is invariant under the operation (1) 
or a property of correspondences which is invariant under (2). There is 
an extensive classical literature on the properties of matrices invariant 
under operations like that of (1), and the effect of Cauer’s theorem is to 
make these results all available to formal realizability theory. 


9.31* It is worth observing here that we are already well set up to use 
Cauer’s theorem: 
_ Lemma: If L is a linear correspondence satisfying P1,--- , P7, then 
the correspondence L w of 9.2 also satisfies Pl, --- , P7. 
Proof: Let M = Ly. FP1 and P2 for M are obvious, with Ty = T,. 
By definition of M/, 


Vulp) = WVi(p) = WV, 
Ku(p) = W'K.i(p) = W’K, 
Vuolp) = W°Vie(p) = WV 
K x0(p) = W’Kio(p) = W'Ki 


where W'S for a manifold S consists of all n-tuples W~’v, where veS. 
Hence in P38, » 


I 


Vulp) Vu = WV: 
K u(p) = Ky = W’K, 
for fixed manifolds Vy , K x as defined. 
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Now if veVio , then 
(v, k) = 0 


for every keK, = (Vi). Then, however, by direct calculation from 
Section 7.0, 


(Wy, W*k) = 0, 


where W* is the adjoint, i.e. transposed conjugate matrix of W. But 
because W is real, W* = W’. Hence if veVio, then 


(W~v, k) = 0 
for every keW’K, = Ky. Hence 
Ky = (WV 10)" sae (Vuo(p))’. 


By this and its dual, P3 is completed for M. 
The reamining postulates for M follow from those for L by the simple 
equality 


(v, k) = (W~'v, W’k) 
already established, combined with Ty = [z. 


9.32 For fixed Z(p), the matrices (1), as U ranges over a group, form an 
equivalence class. Classical matrix theory treats of such equivalence 
classes. This author’s predilection is to regard this theory from a geo- 
metrical point of view. In part this prejudice may be justified by the 
ease with which that slightly more general object, a linear corre- 
spondence, can be treated by geometrical methods. In any event we shall 
begin our program of proofs with a brief introduction to the geometrical 
approach. 


X. GEOMETRICAL PRELIMINARIES 


10.0* We now wish to consider V and K as complex n-dimensional 
linear spaces} respectively of voltage vectors v and current vectors k. 
The distinction here is in point of view. A vector v is regarded as an 
absolute geometrical object; an n-tuple {[v] = [a1 , --- , a,] 1s regarded 
as a set of coordinates for the vector v, relative to some coordinate basis. 
Given a fixed coordinate basis, there is a one-to-one correspondence 
between vectors v and the n-tuples [v] which represent them in that 
basis, a correspondence which preserves the operations of vector algebra. 
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10.01 The effect of attaching a geometric identity to vectors, rather than 
to n-tuples, is to make it possible to choose coordinate bases freely and 
as convenient, without elaborate constructions or even interpretations. 
We can then discuss properties of n-tuples (and other objects, e.g. 
matrices) which are invariant under the kind of operations exemplified 
by (1) and (2) of Section 9 as properties of a single geometric object, 
rather than as properties shared by an extensive class of concrete ob- 
jects which are converted into each other by the group of operations. 
10.1 This change in point of view need not change formally anything we 
have said to date; it simply erects a conceptual superstructure, or pro- 
vides a conceptual foundation, depending on the reader’s personal 
attitude. : 

We shall support this statement by going through the important ideas 
of Sections 4, 6, and 7 and examining their geometrical meanings or 
counterparts. It is convenient to consider first and at some length the 
notions of scalar product and complex conjugate. The geometric struc- 
ture will then be complete enough to permit a rapid survey of the 
remaining ideas. 


10.11* The geometrical counterpart of the scalar product introduced in 
7.0 is a numerically valued function ¢ = o(v, k) of two vector variables. 
Its first argument v ranges over V and its second argument k ranges over 
K. The function o(v, k) is linear in v and conjugate linear in k: 


a(au + bv, k) = ac(u, k) + bot, k), 
o(v, ak + bl) = do(v, k) + bo(v, £). 
We denote this function o(v, k) by the simple bracket notation (v, k). 


10.12 With this scalar product, the geometry of V and K is that of a 
space K and the space K* = V of conjugate linear functionals over K. 
This is analogous to the real geometry of space and conjugate space 
discussed at length in Halmos’. In fact, in the introduction to Chapter 
III of Halmos®, the modifications introduced by the conjugate linearity 
of (v, k) over K are treated in detail. 


(1) 


10.13* Because of its importance, we quote here a paraphrase of the 
results covered in Halmos’®, par. 12. 

(i) If f(v) is any numerically valued homogeneous linear function of 
veV, then there is a unique vector ky;eK such that . 


fv) = (, ky) 
for all veV. 
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(ii) If g(k) is any numerically values homogeneous conjugate-linear 
function of keK (i.e., if g(k) is linear in k) then there is a unique v,«V 
such that 


g(k) = (vy, k) 
for all keK. 


10.2* The annihilator (Vi)° of a manifold V; € V is, as in 7.1, the set of 
all keK such that 


veV; implies (v, k) = 0. 


10.21* It is shown in Halmos? that to each basis v', --- , v” in V there 
exists a unique dual basis k’, --- , k” in K such that 
(ak ie bres (2) 


where 6,; is the Kronecker symbol: 6,; = Oifr ¥ 8,6, =1,l<r,s<n. 
10.22 Tf 

[v] = [a1 , ~~ On| 

[K] — [b,, Mets bn] 


are the n-tuples representing v and k relative to a pair of dual bases, 
then it is easily computed from (1) and (2) that 


(v, k) = 3 Orb. | (4) 


(3) 


Therefore the concrete scalar product of 7.0 is indeed the geometric 
scalar product here considered, when we restrict our pairs of bases in 
V and K always to be dual in the sense of (2). 


10.23* We shall use the words “coordinate frame” or simply ‘‘frame”’ 
to denote a pair of dual bases in V and K. Any basis in V (or K) specifies 
a frame by the uniqueness result quoted above. 


10.24 We shall henceforth deal always with coordinate frames, in fact, 
ultimately, real coordinate frames, rather than arbitrary pairs of bases. 
This means in classical language that we are considering as ‘‘ geometrical 
properties” all properties which are preserved under the group of 
linear transformations which leave the bilinear form (4) invariant. 
The properties related to physical realizability will turn out to be 
invariant only under the subgroup of real linear transformations pre- 
serving (4). 
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10.3* Conjugation is an operation which to each veV associates a 
vector } uniquely determined by v with the properties 


Ti 
lI 
SS 


? 


eee ; (5) 
(au + bv) = aa + bi, 


where -a and b are any complex numbers and 4, 6 their conjugates. 


10.31* Given any such conjugation operation in V, and given any 
keK, define a function gi(v) by 


gv) = (@, k) (6) 
for veV. Then g;(v) is linear in v, by (5) above and (1) of 10.11. There- 
fore, by 10.138, there is a unique vector keK such that 

ge(v) = (2, ).  %) 
10.32* Directly from (1) of 10.11 and (6) above, if 7 = ak + bé, then 

gi(v) = ag.(v) + bgi(v). 
From (7), therefore 
(v, 7) = av, k) + bi», 6) 
for all veV. Comparing this with (1) of 10.11, we see that 
j = ak + 60. (8) 
The second item of (5) above then holds for vectors keK. 


That & = k follows easily: We have from (6) and (7), written for the 
vector k, that 


(@, E) = (, f). (9) 
We also have, by writing (6) and (7) for vectors and k that 
G, k) = @, B). 
Taking complex conjugates of these two numbers, and using 9 = v 
from (5), we have 
(v, k) = @, ). (10) 
Then (9) and (10), which hold for all veV, identify & and & by 10.13. 


10.34* We have now showed in (5), (8) and (10) that this complex 
conjugate satisfies the formal properties of the conjugate for n-tuples 
introduced in 7.2. 
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10.35. The abstract scalar product of 10.11 turned out in the end to be 
no more than the concrete one of 7.0 when we restrict our attention to 
n-tuples derived from vectors by the use of coordinate frames. In a 
similar way, it is not hard to show that there always exists a coordinate 
frame in which the abstract conjugation now introduced has the form 
of 7.2. This will be done in the Appendix (20.2). 


10.36* Our need for writing out the components of vectors has now 
almost vanished. Henceforth we shall use subscripts to denote particular 
vectors, e.g. v,, rather than components. 


10.4* A vector will be called real if it is equal to its own conjugate. 
A manifold will be called real if it contains with each vector also the 
conjugate of that vector. V and K are then real. A basis will be called 
real if it is made up of real vectors, and a frame will be called real if its 
bases are real. Any frame in terms of which our conjugation operation 
takes the form of 7.2 is real by definition because its basis vectors in 
that frame have components which are 0 or 1. The vector 0 is real, 
similarly. 


10.41* The basis dual to a real basis is real, for if 
(v,, Ks) = Srey 
then by (10) of 10.3 and the hypothesis that v, = 0; , we have 
(ve, hs) = dre = See 


so the k, satisfy the same equations as the k. The uniqueness of the 
basis dual to v1 , --- , v, then proves that k, = k,, 1 <s <n. 


10.42* Any vector v can be written 
v= 1 + ibe 


where v and v2 are real. Namely 


1 
n= 5 + dy, 


1 2 
v2 = 5 (v — 9). 


10.5* Itis shown in Halmos’, par. 34, that if veV, keK are represented 
by [v], [k] in some coordinate frame, and by [v]; , [k]: in some other frame, 
then there is a nonsingular matrix [W], which (a) depends only upon the: 
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two frames, and (b) relates these n-tuples as follows: 
[oh = (WY, 
[k], = [W}*[k]. 


It is easy to show that if [W] has real elements, so that [W]* = [W]’, 
then the two frames involved above are either both real, or else neither is 
real. Also, conversely, if both frames are real, then necessarily the [W] of 
(11) has real elements and [W]* = [WY]. 


10.6* Some further important geometrical notions must be mentioned 
before we proceed. 

If V; and V2 are disjoint linear manifolds in V—i.e. linear manifolds 
having in common only the single vector 0O—we write 


, (11) 


Vi © Vo 
for the linear manifold consisting of all vectors v = 1, + ve, where 
vieV;, 2 = 1, 2. The circle around the plus sign is used to denote the 


disjointness of Vi and V2. 
It is shown in Halmos’, par. 19, that if 


V=Vi © Ve (12) 
then 
K = Ki ® K, 5) (13) 


where K; = (V2)°, Ks = (V;)’ and the dimension of K; is equal to that of 
V:,2 = 1, 2. We call (13) the decomposition dual to (12). We some- 
times write K; = V# to denote the K; dual to V; in the decomposition 


(13). It is shown in Halmos’, loc. cit., that there exists a basis 11, °°: , 
vz in V and its dual ki, --- , kn in K such that, if r is the dimension 
of Vi ) 
%1,°':,0-, isa basis for Vi 
Vrgi,°**, Un isa basis for Ve 
: ; (14) 
ki,--:,k, isa basis for Ky 
Kra1, +++, kn is a basis for K,. 
Furthermore, if 0, , +++ , v, is any basis in V satisfying the first half 


of (14), its dual basis satisfies the second half, and dually. 
We shall show in the Appendix that if any one of Vi, Ve, Ki, or 
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K, is real, then they all are, and that in this case the bases (14) can be 
chosen to be real. . . 
Similar considerations apply to decompositions into more summands: 


if 
V=Vi@ V2 ®-::: © Van 
then 
K=K,@K.@::: @K,, 
where 


vi=K=Nvi=a(d vi). 
Tt jAt 


XI. GEOMETRICAL CORRESPONDENCES 


11.0 With the geometry of V and K now in hand, we consider the 
geometric aspects of our network theoretic concepts. 

The definition in Section 4 of general 2n-pole describes a concrete 
thing and stands unaltered in our geometric view. The definitions in 
6.11 of the terminology typified by ‘‘N admits [v, k] at frequency p”’ 
are unchanged except that we should now explicitly indicate that we are 
discussing concrete n-tuples of complex numbers by placing brackets 
around the vector symbols, thus: [v], [k]. In other words, a 2n-pole is 
described by a concrete relation between n-tuples. 


11.1* All of the postulates Pl, --- , P7 are stated in a language which 
now has been given an absolute geometric meaning. In this meaning, 
Pi and P2 describe a geometrical linear correspondence between vectors 
veV and keK. This is the geometric counterpart of the concrete 
notion of a linear correspondence between n-tuples. 


11.11 An impedance matrix, as in 6.3, describes a particularly tightly 
knit linear correspondence, namely a linear function from K to V. 
The geometrical counterpart is an impedance operator which for each 
» is by definition a linear homogeneous function which assigns to each 
vector keK a unique v = Z(p)keV. That is: an operator is a functional 
relatioriship between vectors and as such has a geometric identity. 


11.12 It is easy to provet that, given an impedance operator Z(p), 
and given any coordinate bases in V and K respectively, there is a 
matrix [Z(p)], with elements Z,.(p), 1 < r, s < n, such that relative to 
these bases the coordinates k, of a vector k and the coordinates v, of 
v = Z(p)k are related by (7) of 6.3. We call [Z(p)] the matrix of Z(p) 


* Technical paragraph as explained in Section 2.91. 
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relative to the given pair of bases. A strong analog of this observation is 
contained in the following lemma. 


11.138* Lemma: (i) Let LZ be a geometrical linear correspondence. Fix 
any real coordinate frame and let [Z] be the linear correspondence 
whose paired n-tuples are 


[[ol], (Al, 


where 
[v, kleL(p). 


(ii) Alternatively, let [Z] be a (concrete) linear correspondence be- 
tween n-tuples. Interpret the n-tuples related by [Z] as representing 
vectors in some real coordinate frame. Let LZ be the geometrical cor- 
respondence whose pairs, expressed as n-tuples in this frame, are those 
of the concrete correspondence [LZ]. 

In either case, (i) or (ii), the geometric correspondence L satisfies the 
geometric postulates Pl ,---, P7 if and only if the concrete corre- 
spondence [JZ] satisfies the concrete forms of these postulates. 

The proof of this lemma consists essentially in reading the postulates 
carefully. We shall not reproduce it. 


~ 11.2 Our position is now this: We have on the one hand geometrical 
objects, vectors v, k, operators Z(p), Y(p), and geometrical correspond- 
ences L. On the other hand, we have concrete n-tuples [v], [k], matrices 
[Z(p)], [Y(p)], and linear correspondences [L]. Given any pair of bases 
in V and K, in particular, given any coordinate frame, each geometric — 
object generates a corresponding concrete object which represents it 
relative to those bases or that frame. Conversely, given a concrete ob- 
ject [E], we can choose a frame in V and K and find that geometric object 
£ whose coordinates in the chosen frame are given by [é]. 


11.21* The concrete object, linear correspondence, defines a linear time- 
invariant 2n-pole by 6.21. To complete the picture, we might say that a 
geometrical correspondence L defines a Cauer class of 2n-poles. 


_11.22* This terminology is motivated by the following observation: 
if [ZL] and [LZ]; are linear correspondences representing L in two distinct 
real frames, then there exists a real nonsingular matrix [W] relating the 


([v], (el elZ](p) 
and the 
oh , (ehJe[L].(p) 


* Technical paragraph as explained in Section 2.91. 


FORMAL REALIZABILITY THEORY—I 249 


by the formulas of 10.5. This means that [LZ] and [ZL]; are related like the 
[Z] and [Zw] of 9.2. The 2n-pole associated with [Z]; therefore is a Cauer 
equivalent of that associated with [Z]. 


11.23 The observation of 11.22, combined with (ii) of 11.13, gives an 
alternative proof of 9.31. This proof is deceptively free of calculation, 
but of course the calculations are concealed in the extensive geometrical 
developments of Section 10, many of which are there offered on faith. 


XII. THE FUNDAMENTAL LEMMA ; 
12.0 This section is devoted to the statement, and the proof in part, of a 
lemma which, on the face of it, looks like an exercise in manipulating the 
postulates. In fact, the content of the lemma, and most of the details of 
its proof, are essential in what follows. To postpone them would force us 
into needless duplication of effort. 

Lemma: Let L be a geometrical linear correspondence satisfying P1, 
P2, P4, P5(I), P6(1), P7(1) and the following weak form of P3(I): 


P3’(I): If peP',, then Kz(p) = Ki D (Vu(p))”. 


Then there is a frequency domain I, C I,, differing from IT, by a 
finite set, such that L satisfies all of the postulates for pel’, 

The statement of the dual result is evident and will be omitted. 

The proof that LZ satisfies P3 will be gtven in this section. Verification 
of the remaining postulates will follow in paragraph 16.6. 

We assume that the properties of positive real (PR) functions are 
known. They are summarized for later use in Section 15. We make 
occasional advance references thereto. 

To the proof: 


12.01 First, K, is a real manifold and for pel, 
Ki € (Vin(p))’. (1) 
This, with P3’(I), gives P3(I) for L. 

Proof: Ky is real, as in 7.51. Consider now a peD, and a veVin(p); 
then [v, O]eL(p). Consider any real jeK,; then there is a weVi(p) 
such that [u, j]eL(p). Now 0 and j are real. Hence by P5(I) 

(v, j) = (u, 0) = 0. 
Therefore any real jeK, has a vanishing scalar product with every 
veVio(p). Since K, is real, it is spanned by real 7 and (1) follows. 


12.1 By the dual of 7.7, if we know that 
[v, kJeL(p), 
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then the value of (v, k) is determined by k. This makes it possible to 
state P6(I) and P7(I) for LZ (we take P6(I) to include the hypothesis 
that T, has a finite complement). 


12.11 If keK, , then J;(p) is PR. 
Proof: if k is real then 
Jip) = (v, k) = @, bi), (2) 


where, of course, [v, k]eL(p). Then however [d, k]eL(p), by P4. Hence 
by 12.1, (2) gives us 


E@ = J,(p). 


From this and P6(I), P7(I) we conclude that J:z(p) is PR for any real 
keK, . 
Now, given any keK, , we have keK, by 12.01. Then 


k= ky + tke 


where k, and ke are real and in K,, since K, is a linear manifold (see 
10.42). Let 


[vr , krJeL(p), 
r = 1, 2. Then we have (P2) 
[vy + ive, kleL(p). 
Then 
Ji(p) = (v1, ki) + (2, ke) + tor, ke) — t(v2, ha). 
Now by Pd(1), (1, ke) = (v , ki). Hence | 
Ji(p) = (v1, ki) + (v2, ke) (3) 


for any pel’, . Since each summand in (8) is a PR function, it follows 
that Jz(p) is PR for any keK. 


12.12 Let K; be the set of all vectors keK, such that 
Ji.(p) = 0 for every peI',. 


Notice that we do not assert that K, is a linear maanifolds 
If keKy then keK, and (8) above applies. Then 


(v1, Ia) + (v2, Ix) = 0 


and, using this and the PR property of each summand, we conclude 
that k, and ke are in Ky. 
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12.13 We wish now to show that K; C Ky,(p). Consider a real 7eK,, and 
a real keK,. Let 


[w(p), jleL(p), 
[v(p), kleL(p). 
‘Then, given any real A, by P2 
[u(p) + dv(p), J + MieL(p). 
Then, because keK, , . 
(u + dv, j + AK) = (u, J) + AQ, J) + ACY, Fk). 
Since 7 and k are real, by P5(I) this can be written 
(w+ 0, J + AK) = (uJ) + AO, j). (5) 


Choose any 7; such that Re(p1) > 0. Then P7(I) implies that the left 
side of (5) has a non-negative real part at p = p,. The right side, by 
suitable choice of A, can have any chosen real part unless 


Re(v(pi), j) = 0. (6) 


Hence P7(I) implies (6). Now (v(p), 7) is a rational function, by P6(I) 
applied to the other members of (5). By (6), this rational function has a 
vanishing real part throughout the right half p-plane. Hence it is an 
imaginary constant: 


(4) 


(o(p), j) = ta. (7) 
Then 
(x(p), i) = @(p), i) = —ta. (8) 


But [v(p), k]eL(p), so [v(p), k]eL(p) by P4. Since also [v(p), kleL(p), 
by 12.1, we have from (8) that 


(v(p), j) = —ta. 
Comparing this with (7) written for p, we have a = 0 and 
(v(p),j) = 90 for peD2. (9) 


Now v(p) was determined by (4) wherein k is real. For any keK, , 
k = ky + tke, where k; and /» are real and in K, (12.11). A correspond- 
ing v(p) satisfying (4) can be written 


v(p) = u(p) + wep); (10) 
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by P2, where [v,(p), k,-]eL(p), r = 1, 2. Then (9) holds for each of 
0:(p), ve(p) and therefore also for the v(p) of (10). 

We have showed now that for any pel’, and any keK,, the v(p) of 
(4) has a vanishing scalar product with every real jeK,. Since Ky, 
is spanned by real J, 


| v(p)e(K1)” = Vio. (11) 

12.14 By (11), 
[v(p), OleL(p). 
Comparing this with (4), and applying P2, 
[v(p) — v(p), k — O] = [0, k]eL(p). 
Since k is now any vector in Ky, we have 
Ki © Kn(p) C Ky (12) 

for every peI',. 


12.15 We can now also show that V.(p) C (Ki). We return to 12.13 
and read (9) thereof as originally derived for real 7 and k. Applying 
P5(I), we have from (9) that 


(u(p), k) = 0 for pel, . (13) 


By the argument immediately following (9), (13) also holds for any 
keK, , provided 7 is real. As in 12.11 any jeK, can be written 
gj = ji + te, where 7; and je are real, and the corresponding 

u(p) = w(p) + tu2(p) 


where [u,(p), jrleL(p). Therefore, finally, (13) holds for any w(p) 
satisfying (4)—i.e., any u(p)eVi(p)—and any keK, . Therefore 


Vi(p) & (K)” (14) 
for any peI',. 
12.2 We now fix our attention on a specific real poeT, 
12.21 By P4, if 
[v, k]eL(po) 
we have also 
[3, kleL(po) = L(p). 


In particular, K1o(po) is real. 
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12.22 We ean now show that K;, is a real linear manifold. Consider a 
real keKpo(po). Then [0, k]eL(po) and by 12.1 
J (Po) = 0. 


Then by 12.11 (and 15.12), Ji(p) = 0, so keK,.. Since Kio(po) is spanned 
by real k (12.21), we have 


Ki(po) © Ky. 
Comparing this with (12) gives us 
K 10(po) = K, ~ (15) 


Since Kio(po) is a real linear manifold by definition and 12.21, we see 
. that K, is. 


12.3 Let us now write, by (12) and (15), 
K, = K, (3) K, (16) 


where Kz is an arbitrary fixed manifold disjoint from K, and with it 
spanning K, . All three manifolds are real (12.21, (15), 10.6). 
Choose a K; disjoint from K, such that . 


K=K,;0K,@K,. (17) 
Let the decomposition of V dual to (17) be (10.6) 

V=V3:0@ V2 @ Vi. 
Then V; = (Kz @ K,)° = (Kz)° = Vio by 12.01. Hence 


V=Vn@V: @ Vi. (18) 
By (14) and the definitions, 
Vio € Vilp) & Vin © Ve. (19) 
12.31 Consider a real po . Then by P3’(I), (15) and (16) we have 
K1o(p0) © Ki(po) | Ke © Kro(p»). (20) 


This is now an expression dual to (19). We shall prove next that, given 
any keK1(po)MK2(= Ke), there is a unique veeVz(po)NVe2 such that 


[vz , keL(po). (21) 


Dually, given any veVi(po)NV2, there is a unique k,eK1(po)MKe such 
that 


[v, ke]eL(po). 
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The proof is a standard one in algebra and depends only. upon P2, 
(19), and (20). 
Proof: Given keK,(po)MKz , there is some veVz(po) such that 
[v, kleL(po). (22) 
By (19), then, 
Vv =U) + 
where weVio, v2eV2. Then 
[v0 , O]eL(po) 
so, applying P2 to this and (22), 
[v — »,k — O] = [ve, kleL(po). (23) 


Hence meVi(po)MV2 and v, = v2 satisfies (21). Suppose now 
v3¢V (po) V2 and 


[vs , k]eL(po). 
Then using this with (23) and P2 
[ve — v3, OleL(po). 


Hence (vz — v3)eVio. Now Vi(po)MV2 is a linear manifold and contains 
v2, v3. Hence 


(v2 — vs)eViol Vz (po) V2 = 0. 


Therefore v2 = 03. 
The dual argument completes the proof. 


12.32 The argument actually exhibited in 12.31 uses only P2 and (19), 
hence the v; of (21) is unique whether or not po is real. Indeed, this is 
true even when keK,. | 

12.383 The result of 12.31 establishes a bi-unique linear mapping between 
K, and Vz(po)MV2 . Hence these two manifolds are of the same dimen- 
sion. Since Kz and V2 = K*# are of the same dimension by construction, 
it follows that 


Vi(po)NVe => V2 
and, by (19), that 
Vi(po) = Vio @ V2. 


12.4 Let us now introduce a real frame in V and K which provides real 
bases in K, , Ke , K3 and in Vi, V2, Wz of (17) and (18). Let ki, +++ hm 
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be the basis vectors spanning K,. By 12.32, there are unique vectors 
ui(p),-°: , Um(p) in Ve such that 


[u-(p), krleL(p). 


Let, **+ , Um be the (real) basis vectors in V2 dual to the ky , +++ , km : 
| Oe KS 06 | i ee ae (24) 
Since the u,(p) are all in V2 we have for each pel, 
us(p) = Da arel(p)or (25) 
where the coefficients a;;(p) are calculated by (24) to be 
| Ars(p) = (us(p), kr). (26) 
12.41 Because the k, are real, P5(I) implies that 
Asr(p) = (ur(p), ks) = (us(p), hr) = Grs(p). (27) 


By the reasoning just following (8) and by the uniqueness of the 
us(p)eVe , Since V2 is real, we have u,(p) = u,(p). Then 


Gre(p) = (ua(p), Ix) = (tal), br) = rel). 
12.42 We have by P2 that 
[ur(p) + dus(p), br + AMs|eL(p), (28) 
for any \. The identity 
(up + Us, ky + ks) — (ur — Us, Ky — ky) 
= 2(u,, ks) + 2(us , kr) 


holds in fact for any vectors u,, us, kr, ks . Using (27), (28) and P6(I), 
it exhibits a,.(p) as a rational function. 





(29) 


12.5 Consider the m X m matrix [Zi(p)] whose elements are the a,;(p). 
the s-th column of this matrix consists of the components of u,(p). 
The rank of the matrix is by definition the dimension of the space 
spanned by these columns. 


12.51 Now the rank of [Z;(p)] can be expressed in terms of the vanish- 
ing or not of its various minor determinants. There are finitely many 
such minors and each is a rational function. Each is either identically 
zero or else vanishes at only finitely many points. Hence the rank of 
[Zi(p)], except at these finitely many points, and at the p in the comple- 
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ment of T,, is a constant. We call this constant the nominal rank of 
[Z1(p)]. . 
12.52 Let I’, consist of all pel, where [Z,(p)] has its nominal rank. 
Then I’, has a finite complement. By the reality result of 12.41, if 
pel’, then pels 

It is clear that at any pel’, the rank of [Z:(p)] does not exceed its 
nominal rank. 


12.53 By construction, the vectors w:(p), «++ , Um(p) all lie in Vi(p)NV? . 
By the reasoning of 12.33, at any real poeI', they span V2. Hence the 
nominal rank of [Z,(p)] is m. Therefore, for any pel", , (Z:(p)] has rank 
m and the w(p), «++ , Um(p), lying in V2, still span V2. Therefore for all 
pel, 
Vi(p)NV2 = Vo . 

By (19), then, 

Vi(p) = Vi © V2 = Vi, (30) 
a fixed manifold, for all pels. 


12.54 It is clear by its construction (cf. Halmos’, par. 26) that [Zi(p)] 
describes the mapping of 12.32 from Kz to Ve = Vi(p)NV2 by 


[ve] = [Z1(p) A]. 


Here the m-tuples [v,] and [k] are the components of », and k relative to 
the bases now available in V2 and K.. 


12.55 We repeat 
K; © Kw(p) © Kr = Ki © K. (12) 


Fix a pel, and a keKz9(p){Ke . Then [0, k]eL(p). Since 0eV2 , it follows 
from 12.54 that [Z:(p)] annihilates k. Suppose m + 0. Since the rank of 
[Zi(p)] is m, it follows that k = 0. Hence for pel, 


| Kio(p) NK. = 0. 
By (12), then, (31) 
| Kio(p) =K,= Ky, 


a fixed manifold. This, with the result of 12.58, proves that L satisfies 
P3(A), when m + 0. 
If m = 0 then V2. = 0, Ke = 0 and (81) follows from (12) and (16). 


12.56 [Z:(p)] is of dimension m and rank m for any pel, . Therefore 
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the correspondence of 12.32 and 12.54 between V2 and Kz is bi-unique 
for any pel, . This extends 12.31 to any pel’, . 


12.57 If m = 0, ie., if Ve = Ke = 0, then Vio = (Ki)? and the fact 
that L satisfies all the postulates is trivial because all scalar products 
(v, k) for veVz = Vio and keK, = Ky are zero. If m ¥ 0, we have yet 
to show that L satisfies P5(A), P6(A), P7(A). 


12.6 Since now L satisfies P3, 7.7 as given is applicable and we find 
(with 12.1) that if pels and 


[v, k]eL(p), 
then (v, k) is fixed by either v or k. Furthermore, 
(vk) = (v +m, k + ko) 
for any veVio , KoeK x . 


12.61 If pel, and [v, kleL(p), then veV,, keK,. By (30), (31), and 
(16), therefore, there exist veVio, koe€Kio such that u = v — weVe, 
yi _ (k = ko) Ke . Then by P2 


[u, jleL(p). (32) 


By 12.6, then, any value assumed by a scalar product (v, k) with 
[v, k]eL(p) is also assumed ‘by a product (u, 7), where (32) holds and 
weV> : 7eKy ‘ 


XIII. SUFFICIENCY OF THE POSTULATES 


13.0 We suppose that L satisfies the postulates of 12.0. Then the results 
of Section 12 are applicable. The ones of first importance are contained 
in the facts from (15), (80) and (31), that 

Vi = Vw © V2, 

K, = K, ® Kio ) 


where the choice of K. was governed only by the requirement that the 
second of these formulae hold. 


13.01 Considering K, and V2 as separate spaces, V2 = K by 10.6. 
Let M be the geometrical linear correspondence between them with 
frequency domain I, and pairs described by 12.31 and 12.56 (or 12.54). 
That is, as vectors in V2 and Ky 


[v, kel (p) 
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if and only if, as vectors in V and K, 
[v, k]eL(p). 
13.02 In the real frame of 12.4 let us renumber the basis vectors so that 
V1, °°" 4 Ur span Vio, 
Ur4t,°°* » Urtm Span Ve, 
Vite 5G Un span Vi ; 
Then 
ky, +++, Kram span Ke, | 
Krtm4i,**t,kn span Ky. 
We say that such a frame reduces L. 


13.1 Let us now interpret the s-th components of [v] and [k] in this frame 
respectively as the voltage across and the current in an ideal branch 
B, of a2n-poleN, 1 <s <n. 

By construction, the vectors veV, in this frame have components 
Grim = *'' = a, = 0, since, -*: , %4+mSpan V,. At the same time, 
the components b;4.m41,°°* , bn of [k] may be chosen arbitrarily without 
altering the fact that [[v], [k]]e€{Z](p) because of 12.06. Therefore, the 
ideal branches 6,4m41, °°* , Bn can each be realized physically by a short 
circuit. 

In a dual way, since ky41, +++ ,kaspan K,, any keK, has components 
bi, °°: , 0, all zero in our chosen frame. Furthermore, the components 
a1, °°, a, Of [v] can be chosen at will. Hence the ideal branches 
Bi, °°* , Br can each be realized physically by an open circuit. 

Let N, now be the 2m-pole whose ideal branches are B41, °°* 5 Brim 
Let the pairs [[v], [k]] admitted by Ni at each pel, be the [fv], [A], 
where [v, k]eM(p) (13.01). The representation just found for N shows 
that N is physically realizable if and only if N, is. 


13.11 The matrix [Z,(p)] of 12.54 is the impedance matrix of the 2m- 
pole Ni. 


13.12 We now show that [Z;(p)] is a positive real matrix. The displayed 
formulae of 12.41 show (ii) and (iii) of 1.1, and 12.42 shows (i). Now 
suppose that [v, k]leM(p). Then, as vectors in V and K, [v, kleL(p) by 
definition of M(p). Then, however, if k is fixed 


Ji(p) = (v, k) 
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is a PR function (12.11). Regarding v and k in V2 and Ke let 
[Dra nr) Dr+n] = [k]. 
Then by (1) of 7.0 


m 


mcr k) = 2d oD det(p) DezrDere 


=1 
and this has a non-negative real part if Re(p) > 0. This is (iv) of 1.1. 


13.2 We can now prove the lemmas 8.1 and 8.2. Given a linear cor- 
respondence [L] which satisfies Pl, --- , P7 by 11.13 we can interpret 
[ZL] as the concrete correspondence representing a geometrical cor- 
respondence L in some chosen real frame, and L satisfies P1, --- , P7. 
Then by the results in 13.01-13.12 there exists a real frame in which the 
representative [ZL], of Z has the properties claimed in 8.1 and 8.2 for Dy . 
But we saw in 11.22 that [Z] and [Z]; are related by a real matrix W like 
the Z and Ly of Section 8. Q.E.D. 


13.21 With the proofs of 8.1 and 8.2 we have reduced the sufficiency 


claimed for Pl, --- , P7 in 8.0 to the sufficiency of positive reality of 
[Z(p)] claimed in 1.1, by the argument outlined in 8.5. 


XIV. OPERATOR-VALUED FUNCTIONS OF p 


The next three sections are directed principally toward the proof of the 
matrix theorem of 1.1. They do however, contribute to 12.10 and to 
the necessity proof. 


14.0 We continue to use the geometric language. The reader who re- 
gards this as unduly pedantic is free to place a concrete interpretation 
upon every argument, for all of the arguments are either frankly based 
on matrix representations or upon the three identities: 


14.01 (Zj, k) = (Z*k, 7) for all j, keK. 

1402 Z& = (Zk) for all keK. 

14.03 Z’ = (Z)* = (Z*) 

14.04 These identities are obvious for matrices using 7.0 and 7.2. 
Geometrically, the first and second define, Z* and Z, and the third 


defines Z’ in two ways. The equivalence of these two ways is a theorem 
based on (10) of 10.33. 


14.05 The symbol Z will always denote an impedance (operator, matrix, 
scalar), and Y will always denote an admittance. An impedance oper- 





260 THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1952 


ates from K to V, an admittance dually. The operators in Halmos” 
are physically dimensionless, in that they operate, e.g., from V to V. 
This difference is scarcely noticeable. 

We shall regularly omit the duals to concepts or proofs given in terms 
of impedances. In doing so, we adopt the rule that the dual to an 
expression 


(Zk, k) 


(v, Yv). 


14.1 An operator is called symmetric if Z = Z’. Such operators have 
three useful special properties: 


14.11 If Z is symmetric and j and & are real, then 
by (10) of 10.33, 14.02, 14.01, 14.03, and hypothesis. 


14.12 Let k = ky + tke , where k; and kz are real (10.42). If Z is symmetric 
then 


(Zk, k) = (Zky ) ky) -+ (Zke ’ ke), 
for, by 14.11, 
(Zk, ; tke) = —1(Zky ; ke) —1(Zhe ) ky) 


at (Z (tke), ky). 


lI 


(Cf. the similar identity in 12.11.) 


14.13 The symmetric operator Z is completely defined by the quadratic 
form 


(Zk, k) (1) 


as a function of real keK. For 14.11 permits the formula (29) of 12.42 
in any real frame, where u, = Zk, . The matrix elements of [Z(p)] in 
that frame are then defined by that formula in terms of values of (1) 
for real k. 

The form (1) specifies any Z (symmetric or not) if k is allowed to 
range over all of K (Halmos’, par. 53). 


14.2 Let Z(p) now be an impedance operator depending on p. We say 
that po ~ © is a pole of order m of Z(p) if 


tk) = lim (p — po)"(Z(p)k, ke) (2) 
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exists for every keK and is not identically zero. By 15.13, this limit 
£(k) defines an operator Ry , the residue* of Z(p) at po , by 


(Rok, k) = €(k) for keK. 


The changes in (2) required to define a pole at p = © are obvious. 
14.21 A pole po of order m of Z(p) is a pole of some matrix element of 
[Z(p)], of order m, in any frame, and no element of [Z(p)] has a pole at 
po of order exceeding m. For the elements of [Z(p)] are defined by the 
values of (Z(p)k, k), by 14.11 and Halmos? loc. cit. 


XV. POSITIVE REAL FUNCTIONS 
15.0 Let f(p) be a scalar function of the complex variable p. Following 
Brune’ we define f(p) to be positive real if 


(i) f(p) is a rational function of p, 


(ii) f(p) = 1), 

(iii) Re(p) > 0 implies Re(f(p)) > 0. 

The property (i) of being rational is of course on a quite different 
level of ideas from the other properties, but it saves words later to in- 


clude it specifically in the meaning of positive real. 
We abbreviate the words positive real to PR. 


15.01 The open region of the complex plane consisting of all finite p 
such that Re(p) > O—the right half plane—we denote by Tx. . 


15.1 Brune, loc. cit., established a number of properties of PR functions 
f(p) which will be useful to us here: 


15.11 f(p) has no poles in Ty. 
15.12 If Re(f(p)) = for some pel’, , then f(p) = 0 for all p. 


15.13 If it exists, ——~ is PR. 
i) 


15.14 If f(p) has a pole at p = po, it has one at p = po. 
15.15 If f(p) has a pole at p = ia» , that pole is simple and 


2p 
fp) = ota + fap), 


where 7 > 0, and fi(p) is PR. 


* Properly, Ro is a residue only when m = 1. There is no convenient name 
available for general m. 
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15.16 If f(p) has a pole at p = ~, that pole is simple and 


S(p) = pr + filp), 
where r > 0, and fi(p) is PR. 


15.17 We shall use all of these in the next section, save 15.13. Our aim 
is to prove properties analogous to 15.11, --- , 15.16 for PR matrices 
and operators. 

The reader familiar with the Brune process’ for realization of a 2-pole 
will remember the importance of the properties 15.11, --- , 15.16 for the 
success of that process. Correspondingly, we must establish the analogs 
of these properties to implement the general Brune process for 2n-poles. 


XVI. POSITIVE REAL OPERATORS 


16.0 An operator Z(p) from K to V will be called positive real (PR) if in 
some real coordinate frame the matrix [Z(p)] is a PR matrix in the sense 
of 1.1—that is 


(i) [Z(p)] has rational elements Z,;(p) 


(ii) Zpa(p) = Zro(B) 
(ili) Zre(p) = Zer(p) 
(iv) For any real keK and any pel; 
Re(Z(p)k, k) > 0. 
We intend in this section to establish for PR operators the properties 
listed below. By subtracting 0.9 from the designation of each. property 


one obtains the designation of the analogous property of a PR scalar 
function, stated earlier. 


16.01 Z(p) has no poles in T,. 

16.02 If Re(Z(p)k, k) = 0 for some peT'; , then Z(p)k = 0 for all p. 
16.03 If it exists, Z(p) = Y(p) is PR. 

16.04 If Z(p) has a pole at p = po , it has one at p = fin. 

16.05 If Z(p) has a pole at p = tw , that pole is simple* and 


2p 
oe rea i Zi(p), 


where PF is real, symmetric, and semi-definite, not zero, and Z;(p) is PR. 


Z(p) = 


*i.e., of order one. 
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16.06 If Z(p) has a pole at p = , that pole is simple and 
Z(p) = pk + Z(p) 
where Rk = R' = R, R > Oand Z,(p) is PR. 


16.07 There is property of rational scalar functions f(p), whether PR 
or not, that is essential in the Brune theory: the existence of a finite 
integer, the degree of f. Each step in the Brune reduction of f(p) leaves 
an unreduced portion which is of lower degree than the function upon 
which the step was performed. The finiteness of the original degree of f 
then guarantees the termination of the process in finitely many steps. 

There exists also for rational matrices (and operators) a concept of 
degree. This degree plays the same role in the general Brune process for 
2n-poles as the degree of a scalar function does in the process for 2-poles. 
To define this degree and develop its properties requires an excursion 
into classical algebra. Since we shall not need these ideas until Part IT 
we defer further discussion of them to that part. 


16.1 If Z(p) is PR it follows at once that the matrix [Z(p)] is PR in any 
real frame. 
_ Proof: Two such matrices are related by 


[Zh = (UZ) [UY 


where U is real, by 11.22 and the argument in 8.6. The PR properties 
of [Z(p)] are obviously preserved by this operation. 


16.11 If Z(p) is PR, then 
Z(p) = Z'(p) = Z*(p) = Z(p). 
Proof: Use 16.0 and 14.03 in a real frame. 
16.12 If Z(p) is PR, then for any given keK the function 
Ji(p) = (Z(p)k, k) 


is a PR scalar function. It follows that the limitation in (iv) of 16.0 
to real k is a simplification, not a restriction. 

Proof: J;,(p) is independent of coordinate representation. By use of a 
real frame, (i) of 16.0 implies (i) of 15.0. 

By 14.01 and 16.11 


Jip) = (Z*(p)k, k) = (2(p)k, k) = Je(p). 
This is (ii) of 15.0. For any k, 14.12 and (iv) of 16.0 imply (ii) of 15.0. 


264 THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1952 


16.13 Conversely to 16.12, if Z(p) is symmetric and J;(p) is PR for 
every real k, then Z(p) is PR, and J;(p) is PR for all k. 

Proof: J(p) is rational so (i) of 16.0 holds in any frame by 14.13. 
Clearly (iv) of 16.0 holds. 

Now for real k, by (10) of 10.33 and 14.02 


Jp) = Inlp) = (Z(p)k, ki). 


Hence Z(p) = Z(p) by 14.13. This is (ii) of 16.0, and (iii) there holds 
by hypothesis. 


16.2 Proof of 16.01: By 15.11 and 16.12, J;(p) has no poles in T,. 
This is 16.01 by the definition 14.3 of pole. 


16.21 Corollary: Any PR Z(p) can be considered as defined throughout 
Ty: for any k, J.(p) is defined throughout T, by 16.2. For each p, 
as a function of k, J:.(p) defines Z(p) (14.138). 

16.3 Proof of 16.03: In any frame [Z'(p)] = [Z(p)|” = [Y(p)] consists 
of rational elements, by direct calculation of the inverse matrix. In a 
real frame [Y(p)] = [Z~'(p)] is symmetric and real for real p by the same 
argument (both facts are also deducible geometrically). Hence we have 
the duals of (i), (ii) and an) of 16.0 for Y(p). Clearly Y(p) is defined 
throughout IT, . 

Now suppose that for some veV and some pel, we have 


Re(v, ¥ (pov) < UV. 
Then there is a keK such that v = Z(po)k. Therefore 
Re(Z(po)k, k) = Re(Z(po)k, k) < 0. 


Since this is impossible, we have the dual of (iv) of 16.0 for Y(p) and 
Y(p) is PR. 
16.4 Proof of 16.04: This is immediate from 15.14, 14.3, and 16.12. 


16.5 Proofs of 16.05 and 16.06: Suppose Z(p) has a pole at p = ta. 
Then (Z(p)k, k) does and that pole is simple by 15.15 and 16.12. Then 
_ by 14.3 we can write 





Dp es a it 
P — Wo 


where Z)(p) is regular at p = tw. Now Z(p) has a pole at p = —ta 
by 16.5, so a similar argument gives 


Ri + “Z(p), (1) 





1 1 
DONS 5 ag a, 
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where Z:(p) has no pole at iwy or —iw) . The symmetry of Z and linear 
independence of the terms above then imply the symmetry of Ry, R; 
and Z,(p). 

For any keK, now, 


: (Rok, k) -}- ! 


p— Lar 8) 10 








(Z(p)k, k) = (Rik, k) + (Zi(p)k, k). 


Applying 16.12 and 15.15, 

(Rok, k) = (Rik, k) > 0 
for all k. Hence Ryo = Ri = R (say) and A is semi-definite. Also, 
(Zi(p)k, k) appears as the residue fi(p) in 15.15 and is therefore PR. 
Then Z,(p) is PR by 16.18. With Ry and R, identified, (1) above is the 


expansion given in 16.05. We have now proved all of 16.05 save the 
reality of R. But 


2p 
p +a t 
is PR, by 16.18, hence is real for real p. Therefore FR is real. 
The proof of 16.06 is similar. 





16.6 To prove 16.02 we appear to digress somewhat, by first com- 
pleting the proof of the fundamental lemma of 12.0. It was established 
in Section 13 that the matrix [Z,(p)] describing M(p) in the chosen 
basis is PR. The case in which it is nonsingular (i.e., m # 0, ef. 12.56, 
12.57) remains to be examined. 


16.61 If [Z:(p)] is nonsingular then its inverse is PR (16.38). Then for 
any veV2, 


(v, k) = (v, Y@)v) | (2) 


is PR (16.12 dual). By 12.61, for any weV,, the values of the function 
F,(p) are the values of (2) for some veV2. Hence F.,(p) is PR. This is 
P6(A) and P7(A) for L. 

16.62 To settle P5 for L in 12.0, consider pel, and 


[v, kleL(p), — [u, JleL(p), 
where wu and v are real. Then, say, 
v=eutnu, 
where eV 19 , n€V2. But then 


v=D=H + YY 
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and, because Vio and V2 are real, i) = vo , d1 = v4, , and these vectors are 
real. Using similar reasoning for wu, 


0,7) = @,V¥(p)m), (wu, k) = (um, Y(p)u), (3) 


by 12.61. The equality (u, k) = (v,7) now follows from (8) and the duals 
of 16.11, 14.11. Hence we have P5(A) for Z and 12.0 is proved. 


16.7 We now prove an important 

Lemma: Let Z(p) be a PR operator from K to V. Let I, be the set 
of » where Z(p) is defined and has a rank equal to its nominal rank. Let 
L be the correspondence with domain I, and pairs 


‘([Z(p)k, F), keK, . 


Then L satisfies Pl, --- , P7. 
Proof: L satisfies Pl and P2 (6.8). T', satisfies P4 by the argument of 
12.52. Then L satisfies P4, for by 16.11 


Z(p)k = Z(p)k. 


L satisfies P5(I) by 14.11 and 16.11. Tz, satisfies P6 by 12.52. Then ZL 
satisfies P6(I) and P7(I1) by 16.12. The fundamental lemma, 12.0, 
now proves that L satisfies all the postulates. 


16.71 We call a correspondence satisfying all the postulates PR. 


16.72 Proof of 16.02: Suppose Re(Z(po)k, k) = O for some poeD,. 
Because this function of p is PR (16.12) we have 


Ji(p) = (Z(p)k, k) = 0. 
Hence keK, = Ky (12.12, 12.55). Hence [0, k]eL(p) for every peDz. 
That is 
: Z(p)k = 0 for peDz. 
16.73 Corollary: If Z(po)k = O for some poeDs, then Z(p)k = 0. 


For the hypothesis here implies that of 16.72. This is the analog of 
15.12; the result of 16.02 is stronger. 


16.8 An important consequence of 16.7 is the 

Lemma*: If Z(p) is PR and of rank m, then there exists a real coordi- 
nate frame in which the matrix [Z(p)] is an m X m nonsingular PR 
matrix [Zi(p)] bordered by zeros. 


* Proved by Cauer®. 
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Proof: Consider the PR correspondence L defined by Z(p). Then 
Vw = 0, because Z(p)0 = 0 for every pel’, . Consider the real frame 
of 13.02. [Z(p)] in this frame takes any of ky4mii,°+°* , kn into 0 because 
these span K,»). Within K2 , [Z(p)] must describe the same operation as 
the [Z,(p)] of 12.54. Because [Z(p)] is symmetric the lemma follows. 


XVII. THE JUXTAPOSITION OF CORRESPONDENCES 


17.0 This section and the next will consider ways of constructing new 
correspondences from old. This will provide the basis of the necessity 
proof of Section 19. 

17.01 It is obvious that if two physical networks are set side by side 
and their accessible terminals regarded as the terminals of a single 
larger network, that enlarged network is again a physical network. 
This is the gist of the present section. 


17.1 Suppose that 
V=V,@0 V2, K=K, @ Ky, 


where K; = Vj and all spaces are real (10.6). Let #, project on V 
along V2. (Halmos’, par. 33) and FE, = 1 — EF, project on V2 along V; . 
Then Ej} projects on K; along K;, 7 * i (Halmos’, loc. cit.). It is 
easily verified that Z; = E;, E? = Ej , from the analog of 14.02 for 
dimensionless operators. 

Considering V; and K; as separate spaces, let L; be a geometrical 
linear correspondence between them with frequency domain I; ,z = 1, 2. 

Consider the correspondence ZL between V and K defined by 


(i) The frequency domain Tr, = TiN. 
(ii) [v, kleL(p) if and only if (Ew, E*keLi(p), 7 = 1, 2. 
In (ii), of course, we regard Hw and Ek as elements of V;, K;. 
17.11 L so defined is called the juxtaposition of Z, and Ip. 
17.2 Lemma: L is PR if and only if each of LZ, and Lz is PR. 


17.21 Proof of “if”: It is clear that ZL satisfies Pl and P2. Further 
notation is now simplified if we put Li = M, ZL. = N. Consider the 
manifolds 


Vi © Vy, Virco © Vuo, Ky @ Ky, K wo ® Kyo, 


where Vy, C V; is the manifold of voltages admitted by I, = M con- 
sidered, as a correspondence between V; and K; , and Vo the manifold 
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of voltages veVi such that [v, O]eZi(p) for all pel. Dual definitions 
for Kay, Kano, and symmetrical ones for Vy, -:-, Kero need not be 
repeated. 

It is clear from these definitions that the four manifolds above are, 
in the order listed, the manifolds 


Vi, Vo, Ki, Kio 
for L. Now, for example, 
(Kio)? = Kao © Kyo)” = (K 100) (Kao)? 


by 10.6. This last manifold, in V, is (Va ® V2) N (Vw @ Vi), byP3 
for M and N; and by 10.6. But by direct calculation 


(Var ® V2)N(Vw @ Vi) = Vu © Vw = Vi. 


The dual of this result then completes P3 for L. 
P4 for L is immediate because the FE; and Ej are real. 
The duality of the decompositions of V and K implies the identity 


(v, k) = (Ew, Efk) + (Ew, Eyk) 


(that is iL, = HE, = 0, and dually. This is Halmos’, par. 33). All 
of P5, P6, and P7 for L follow at once from this identity. 


17.22 The “only if” of 17.2 is a special case of the result of Section 18. 
Tes nroof will he deferred ta 18 4 . 


17.23 It is obvious that the notion of juxtaposition and the lemma of 
17.2 extend to juxtapositions of more than two correspondences. 


17.8 Even without the ‘‘only if” part of 17.2, we have enough for the 
following characterization of PR correspondences: 

Theorem: A correspondence L is PR if and only if it is the juxtaposi- 
tion of 

(i) a correspondence defined by a nonsingular PR matrix between 
a Vi and aK, = Vr, 

(ii) a correspondence consisting of short: circuits: that is of pairs 
(0, k] for all keKe and all p, 

(iii) a correspondence consisting of open | ‘circuits: that is, of pairs 
[v, 0] for all veV3 and all p. 

Proof: If Lis PR, the decomposition indicated is that of 13.1, 13.11, 
13.12. If Z is the juxtaposition indicated, then it is PR by 16.6 and the 
“if” in 17.1, provided the short and open circuits are PR correspondences. 
The verification of the postulates for these latter is easy and will be 
omitted. 3 
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17.31 The labor of considering PR correspondences instead of matrices 
has yielded the disappointingly simple result of 17.8. We have already 
been warned of this, however, by our knowledge of the properties of 
physical networks (2.9). 


XVIII. THE OPERATION OF RESTRICTION 


18.0 In addition to juxtaposition, which is an operation on correspond- 
ences clearly motivated by physical considerations, there is an operation, 
here called restriction, which has important use in the next section. 
There the physical meaning of the operation will become clear. 


18.1 Let Vand K = V* be a pair of dual spaces. Let U and J = U* be 
another pair. Suppose that C is a given fixed linear operation from J to 
K: given any jeJ, there is a unique k(j)eK, written 
kQ@) = Cj, 
such that if k, = Cj,, r = 1, 2, then 
ayky + ake = C (aiji =r ajo) 
for any complex scalars a; , dz. 


18.11 Let (v, k) denote the scalar product between Vv and K, and 
(u, j)2 that between U and J. Given C, and any veV, let us Guid that 
unique vector u(v)eU for which 


(u(v), 2 = (, Ci) (1) 


for every jeJ. That such a vector u(v) exists and is unique follows from 
10.13 when we notice that the right-hand side of (1) defines a function 
conjugate linear in 7. Now for fixed j, the right-hand side of (1) is linear 
in v, hence so also is the left side. That is, there is a linear operation C’* 
from V to U such that 


uv) = 
The following chart illustrates the situation: 


VK 
CF] {Cc 
U j 
18.12 We suppose now that C takes real j into real k, i.e., that C is 
real. Then by (1) 


(Cv, je = (C*v, je = (v, Ci = (8, Chr. 
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By comparison with (1), we have 
C*y = CB, 
Hence C* also takes real vectors into real vectors and is real. 
18.2 Now let LZ be a PR correspondence between V and K. We define 
one, say M, between U and J, as follows: For each peI',, let M (P) 


consist of all pairs 


— [u, J] 


such that u = C*v and 
[v, CjleL(p). 
This definition can be illustrated by enlarging the chart of 18.11: 


L 
V— K: 


ct | {Cc 
U Vi J 
The w’s corresponding to jeJ can be constructed by going around 
through C, Z and C*. This then defines a direct mapping from J to U. 


18.21 We call the M defined by 18.2 a restriction of L, since its pairs 


are images under C* and C™ (which is not defined over all of K) of a 
restricted set of pairs drawn from L. 


18.22 Clearly there is a dual operation defined by an operator D from 
U to V. We might distinguish the operation of 18.2 by calling it a 
current restriction, its dual by calling it a voltage restriction. 
18.23 The restriction M of LD is defined by lists M(p) which exist for 
any pel',. The frequency domain of M has not yet been specified, 
however. 
18.3 Theorem: If L is PR, then there is a frequency domain I» for M 
such that M is PR. 

Proof: P1 and P2 for M are evident at once, for any peC',. The 
remainder of the proof is divided among 18.31, --- , 18.37 below. 


18.31 For P3, let Ja be all zeJ such that CyeK,. Then, given jeJx, , 
for each pel’, there is a v such that 
[v, CileL(p), 
whence 
[C%v, jleM (p). 


FORMAL REALIZABILITY THEORY—I 271 


Therefore J a(p), the space of currents admitted by M at frequency 7p, 
coincides with the fixed J» at each peI', . 
Clearly J ar is a real linear manifold. 


18.382 Consider now U so(p): if [u, OJeM(p), then there is a v such that 
u = C*y and 


[v, CO] = [», O]eL(p). 
Hence veVio(p) = -Vio for each peF',. Therefore, for each pe, , 
U so(p) CC*Vy. (2) 


Now suppose, conversely, that pe!', and veVio = Vuo(p). Then 
[v, O]eL(p). Now 0 = CO, so [v, CO]eZ(p). Hence [C*v, O]eM(p), so 
C*veU ao(p). This proves the inequality opposite to that of (2), so for pel’, 


U so(p) = C*Vip = Uso, (3) 
a fixed space. 
18.33 Now consider (U yo)’. If je(U so), then 
(u, i). = 0 
for every weU yo . That is, by (8), 
(C*v, je = (v, Cj) = 0 


for every veVzo. Therefore Cje(Viv)’ = Kz, and jeJ uw by 18.31. That 
is, we have proved 


Ju 2D (Usw)’, 
and, combining 18.31 with this and (8), 
Ju(p) =JuD (U s0(p))° a (U a)’. (4) 


This is the weak form P3’(I) of 12.0 for M. It is as far as we can go 
with P3 at the moment. 


18.34 Consider P4. If for peI', we have 
[u, jleM! (p) 


then [v, CjJeL(p) and wu = C*v. But then [d, CjJeL(p) and a = C*o, by 
18.12. Then however 


[a, jl eM (p) 
by definition of M. This is P4. 
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18.35 Consider P5(1): if 


[ur , jrleM (p), 
where j, is real, r = 1, 2, then 
(u, ? jaya ma (C*, ’ je)a os (v ? Cjs)s ’ (5) 


where [v,, Cj,JeL(p). Since Cj, is real | 

(1, Chr = (2, Chi) 
by Pd (1) for LZ. This with (5) for r ¥ s proves P5(I) for M. 
18.36 Fix a jeJ » and for each pel’, a u(p) such that 


[u(p), jleM (p). 
Then u(p) = C*v(p) and 


[v(p), CyleL(p), 


for some v(p). Then as in (5) above 


(u(p), j)2 = (v(p), Cj) . 

P6(I) and P7(I) for Z then imply that P6(I) and P7(I) hold for M, 
using I’, for Ty in P6. 
18.37 We now have M satisfying the hypotheses of 12.0. Therefore 
there is a Ty such that M satisfies all the postulates. This is 18.3. 
18.4 Proof of “only if” in 17.2: Suppose that LZ between V and K is 
the juxtaposition of LZ; between V; and K,, Ll. between V2 and K,. 
Let, say, U = Viand J = K,. Let C be the identity map from K, to K: 
if je} = K,, then C7 is just 7 considered as a vector in K. Then C is 
real. It is easily computed that C* is F,. 

Consider the restriction M of L based on this C. Its pairs for 
pel wCI, are all the pairs [u, 7] such that 7 = H*7eK, and u = Lp, 
where 


[v, jleL(p). (6) 
But then 
[u, j] = [Ev, £*)] 


and this is in Li(p) by (6) and the definition of juxtaposition. Therefore 
the list M(p) is contained in Iy(p). 
Suppose that [u, jJeZi(p). We have [0, OJeZ2(p) so by P2 and the defi- 
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nition of juxtaposition 
[w, gleL(p). 
But then 7 = E*j, u = Eu, and by definition of 


[u, jleM(p). 


Therefore for every pel x, M(p) = Ini(p). Therefore there is a fre- 
quency domain (IT ,,) for ZL; such that Ly; is PR. 


XIX THE NECESSITY PROOF 


19.0 Fortunately for this section, those parts of network theory which 
we require have recently been very succinctly stated by J. L. Synge”. 
We shall paraphrase them here, referring the reader to the source” for 
details of definition. 


19.01 First, we observe that in Cauer’s definition’, which we shall 
repeat in detail below, an ideal transformer with m windings is a 2m-pole 
whose terminal pairs are the termini of the respective windings. 

A system of m coupled coils is a 2m-pole with similarly defined terminal 
pairs. 


19.02 Given a 2n-pole N which is a finite passive network, let us adjoin 
ideal transformers as in Figure 1(b). We then draw the ideal graph of 
this network. Adjoin to the graph ideal generator branches y1, :-- , 
Yn 5 Yr between 7’, and 7’ ,l1<r<_n. Let 8, be the ideal branch repre- 
senting the transformer winding between 7, and 7,,1 <r < n. Enu- 
merate the remaining branches of the graph Bnii, --- , Bo. 


19.03 The branch 7, is in a mesh with 8, and no other branches. Let us 
call this the r-th external mesh. Any basic set of meshes must include 
each of these. 


19.04 Let &, ---, ¢, be the currents in the generator branches, 
ky, +++ , ky the currents in the branches 8; , --- , 6) and 
[¢] = [fs , BO Sany fa, kn, ee » ke), [k] = [kK , ee » Ky). 
Let w,, --+ , Wa be the voltages across the generator branches, 
V1,°°* , % the currents in the 61, --: , 6, and 
[w] = [wi, +--+, Wn, M1, °° , Vl, [v] = [ro , +++, Vl. 
19.05 Let us choose a basic set of meshes, let 71, --- , 7s be the respec- 


tive mesh currents, and 


[J] i [ji ; ae , Jal. 
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Let 


[u] < [uw , oe 5 Us| 


be the s-tuple of mesh voltages. We suppose that ji, --- jn, U1, °°° 5 Un 
refer respectively to the n external meshes. (Cf. 19.03.) 


19.06 The results of Synge” can now be stated as follows: 
There exists a real constant matrix [C,] of s columns and b + n rows 
(having, in fact, elements which are +1, —1, or 0) such that for any [7] 


(¢4] = [Cll] (1) 
is a set of branch currents satisfying Kirchoff’s node law, and for any [w] 
[u] = [Ci]’[w] (2) 


is a set of mesh voltages satisfying Kirchoff’s mesh law. Furthermore, 
given any [£] which satisfies the node law, there is a [j] such that (1) 
holds. 


19.07 If we interpret the [f], [yj], etc., as representations in real bases 
then [C,] is real and [Ci]’ = [C,]*. 


19.08 The matrix [Ci] has the form 


C2} 0 
(CG) Se 
0; Cc 


where [C2] is an n X n diagonal matrix (having diagonal elements +1, 
in fact). 

Proof: By construction, j:, -+- , jn are mesh currents in the external 
meshes. These are then equal, save for sign, to the currents fi, +: , ¢n 
in the generator branches. 


19.09 By 19.08, (1), and the definitions in 19.04, 


[A] = [Cll], [uw] = [C) to, 
and by 19.07, [C]/ = [C]*. 


19.1 Let us suppose that we have enumerated the branches Bn41, 
-++ , yin 19.02 in such a way that Bn41,--+ , 8. are all the two poles in the 
graph, Bci1, °°: , 8a are all the branches containing coils which are 
magnetically coupled, and Bay1, --* , 8 the remaining ideal branches of 
ideal transformers. 

Let [Za(p)] be the (d — n) X (d — n) impedance matrix relating the 
voltages across the branzhes 8,41, ++: , 8a to the currents in them when 
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we consider the individual two-poles and the system of coupled coils 
as separate unconnected networks. Then [Za(p)] is composed of a 
(c — n) X (c — n) diagonal matrix in the upper left field and a (d — c)-X 
(d — c) matrix in the lower right, with zeros elsewhere. 


19.11 The diagonal part of [Za(p)] has elements drawn from the follow- 
ing list: 

(i) f(p) = p 

(ii) f(p) = dp 
(ili) f(p) = Ap 
where p, 6, \ are non-negative constants, possibly different for each 
branch. 


19.12 It is shown in texts on electromagnetic theory that the matrix 
representing a system of coupled coils is of the form 


pl), 
where [G] is a real, constant, symmetric, and semi-definite matrix. 


The lower right field of [Ze(p)] is then such a-matrix. 


19.13 It is obvious from this description that [Za(p)] is PR. It therefore 
describes a PR correspondence between (d — n)-tuples of current and 
voltage. 


19.2 We must at last consider ideal transformers in detail. Let Vi and 
K, be m-dimensional spaces represented as aggregates of m-tuples. 


Let pi, p2, °° * ; Pm be m real numbers. Let Vz consist of all m-tuples 
[a] = [a , --: , GmleVi such that 
ee a he ee 
P1 p2 Pm 


We interpret these relations as follows: 
(a) If any p, = 0, then a, = 0 
(b) If any two p,, ps are not zero, then 
Ge _ as 
Pr Ps 


(c) If only one p, + 0, then a, is arbitrary. 
Let Ky consist of all m-tuples [b] = [bi , --- , bmleKi such that 


pibi + probe ae Sees Pndin = 0. 


V,and K, are linear manifolds. 
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Let [Zr] be the concrete linear correspondence defined by the list 
[L7](p) which consists for each complex p of all pairs [[a], [b]] where 
[aleVz, [b]eKr. 

The correspondence described: by (L7] is what Cauer’ defines as an 
ideal transformer. He shows, loc. cit., how it can be defined as the 
limiting case of a physical transformer. 

There is also a dual kind of device, described by a correspondence 
admitting all [b]eK, for which 


and all [a]JeVi for which 
MG -+ e's + Asalloe = (0. 
This also is an ideal transformer obtainable as a limiting case of a 
physical one. 


19.21 The correspondence Lz is PR. 
Proof: We observe that Vr = (K,)’, for let [a]eVr, [b]leKr, and let 
t be the common value of the a,/p, . Then 





(a, b) = Za,b, = tZp,b, = t(Zp,b,) = 0. 


The postulates are now all easily proved. We omit the details. 


19.8 Let V and K be b-dimensional spaces. We interpret the b-tuples 
[v] and [Kk] of 19.04 as representing vectors veV, keK in a real frame. 
Let L be the correspondence between V and K formed by juxtaposing 
(i) the correspondence described by [Za(p)] relating components with 
indices in the range n + 1 to d, 
(ii) the several correspondences described by ideal transformers, 
relating components with indices in the ranges 1 to n and d + 1 to b. 
L is PR because it is the juxtaposition of PR correspondences. 


19.31 Let U and J be s — n-dimensional spaces. We interpret the [uw] 
and [7] of 19.04 as representing weU, jeJ in a real frame. 


19.382 Let C be the operation from J to K whose matrix in our chosen 
frames is [C]. Then C* operates from V to U with the matrix 

[C]* = [C]’. By these definitions, C is real. Let M be the correspond- 
ence between U and J obtained by restricting L with C. Then there is a 
frequency domain I’ y such that M is PR (18.3). - 


19.4 By 19.09, [/] in our chosen frame is the correspondence estab- 
lished between mesh currents and mesh voltages by the network of the 
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2n-pole N. When this network operates as a 2n-pole, the only mesh 
voltages which are not zero are those relating to the external meshes, 
since there are no internal sources of voltage. We must now account 
for this. 


19.41 Let V2, K2 be n-dimensional spaces. Choose a real frame and let 
D be the operation which takes 


[a1, °-+ , GnleVe (3) 

into 
fai, +--+, an, 0, -°+, O]eU ~ (4) 
in the frame of 19.381. Then D is real and D* in the chosen frames takes 
[br, +++, OsleJ (5) 

into 
[br +++, DnleKe. : (6) 


19.42 We interpret the n-tuples (3) and (6) as voltages and currents 
in the external meshes of N. Their relations to (4) and (5) are con- 
sistent with this interpretation. 

Let us restrict M by D, to get a correspondence M, between V2 and 
K, . In our chosen frame, the passage to [| corresponds, by (8) and 
(4) of 19.41, to considering mesh voltages in N which vanish for every 
internal mesh, and, correspondingly letting the mesh currents adjust 
themselves to this situation. We of course observe only the external 
mesh currents (6). 


19.43 M was PR. So, therefore is M, (18.3 dual). Since [14,] is the 
correspondence established by the physically realizable 2n-pole N, 
the necessity of Pl, --- , P7 for formal realizability is established. 


XX. APPENDIX TO PART I 


20.0 We must prove 7.22 and those assertions of 10.6 which are not 
covered in Halmos’. These concern reality. 


20.1 Let Vi be a real manifold and 
V=Vi@V:, K=K, © K, 


where K, = (V2), etc. The basis (14) of 10.6 exists by Halmos’, par. 19. 
We show that it can be chosen to be real. We have linearly independent 
vectors 


Vp, °°' 4 Ury Urg1, °° 5 Vay 
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where the first r span Vi, the last n — 1, Vo. Let 
Vs = Us + 1; , l<s<n, 
where u; , w; are real (10.42). Since Vj is real and a linear manifold, 
Us = 3% +0)eVi, Los<r, 
and, similarly, w.eVi, 1 < s < r. Among the 2n real vectors 


U1, Ua, °c? Ur, Wi, et? y Wry Urq, 6° * y Un, Wry, 2°? Wn, (1) 


the first 2r are in V; , and they span V; because the v,, 1 < s < r, can 
be constructed from them. The whole list (1) spans V, because from it 
all the »,, 1 < s < n, can be constructed. Since the v,eV2 do not use 
in their construction any of the first 2r vectors (1), it follows that the 
last 2(n — r) vectors in that list must contain a set spanning V2. The 
reality of the vectors (1) then establishes the existence of a real basis, 
say, 

Te ee cc (2) 
which provides a basis in V; and V2. 
20.11 We now have 7.22. The unique dual basis 

ky, °°, kn 

to (2) is real by 10.41. Hence all of Vi, V2, Ki , Ke are real. The proof 
of 10.6 is then complete. 


20.2 If in a real basis (2) (dropping primes) 


V = AW + Ave + +++ + Anda, 
that is, if 
[vo] = [a , +++, ar, 
then by (5) of 10.8 
d= G1 + +++ + Gna, 
hence 
[6] = [a1, +--+, Gal. 


The geometrical conjugation of 10.3 is therefore simply the concrete 
one of 7.2 in any real basis. This proves the remark of 10.35. 
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An Application of Boolean Algebra to 
Switching Circuit Design 


BY ROBERT E. STAEHLER 
(Manuscript received January 10, 1952) 


This paper discusses the application of switching (Boolean) algebra to the 
development of an all-relay dial pulse counting and translating circuit em- 
ploying the minimum number of relays. An attempt is made to outline what 
appears to be the most promising method of obtaining beneficial results from 
the use of the algebra in the design of practical switching circuits. 


INTRODUCTION 


The demands made upon telephone switching systems in regard to im- 
provements in handling capacity, speed, flexibility and economy are con- 
tinually increasing. In order to meet design objectives enabling the 
fulfillment of these demands, switching circuits have of necessity become 
more and more complex and intricate. As certain types of relay switch- 
ing circuits increase in complexity, the problem of control and output 
contact network design becomes more and more laborious and time con- 
suming. This is especially true in those circuits in which an attempt has 
been made to achieve the ultimate in efficiency and economy in that the 
number of relays used therein approaches the absolute minimum neces- 
sary to provide the required number of distinct output combinations. 
In this type of near-minimum combinational or sequential relay circuit 
there are numerous parallel control and output contact paths which 
thread through the same relays repeatedly, thereby causing the indi- 
vidual relay contact loads to become relatively large. Thus the designer’s 
problem becomes that of first developing a workable control and output 
contact network and then manipulating and minimizing contacts within 
that network so that the maximum number of contacts used on any one 
relay is within that permissible on any commercially available relay 
having the necessary speed characteristics. 

Even in those combinational and sequential relay circuits which are 
not near-minimum and therefore probably have fairly light individual 
relay contact loads, there are, of course, advantages to be gained by 
using the least number of contacts possible. Although the initial cost per 
additional contact (assuming that a few added contacts per relay will 
not impair the relay speed or space characteristics to an extent that the 
circuit requirements are not met) is almost negligible, there are other 
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economic savings possible. Since each contact must be connected to the 
remainder of the contact network, minimizing contacts and consequently 
soldered connections means a saving in wiring time and labor. Further- 
more, if the designer will manipulate the contacts so that the relays can 
be chosen from a comparatively few standardized codes, which are in 
large demand, it is possible to avoid the expensive stockpiling of numer- 
ous special designs having only a limited demand. In addition, using 
the least number of contacts minimizes the focal points of most relay cir- 
cuit failures which are the contacts themselves (i.e., dirty or worn con- 
tacts). 

It might also be noted at this point that electronic combinational or 
sequential circuits usually require electronic gating networks to perform 
functions which are completely analogous to those of relay contact net- 
works. Hence, the same problem of minimization exists. However, in 
electronic circuits, gate minimization is even more advantageous since 
the cost per additional electronic gate is much higher than the cost per 
additional] relay contact. 

It is rather obvious that the multiplicity of paths in most seibiie 
tional and sequential circuits can cause their design to become an ex- 
tremely difficult and time consuming problem if the contact paths are 
developed with the aforementioned considerations in mind. 

The circuit designer’s usual approach to the solution of such contact 
minimization and manipulation problems is that of inspection. The 
method of inspection presupposes a background of considerable experi- 
ence in that the designer must recognize certain contact network arrange- 
ments that may allow further rearrangements and thereby he must 
mentally develop his own rules. In order to check on any of his manipu- 
lations he must repeatedly redraw the network during this inspection 
design process. It is evident that this is often a long and tedious method 
and, depending on the skill of the designer, may or may not result in an 
optimum or even adequate solution. 

Suitable contact network arrangements often appear only after con- 
sideration of several alternative schemes and the rearrangements of the 
network interconnections of these schemes. Realization of this makes it 
quite evident that any means of obtaining and comparing these various 
schemes quickly and with a mathematical accuracy which does not 

“require continuous checking of network paths permits a more rapid and 

complete exploration of the particular problem. Switching algebra, first 
codified by C. E. Shannon!, is the systematic application of G. Boole’s? 

1C, KE. Shannon, A Symbolic Analysis of Relay and Switching Circuits, Trans. 
ATEE, 57, 1938. 


2G. Boole, The Mathematical Analysis of Logic (Cambridge 1847) and An In- 
vestigation of ‘the Laws of Thought (London 1854). 
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“Algebra of Logic” to switching circuits and is just such a means. It is 
a tool which can be used to investigate the complex combinational and 
sequential networks to determine satisfactory contact arrangements or 
reject unsatisfactory ones with a minimum of time and effort. It should 
be emphasized, however, that as with any tool, satisfactory results de- 
pend upon the judgment, ingenuity and logical reasoning of the user. 
Furthermore, as will be evident from the following development, switch- 
ing (Boolean) algebra in its present state is not to be considered entirely 
selfsufficient but, for the most beneficial results, should be applied, when 
warranted, in conjunction with inspection techniques so that the latter 
may fill in any limitations in the algebra techniques which have not been 
completely systematized as yet due to the newness of this field. 

The problem of solving the contact requirements of a minimum relay 
dial pulse counting and translating circuit recently developed as a com- 
ponent of the originating register of the No. 5 Crossbar System will be 
used as a means of illustrating the practical use now being made of 
switching algebra and of indicating exactly where the application of the 
algebra enters the design problem. 


BASIC DIAL PULSE COUNTER REQUIREMENTS 


The primary function of the originating register is to receive pulse 
signals representing digits from a telephone dial or similar calling device 
and to store a record of the digits in a form suitable for use by an external 
circuit. The dial pulse counting and translating circuit, an integral part 
of the originating register, 1s oriented with respect to other parts of the 
register by the block diagram of Fig. 1. The LZ relay is the pulse detecting 
relay. When the subscriber’s switchhook contact is closed due to the 
lifting of the phone, the originating register is connected to the line and 
the L relay is operated. Thereafter it follows the breaks and makes of 
the subscriber’s dial and feeds these repeated dial pulses into the counter. 
After the pulses are counted they are translated to a new code. In switch- 
ing systems it is advantageous to translate from the basic dial ten pulse 
decimal code to a ‘‘two out of five” self-checking code. In this latter code 
any single error within the circuit will result in either one or three relays 
operated in the associated storage circuit rather than two and thus an 
error can readily be detected. The output of the translator is fed via a 
steering circuit to the register or storage circuit. The slow release RA 
relay is the pulse train detecting relay which holds between the indi- 
vidual pulses of a digit and releases, only at the end of the pulse train. 
When it releases it activates the translating circuit and thereby transfers 
the translated code information to the storage circuit. The RAJ relay 
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in operating terminates the output from the translator and simultane- 
ously releases the relays in the counter to prepare it for the next digit. 

Specific requirements imposed by the originating register circuit neces- 
sitate the counting of one to eleven pulses; the use of a driving source 
consisting of a single break-make (or transfer) contact with ground on 
the armature spring; and outputs as follows: 

1. Count of 1 through 10: ground on two of the 0, 1, 2, 4, 7 puta 
leads in the combination corresponding to the count. 

2. Count of 10: ground on the ZO lead. 

3. Count of 11: ground on the 0 lead only (this is a trouble-detecting 
feature). 

In addition, the design of the steering and register-storage circuit 
requires that no output leads be connected together until the second 
pulse is received. Furthermore, each relay is limited to a combination of 
simple make and break contacts not exceeding a total of twelve. This uti- 
lizes the maximum number of springs obtainable on presently available 
relays and also avoids the larger armature gaps imposed by transfers 
which would result in a reduction in the relay speed of operation. Speed 
requirements also do not permit the use of shunt release in the circuit 
operation. 












COUNTER 


RELAY COILS AND 
ASSOCIATED CONTACTS 
FOR PRODUCING 
REQUIRED OPERATING 
SEQUENCE 






TRANSLATOR 2. VIA STEERING 


CONTACTS ON CIRCUIT TO 
COUNTER RELAYS A. 7 REGISTER OR 
FOR TRANSLATING none 


RELAY SETTINGS TO 
DESIRED OUTPUT 
CODE FORM 


Fig. 1—The schematic of a portion of a dial pulse register circuit for counting 
decimal code pulses and translating them to ‘two out of five’’ signals. (In the 
symbolism used in the illustrations a cross indicates a “make’’ contact and a 
vertical bar indicates a “break’’ contact.) 
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THEORY OF A MINIMUM RELAY COUNTER 


The counting circuit under consideration does not contemplate the 
use of any circuit elements other than relays that react to the beginning 
or end of a pulse. Therefore it must establish a distinct combination of 
relays operated or released during and between successive pulses. The 
minimum number of ordinary “two-position” relays, R, required to 
count P pulses can be obtained from the expressions (1) 2P < 2” if 
the counter is to lock up during, or recycle after, the last pulse or (2) 
2P < 2" — 1 if the counter is to lock up after the last pulse. 

The usual counting circuit used for determining the number of pulses 
in a dial train is required to count ten pulses, however there are certain 
advantages in regard to trouble indications if the counter counts eleven 
pulses. In either case the minimum number of relays necessary, accord- 
ing to the preceding formulae, is five. It should be noted that the ease 
with which this minimum number can be attained depends upon whether 
the input is derived from a single, double or transfer contact source. 


DETERMINATION OF OPERATING SEQUENCE 


Having determined that the minimum number of relays necessary is 
five, the first step in design is to develop an operating sequence pattern 
from the resulting 2° or 32 possible relay combinations. These combina- 
tions may be utilized in any order deemed desirable to obtain the 23 dis- 
tinct combinations needed to differentiate between eleven pulses (22 for 
the eleven makes and breaks plus an all-relays-normal combination). In 
this phase of the design switching algebra is not involved. The optimum 
sequence to meet a particular set of requirements can only be determined 
by repeated trials guided by an intimate knowledge of objectives. 

Initial studies, made by Joseph Michal, of various possible sequence 
patterns for a five relay circuit, including those having a three relay 
“ring” followed by two auxiliary relays and those having a two relay 
pulse divider followed by three auxiliary relays, resulted in the conclu- 
tion that the latter approach was the most fruitful. The sequence pat- 
tern adopted is shown in detail in Table I. The pattern is extended 
through 12 pulses, and it can be seen that the nature of the sequence is 
such that this employs all 32 combinations of the 5 relays. Several of 
these are transient and occur during part of a pulse or inter-pulse inter- 
val. Examination of the tail end of the sequence indicates that it will be 
simpler to design on the basis of a full 12 pulses than attempt to block 
at the end of the 11 pulses specified by the requirements. If trouble con- 
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TABLE I 
SEQUENCE OF OPERATION 





























Pulsing Counting Relays 
Relay Relay. Two out of 
Combination} Five Code 
L A B Cc D E 
Seizure 0 1 1 1 1 1 1 
Ist pulse 1 0 1 1 1 1 2 
0 0 0 1 1 1 3 0,1 
2nd pulse 1 1 0 1 1 1 4 
1 1 0 0 1 1 5 
0 1 1 0 1 1 6 0,2 
3rd pulse 1 0 1 0 1 1 7 
0 0 0 0 1 1 8 
0 0 0 0 0 1 9 152 
4th pulse 1 1 0 0 0 1 10 
0 1 1 0 0 1 11 0,4 
5th pulse 1 0 1 0 0 1 12 
1 0 1 1 0 1 13 
0 0 0 1 0 1 14 1,4 
6th pulse 1 1 0 1 0 1 15 
0 1 1 1 0 1 16 
0 1 1 1 0 0 17 2,4 
7th pulse 1 0 1 1 0 0 18 
0 0 0 1 0 0 19 0,7 
8th pulse 1 1 0 1 0 0 20 
1 1 0 0 0 0 21 
0 1 1 0 0 0 22 1,7 
9th pulse 1 0 1 | 0 0 0 23 
0 0 0 0 0 0 24 
0 0 0 0 1 0 25 2,7 
10th pulse 1 1 0 0 if 0 26 
0 1 1 0 1 0 27 4, 7-ZO 
11th pulse 1 0 1 0 1 0 28 
1 0 1 1 1 0 29 
0 0 0 1 1 0 30 0 
12th pulse 1 1 0 1 1 0 31 
0 1 1 1 1 0 32 0 











Total of 25 = 32 combinations used. 
Note: 0 is used to indicate that the relay listed at the head of the column is 
operated, and 1 is used to indicate that the relay is released. 


ditions introduce pulses beyond 12, the circuit will without difficulty 
recycle through combinations corresponding to pulses 11 and 12. 
Table I also indicates the leads which must be grounded in order to 
provide the translations to the “two out of five” and “single lead’’ codes. 
The characteristics of this circuit may be summarized as follows: It 
contains only five relays which is the absolute minimum necessary. It 
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uses all 32 of its available combinations. Its control and translating job 
is complex enough to indicate the need for a considerable number of con- 
tacts and hence the need for extensive contact manipulation to minimize 
and distribute these contacts. 

It is apparent that a great deal of time would be necessary to accom- 
plish this manipulation by inspection methods, thereby indicating the 
need for an additional tool such as switching algebra to assist the oe: 
signers in this task. 


ALGEBRAIC METHODS APPLIED TO CONTROL CIRCUIT 


The sequence of operations of Table I is used as the starting point 
in the application of the algebra. The exact calculations necessary to 
develop the control and translating circuit by this means are shown in 
detail later. However, the individual steps in the solution might well be 
outlined here. First, the design of the control and translating networks 
will be regarded as separate problems. In theory these can be integrated 
together, but the resultant network is likely to be so complex that under- 
standing and maintenance of the circuit would suffer. Each of the two 
networks can be individually considered as a multi-terminal network of 
the single input type. That is, the control network is an associated set of 
contacts which connects a single ground input to the windings of five 
relays, and the translating network is an associated set of contacts which 
connects a single ground input to the six output leads. Since switching 
algebra is directly applicable to two-terminal networks rather than multi- 
terminal networks, the approach to this particular problem is of neces- 
sity somewhat indirect. 

The most satisfactory method of attack is to develop first a two- 
terminal network for each of the output paths of the multi-terminal 
network under consideration. The two-terminal networks can be ex- 
pressed algebraically and manipulated into their simplest form by means 
of the switching algebra theorems to be given later. The individual net- 
works can then be inspected carefully, either in algebraic or circuit form, 
with the objective of combining them in the most advantageous fashion. 
It will be found, in general, that the simplest network configurations do 
not readily combine and that further manipulation is necessary to obtain 
an economical circuit. It is at this point that the algebra achieves its 
greatest utility, since its application permits the simple and rapid chang- 
ing of a given two-terminal network into a large variety of different 
forms with mathematical assurance that circuit equivalence is main- 
tained. Inspection of the networks in the several forms provides clues 
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as to the preferable combining forms and often indicates additional 
manipulations that might be desirable. 

This network development is a combination of mathematics and inte- 
gration by inspection. It is characterized by repeated trials of alternative 
forms and at no stage is there any definite assurance that the optimum 
circuit has been attained. However, the ease of manipulation provided 
by the algebra greatly enhances the probability of designing a better 
circuit than would be possible by inspection alone. In combining the 
two-terminal networks, care must be taken not to introduce “sneak”’ 
- paths which improperly connect outputs together. The algebra usually 
offers means of introducing one or two additional contacts which permit 
combining networks and yet eliminate the adverse effects of the sneak 
paths. 

The above procedure will now be carried out in detail with the switch- 
ing algebra theorems that are used in all the following algebraic manipu- 
lations noted at the margin by the number which corresponds to the 
number of the theorem in the complete listing in Table II. This table is 


TaBLe II 
SwitcHING (BOOLEAN) ALGEBRA* 














Definitions Postulates 
Addition (1) X = 0 or X = 1, where X is a con- 
(+) = AND = Series tact or a network. 
Multiplication (2a) 0-0 = 0 
(-) = OR = Parallel (2b)1+1=1 
(8a) 1-1 =1 
Circuit States (3b) 0 +0 = 0 
0 = Closed Circuit ; (4a) 1-0 = 0-1 = 0 
1 = Open Circuit (4b)0+1=1+0=1 
Theorems 
(la) X + Y=Y+X (8a) X’ +X =] 
(lb) XY = YX (8b) X’X = 0 
(Qa) X++Y¥Y+2Z2=(X + Y)+2Z (9a) O+ X = X 
=X + (Y + Z) (9b) 1-X = X 
(2b) XYZ = (XY)Z = X(YZ) (10a) 1 +X =1 
(8a) XY + XZ = X(¥ + (10b) 0-X = 0 
(8b) (X + Y)(X +2) =X+ YZ (ila) (X + Y’)¥ = XY 
(4a) X +X =X (lb) XY’+ Y=X4+Y 
(4b) XX = X (12a) (X + Y)(X’ + Z)(Y + Z) 
(5a) X + XY = X = (X + Y)(X’+ Z) 
(5b) X(X + Y) = X (12b) XZ + X’Y + YZ = XZ4 X'Y 
(6a) (X)’ = X’ (13) (XX + Y)(X’ + Z) = XZ + X'Y 
(6b) (X’)’ = X (14a) f(X) = A-f(X)aa1. a’=o 
(7a) (X + Y+ 24+ x + A’-f(X)ao, A’=1 
= X’-Y’-Z’. .«.-. (14b) f(X) = [A + f(X)aqo, a’=i) 
(7b) (X-¥+Z+ +++)! [AO ciate) 
SX YO Ae ee 








* Reprinted from The Design of Switching Circuits by Keister, Ritchie and 
Washburn with the permission of D. Van Nostrand Co., Inc. 
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taken from The Design of Switching Circuits by Keister, Ritchie and 
Washburn*. The development of the algebraic expressions from the 
sequence of operations table will be in exact parallel to the methods 
suggested in the aforementioned text. 

The symbolism adopted in the following development is basically that 
of using the notation A for all the make contacts on the A relay, and A’ 
for all the break contacts on the A relay. Contacts or groups of contacts 
in series are related by the symbol of addition (+) and contacts or groups 
of contacts in parallel are related by the symbol for multiplication (-) 
which may or may not be explicitly written, as in ordinary algebra. 
Therefore (A + B’) symbolizes a series contact path that is closed 
when the A relay is operated and the B relay is released, while (A B’) 
symbolizes the parallel contact path that is closed when either A is 
operated or B is released. Switching algebra includes only two numerical 
values, 0 and 1, with the quantity 0 assigned to represent a closed path 
and 1 to represent an open path. For the tabular notation of Table I, 0 
is used to indicate that the relay listed at the head of the column is oper- 
ated and 1 is used to indicate that the relay is released. 

As stated earlier, the present application of switching algebra utilizes 

‘the sequence of operation chart of Table I. The operate and release 
combinations for controlling the A, B, C, D and E relays can be selected 
from this table by observing where each relay to be controlled changes 
state. For example, the operate combination for relay D is relay combi- 
nation 8 and the release combination for relay D is relay combination 24. 
It is not necessary to include the contacts of a relay in its own operate 
and release combinations. Note that the A and B relays which serve as 
a pulse divider can be controlled solely by the Z relay and contacts on 
A and B without reference to C, D, Z. However the C, D and E relays 
are internally controlled by all five counting relays. The development 
of all these control paths uses the following abbreviations: 


g(X) = operating combinations for the X relay 
r(X) = releasing combinations for the X relay 
h(X) = holding combinations for the X relay 


X = make contact on the X relay 


Furthermore as expressed by theorem (6a and 6b) the negative of a 
contact network X is defined as a network which is a closed path under 
all conditions for which X is open, and is open under those conditions 


*D. Van Nostrand, 1951. The Bell Telephone Laboratories Series. 
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for which X is closed. Hence h(X) may be obtained from r(X) by noting 
that h(X) is the negative of r(X). Therefore the entire control path of 
any relay can be expressed generally as 


F(X) = G(X)[X + h(X)] = g(X)EX + (7 (X)) 
Thus for the A relay 


g(A) = L’ + B’ 
r(A) = L'+ B 
h(A) = [L' + BY = LB’ (7a) 


and 


F(A) = (A! + B)(A + LB’) 


= (i + BA + LA + B’) (3b) 

= (L+ A)’ + B’) (12a) 
Also for the B relay 

g(B) =L+A 

r(B) = L+ A’ 

h(B) = (2+ A’! = LA (7a) 


and 
f(B) = (LE + A)(B + L’A) 
Sir AB LB oe A) (3b, 
= (L + A)(L' + B) (12a) 


The schematic forms of the A and B control circuits as represented by 
the above algebraic expressions are shown in Fig. 2a and 2b. Since the 
general requirements of the basic problem specify only one transfer on 
the L relay, only simple makes and breaks on the A and B relays and no 
shunt release paths (to avoid reduction in speed of operation), the combi- 
nation of the above specific circuits is not possible without recourse to 
double windings. Another factor which affects the practical form of the 
circuit is the finite transit time of the L relay armature spring. Switching 
algebra presupposes instantaneous action of relay contacts and in certain 
cases, when the use of a break-make transfer is required, additional con- 
tacts are necessary to cover the open contact interval. The final circuit 
form, conceived by F. K. Low, is shown on Fig. 2c and uses an added A 
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contact. Algebraic equivalence of this circuit with the original is shown 


below. 
f(A) = (L'A + B')(L + A) (Fig. 2c) 


= (L' + B’)(A + BY)(L + A) (3b) 

= (i BU A) (12a) 
f(B) = (A + BYE + A) (Fig. 2c) 

= (L' + B)(A + B)(L + A) (3b) 

= (L' + B)(L + A) (12a) 


Yor the C relay which operates and releases twice in the entire 32 
combination cycle 


g(C) = (A'+ B+ D'+ BY\(A'+B4+D+E£) 
(C) =[(A+t+B4+D4+P)(A+ B'+D'4+ £) 
A(C) = [(A+ B'+ D+ E)(A + B+ D'+ By 
= (A’BD'E + A’BDE’) (7a, 7b) 
and 
f(C) = [(A' + B+ D'+ BY(A'+ B+ D+ B)I 
[C + A’BD'E + A'BDE’ 
= [A'+ B+ (D'+ E’\(D + EB) 
[C + A’B(D’ + E’)(D + EB)] (3a, 3b, 13) 
= [A’+ B+ (D+ B)\(D + ENC + A’'B) 
[C+ (D'+ ED + B) (3b) 
= [(A’ + B)C + (D’ + E’)\(D + BE\IIC + A’B] (3b) 


The schematic circuit which the above represents is shown in Fig. 3a. 
Circuits of this type which use certain contacts more than once can some- 






LB tA) L 3 (A) 
aaa Pai, ial 
2 (a) Lt 

i 

LB (B) - 
ct / 

—~ LA II_ x IL 
are ABER NEE, 


= (L) RELAY CONTACTS 


(b) (c) 


Fig. 2—Pulse divider of counting circuit. 
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times be drawn in bridge form with a gonsequent saving of contacts. 
One method is to manipulate the expression into a form which is known 
to be the series-parallel equivalent of a bridge. However, following 
usual algebraic procedures it is often difficult to recognize where this is 
possible. In the present case a method developed by G. R. Frost (not 
yet published) was used effectively. This resulted in the bridge circuit of 
Fig. 3b which has the series-parallel equivalent: 
fC) = (C+ D+ B)\(D + LA’ + B+ DW + HYD + £) 
[C + AIC + B] 
By use of theorem (8b) this is seen to be equivalent to the previous 
expression for f(C). 
For the D relay which only operates and releases once in the entire 
cycle 
gD) =(A+B4+C+ FE’). 
r(D) =(A+B4+C+ 2) 
h(D) = (A+ B4+C4 B)’ 
— A'B'C'E' 
and f(D) = (A+B+C+ H’)(D+ A'BC'F’). 


AB DE 
oe (c) po a Bc (E) 
x 
tL ae Iq a E Pay 
A = wiz = 
c : (d) . 





Fig. 3—Internal control of counting circuit. 


292 THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1952 


By noting that A + B + (Cis the negative of A’B’C’, this can be 
reduced to 
f(D) =(A+BH+CO4+ B)\(D + A’B'C) (3b, 12a) 


However, this transformation introduces a hazard caused by the 
transit time of A relay contacts in passing from relay combination 9 to 
10. Therefore the original expression will be used for relay D. The control 
path is shown in Fig. 3c. 

For the E relay which operates and locks only once in the cycle 


gf) = (A’ + B’ + C’ + D) 
h(L) = ££ 
and 
fH) = (A! + BY + C' + D) + B) 
(A4’+ B+ C’ + D)E (4a) 


This control path is shown in Fig. 3d. 

Apart from the problem of developing the required contact network, 
the practical problem of what operating power must be given to the 
relays in order to meet speed requirements must be dealt with. Since 
the use of low resistance windings in series with protective external 
resistors is called for to obtain the speed required, it appears that the 
use of two windings per relay might prove advantageous. By operating 
on the low resistance winding while locking on the high resistance wind- 
ing, the current drain may be reduced (thereby saving a fuse) and 
furthermore some code reduction may be made possible as shown later. 
If double windings are used, two of the external series resistors may be 
eliminated by combining the control network so as to make certain that 
only one of the low resistance windings on the C, D, or & relays is 
energized at any one time. This would permit the use of one common 
external resistance with the aforementioned relays instead of three. 

Keeping these practical considerations in mind, further savings may 
be made by combining the control circuits as shown in Fig. 3e. Although 
there is in this circuit a possibility of contact stagger on the A relay 
contacts causing the C and D low resistance windings to be energized at 
the same time, this will not be harmful since, when the stagger occurs, 
both relays are firmly locked operated by their high resistance holding 
windings. 


TRANSLATING CIRCUIT 


The translating circuit is particularly adaptable to switching algebra 
manipulation. Table I shows the combinations which prevail at the end 
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of each pulse and the necessary ‘‘code”’ leads that must be grounded at 
these times. Reference to the block diagram of Fig. 1 shows that the 
output of the translator is not activated until the slow release RA relay 
releases after the last pulse of a digit has been received. Therefore the A 
relay can be eliminated from these combinations since at the end of 
every pulse the A and B relays are either both operated or both released 
and hence only one is needed to indicate the condition of both. Table III 
lists the numerous combinations which must close a ground path through 
to each of the five code leads and the ZO lead. At the conclusion of the 
algebraic manipulation, the A and B contacts may be redistributed 
evenly since they perform interchangeable functions in translation. 
The objective in the design of the translating circuit is to obtain the 
most economical contact network subject to a spring distribution that 


TABLE III 
TRANSLATION 





Counting Relays 
Output Lead 




















Grounded Decimal Pulse 
B C D E 
0 0 1 1 1 1 
1 0 1 1 2 
1 0 0 1 4 
0 1 0 0 7 
0 1 1 0 11 
1 1 1 0 12 
1 0 ] 1 1 1 
0 0 0 1 3 
0 1 0 1 5 
I 0 0 0 8 
2 1 0 1 1 2 
0 0 0 1 3 
1 1 0 0 6 
0 0 1 0 9 
4 1 0 0 1 4 
0 1 0 1 5 
1 1 0 0 6. 
1 0 1 0 10 
7 0 1 0 0 7 
1 0 0 0 8 
0 0 1 0 9 
1 0 1 0 10 
ZO 1 0 1 0 10 





Inconsequential combinations that may be used for simplification 





a Oorcd 
eS Or© 
rocr 
rFPOoOrre 
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fits in with the control contact network. Again, this is a multi-terminal 
network problem and the procedure is to design two-terminal networks 
that combine most readily. Since it is impractical to illustrate all the 
repeated trials that led to the final design, each network will be designed 
separately with the understanding that some of the steps are imposed 
by the form of all networks viewed collectively. 

The procedure adopted for developing the “0” lead network is as 
follows. First set up the miniature table repeating the portion of Table 
III that corresponds to the “0” lead. These parallel combinations should 
then be manipulated algebraically to obtain the greatest simplification 
possible. It is rather easy to apply some of the algebraic rules by observ- 
ing the condition of the relay in the several combinations in the table. 
A simple “shorthand” rule to follow is: if in the table of combinations 
describing a particular two terminal network, all possible combinations 
of certain relays appear in conjunction with a single combination of other 
relays, the network contacts on the former relays may be neglected. In 
other words when 2” different combinations of any number of variables 
m, are identical in all but n columns, contacts on the corresponding n 
relays are not required. This procedure is carried out below. 


B.C DE 

0111 (B + C’ + D4 
ae B’ + C + EF 
100 1 

0 1 0 0 (B+ C’ + £) 
0110 , 
ee ed ee ee 


Thus we have the following algebraic expression for the ‘0’ lead, 
which can be simplified as shown. 


(B+ C4 D/P +C4+ L)\(B+C4 ENC + D’'+ EB) 
[C’ + (B+ D')\(B + E)(D' + £)\(B’ + C + EF’) (3b) 
IO PE BDY(B DOB G ee) (3b) 


This is shown on Fig. 4a. A somewhat different manipulation of the 
equation permits placing the network in the bridge form of Fig. 4b. 
The algebraic equation, given below, can easily be shown to be the 
equivalent of the original. 


[Bae O's BB. EDO CR +p Ca HUB ae Cred) 
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In certain cases the use of theorem (14b), normally employed to reduce 
the contacts of a particular relay to a single make and break, can pro- 
duce simplifications difficult to‘ accomplish otherwise. This is shown 
below, with the theorem applied with respect to relay E' since E’ tended 
otherwise to be heavily loaded. , 


(B+ C+ DB +C+ HYB+C' + EC' + D+ £) 


[E+ (B+ C’+ D’)\(B’ + C+ 1)(B+ C’ + 0)(C’ + D’ + O)I 
[Eo + (B+ C4 D)(B'+C4+0(8+C4 1000 + D’ 4+ 1)) 
(14b) 


(E+ C’ + BD')[E’ + (B+ C)(B+C’ + D*] (9a, 10a, 3b, 9b, 5b) 


By modifying the first factor of the final expression in accordance with 
theorem 11a, this equation can be put in bridge form as shown on Fig. 4c. 


[E+ C’ + BB’ + D')\E’ + (B+ CB + C+ DI 


The above equation uses the same contacts as the previous expression, 
and although the right hand member is in a slightly different form, the 
expression is equivalent to the one obtained earlier. 

When it is known. that output conditions are inconsequential for some 
relay combinations, these inconsequential relay combinations may be 
combined with valid combinations to eliminate contacts in the network. 
Inconsequential means that the output during these particular combina- 
tions does not affect the proper functioning of the circuit. Four such 
combinations are listed in Table III. Only those inconsequential combi- 
nations which will combine readily with the actual combinations, thereby 
resulting in a reduction in the number of contacts, are to be used. Al- 
though the use of all the all-relays-released condition may be helpful in 
certain cases, it will not be used in the circuit under consideration since 
its use makes the requirement that no tie shall exist between output leads 
until the second pulse is received hard to meet. 

With this in mind the “0” lead network is again examined. Note the 
use of another “shorthand” rule which states that if a part of the 2” 
possible combinations is used in closing a path, the negative of the unused 
part of the 2” possible combinations is equivalent to the original com- 
binations. Thus if in the case at hand three of the possible four combina- 
tions of the B and C relays occur in series with the same combination of 
the D and E relays, the expression used is that for the series path of the 
D and E relays plus the negative of the missing combination of the B 
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and C relays. In the following tabulation the combinations below the 
horizontal line are inconsequential. 


BBCODE 
Oa tA (C’ + D’ + BE) 
1 041 1 
1 0 0 1 
010 0 (C +:H’ + BD’ 
0 1 1 =O 
1 1 i 0 
et eee _—-® + D + ® 
0 0 0 0 


The expression becomes: 
(C+ BE’ + BD')(C’ + D’ + BEV B+ D+ E) 


[H+ (C+ 1+ B’D')(C’ + D’ + BO)(B + D+ 0)] 
[E’ + (C+ 0+ B'D’)(C’ + D' + B1)(B + D+ 1)) 
(14b on £) 


[E + (C’ + D')(B + D)IE’ + (C + BYD')(C’ + D' + B)] 
(9a, 9b, 10a, 10b) 


[EZ + (C’ + D’)(B + D)|LE’ + CD' + CB + B'D'| (8b, 9a, 4b, 5a) 
[E+ (C" + D’)(B + DIE" + BC + B'D' (12b) 
[E+ (C" + D')(B + DIE + (B+ D')(B + C)] (13) 


Fig. 4d shows the schematic of the above expression. It is possible to 
put this in a bridge form without other changes because of the manner 
in which the front and back contacts of D are related to the other con- 
tacts. Comparison of all the circuits of Fig. 4 indicates that they all use 
the same number of contacts although final decision should be post- 
poned until all the output circuits are obtained and the ease of combina- 
tion of the different circuits can be compared. 

The procedure for determining the remaining code leads is carried out 
on the following pages. 
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67? lead— 
BCODE 
6122 I~ 
OAR eee a 
010 1 
1 0 0 0 
001 1 
110 1 
000 0 (Cc + D + E&E) 


resulting in [ZH + C + D][E’ + B] which is shown on Fig. 5a. It will later 
be found advantageous, in combining, to include the B’ term in the first 
factor, giving the expression: 


(i+ B’+ C+ D)(E’ + B) shown on Fig. 5b 


“2” lead— 
B CODE 
102 I 
000 1 (Cc + D’ + E’ 
1 1 0 0 
0. 0 1 0 
0011 (Bo + C’ + D) 
110 1 
C 
0 0 0 9——_—— * rae 
hence 
(C + D’ + E")(B’ + C’ + D)(B+ C) 
or 
[C + BD’ + E)INIC’ + B’ + D (3b) 
For later ease in combining, this is changed to: 
[C + BOB’ + D' + EC’ + B’ + D (11a) 


shown on Fig. 5c. 
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2” lead— 
B CODE 
1 0 0 ee ee 
0 1 0 1—— SD + BH + BC 
1 1 0 0 
1 0 1 0O (Bi + C + D’ + B&B) 
001 41 
110 1 (BY + C’ + D) 
0 0 0 0 
hence 


(D+ HE’ + B’C’)(B’ + C’ + D)(B’4+C4+ D'4+ ##) 
[D + (B’ + C’)(B'C’ + E’)ID' + B’ + C + £] (8b) 
which is shown on Fig. 5d. 


““Y Tead-— 
BC 


co) 


a ee 


(C + E) 


oem CO [lH CO FH OC 
or Ooo oO CoO = 
CcOorR;FrFoOoC Yy 
or eS [CO CO CO CO 


hence . 
(B+ D+ E\C+ £) 
or 


E+ C(B + D) | (3b) 


which is shown on Fig. 5e. 


“ZO” lead— 
BBCODE 
1 0 1 O 


hence one has (B’ + C + D’ + £E) which is shown in Fig. 5f. 
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The final contact savings are achieved by combining the various out- 
put paths. The combined translation circuit that appears to be as re- 
duced as possible is shown in Fig. 6. Note that certain forms of the 
individual output paths combine more readily than others. For example 
Fig. 4c and 5b combine more readily with the remaining paths than 
Fig. 4b and 5a. Note also that sometimes it is not the most reduced form 
of the individual output paths that permits efficient combining. This is 











Non LEAD No” LEAD 


E 
x x (b) 
Wor" LEAD 


Nol LEAD 


Fig. 4—The “0” lead of the translating circuit. 
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x { LEAD x { LEAD 
= Te eC a =e sb c¢ 4 
Xone X x X———$X— 
a) (b) 
Wo" LEAD w; B CG Nye LEAD 
D 
Dee ome 
a (d) 
E B D No B or D E NS_~n 
xX x. x 7 LEAD x x ZO_LEAD 
= | Cc = 
* (f) 


(e) 
Fig. 5—The “1, 2, 4,7, and Z0” leads of the translating circuit. 
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exemplified by the use of Fig. 5b rather than Fig. 5a. Although various 
forms of all of the output leads were tested for efficient combination only 
the form used is shown for outputs other than the “0” and ‘‘1”’ leads. 
It is essential to scrutinize the final network for possible sneak paths. 
Sometimes to avoid these sneak paths it is necessary to add one contact 
on one relay to allow savings on others. Here again the inspection tech- 
niques go hand in hand with switching algebra and the need for both is 
obvious. The algebra obtains the various forms which are capable of 











“o”" LEAD 
RAI Non LEAD 
E "70" LEAD 


RAI 4" LEAD 


RAI "4" LEAD 


RAI Nate LEAD 


Fig. 6—Combined translating circuit. 


different degrees of combination very quickly and efficiently. The inspec- 
tion method is then necessary for the actual combination of these forms. 

The additional RA relay contact is necessary to assist in avoiding in- 
terconnections between the output leads until after the second pulse is 
received. The final assignments of either A or B relay contacts are 
chosen to equalize the load on these relays. 


THE COMPLETED CIRCUIT 


The final form of the counting and translating circuit 1s shown on 
Fig. 7. The relays are all double wound to gain the benefits of current 
drain reduction. One additional advantage of using double windings is 
the relay code reduction made possible since now only two codes are 
necessary. One code serves the A relay and one other code serves the 
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B, C, D and EF relays. In comparison to this total of five relays and two 
codes the circuit in present use in the latest crossbar system requires ten 
relays and seven codes. 


AN ALTERNATE DESIGN OF THE PULSE DIVIDER 


To illustrate the application of algebra where the apparatus contem- 
plated puts less premium on contact minimization but more on stand- 
ardization and winding minimization, certain modifications of the 
proposed circuit are considered. 

In the event that new apparatus developments make possible the 
construction of relays that meet the necessary speed requirements even 
though winding impedance is increased, it appears possible (if the pulse 
divider is redesigned) to use only one code having a single winding for 
all five relays. The use of added contacts might be allowable if the new 
type of relay carries more springs than the present relay. 

The redesign of the pulse divider to use single windings can be accom- 
plished by manipulation of the basic algebraic expressions derived earlier 
for the pulse divider. 

Thus for the A relay 


f(A) 


(Li + BY(A + LB’) 
[L’ + BA + BL + B)| (11a) 


By attempting to manipulate the B relay control circuit into the same 
form, one obtains 


(L + A)(B + L’A) 


f(B) = 
= (L + A)(L' + B)(A + B) (3b) 
= [L’ + BI[A + BL] (3b) 
= [E+ BA -- BG BD] (11a) 


The schematics represented by the above algebraic expressions are 
shown in Fig. 8a and 8b. The circuit of Fig. 8c is obtained by combining 
the first two circuits so that only a single transfer is needed on the L 
relay. Note however that it is necessary to make the lower two B trans- 
fers have continuity action to insure proper functioning. Fig. 8d shows the 
pulse divider drawn in conventional form. 


CONCLUSION 


As far as is known, the dial pulse counting and translating circuit 
described herein requires fewer relays than any other circuit with similar 
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ing the minimum number of relays. 


pulse counting and translating circuit employ 


Fig. 7—The complete all-relay dial 
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functions at present employed in Bell System standard switching equip- 
ment. The previous dial pulse counter used in the latest crossbar system 
required a total of ten relays. Thus the present design represents a 
considerable saving in cost and space. To a certain extent this result can 
be ascribed to the use of switching algebra during the circuit develop- 
ment. 

Relay circuits designed on the basis of utilizing a large proportion of 
the possible combinations permitted by the component relays usually 
require heavy spring pile-ups. Since general purpose relays are limited 
inthe number of springs which they can carry, this type of circuit usually 


L B (A) (B) 


bh Foe IK 


ik, A II 7 
oo x = So x tL 
I (a) aoe 





Fig. 8 


entails considerable design effort to make most effective use of the avail- 
able springs. Application of switching algebra to this aspect of the design 
problem can often provide crucial assistance. 

It is recognized that switching algebra, in its present state of develop- 
ment, does not permit complete mathematical statement and manipula- 
tion of multi-terminal networks as represented by the counting and_ 
translating circuits. It does provide, however, facilities in manipulating 
two-terminal networks into a variety of forms from which can be selected 
those that combine most readily. This can result not only in a saving of | 
time, but also in improved circuits which might not be realized by other 
design techniques. Unfortunately the algebra in its present state does 
not indicate when the optimum circuit has been attained. To some extent 
this is caused by apparatus or circuit considerations to which, since it is 
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concerned solely with contact networks, the algebra does not apply. 
Thus, there is still considerable room left for the ingenuity and judgment 
of the switching circuit designer. 

As a result of the experience in designing the dial pulse counter and 
translator, certain observations on the use of the algebra are believed to 
be valid. Although switching algebra may be used in the design of the 
simplest circuits, the most noticeable benefits are obtained by the appli- 
cation of the algebra to the design of those circuits in which the control 
and output paths are complex and interrelated. The particular minimum 
relay counting and translating circuit under discussion is an excellent 
example of this type of circuit. A secondary advantage of the algebra is 
its compact notation and its value as an efficient circuit “bookkeeping’”’ 
method. 
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Interaction of Polymers and Mechanical 
Waves 


BY W. O. BAKER AND J. H. HEISS 
(Manuscript received October 19, 1951) 


New techniques of Mason, McSkimin, Hopkins and co-workers for gener- 
ation of shear waves over the frequency range 2 X 10° to 2.4 X 10° eps 
have been used to study mechanical properties of chain polymers. Polymer 
solids, melts and dilute solutions, representing the main states in which 
plastics and rubbers are fabricated or used, were explored to find the char- 
acteristic relaxation times, rigidities and viscosities of various chemical 
structures. Polyisobutylene, hevea rubber, polydimethyl siloxane, vinyl chlo- 
ride-acetate copolymers and plasticized nitrocellulose were compared with 
polyethylene and polyamides as examples of the range of solid properties 
encountered. 

As melts, several polyisobutylenes, polybutadiene, polypropylene, poly- 
propylene sebacate and poly-a-methyl styrene were investigated as models 
for varying degrees of chain substitution. Chain rigidity in, for instance, 
polytsobutylene, seemed to reflect visco-elastic over-all configurational changes 
up through the kilocycle range, but nearest neighbor interactions took over in 
the megacycle region, leading to moduli of 10° dynes/cm? even for syrupy 
fluids. 

In dilute solution, polyisobutylene, polystyrene, natural rubber and poly- 
butadiene microgel exhibited characteristic dynamic viscosities and rigidt- 
ties depending linearly on concentration. Presumably, this reflects mechan- 
ical properties of tsolated chains. Some possible models were suggested for 
the frequency dependence of such properties. 


INTRODUCTION 


The “equilibrium” mechanics of polymers, the giant molecules of 
plastics and rubbers, have been quite elegantly developed in the range 
of high strains (“kinetic theory” of elasticity—Meyer,’ et al.). However, 
the molecular displacements as these strains, and, indeed, much smaller 
ones, occur, are little understood.” Nevertheless, it is essential to know 
about detailed motions in connecting chemical structure with physical 
properties. Only in this way can there be obtained from the chemical 
industry compositions which will serve properly in telephone apparatus. 
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Other studies have treated one way of getting at these mechanisms. by 
relating stress relaxation, creep, viscosity, etc. to a distribution of 
molecular relaxation times (and energy barriers), as originated by 
Kuhn.** * Another approach is to strain polymers with periodic waves 
over a very wide spectrum of wavelengths, eventually going to fre- 
quencies comparable with those of the thermal vibrations of significant 
groups or segments in the macromolecules. The resulting dispersion or 
resonance phenomena can then be examined. Hence a mechanical radi- 
ation field can interact with the masses of elementary structural units, 
as the usual electromagnetic field interacts with atomic and group 
charges. In general, direct interpretations of this kind must be done 
with shear waves, and, at least, not only with longitudinal or ultrasonic 
_ Waves. 

This kind of study is now proceeding using waves generated and fol- 
lowed by piezoelectric crystals connected in as actual electromechanical 
circuit elements (A. M. Nicolson, 1919). Recent schemes of Mason and 
co-workers cover the frequency range from 10 X 10° to 60 X 10° eps, 
as reported in the paper by Mason and McSkimin in the last issue, 
while a tuning fork method used by I. L. Hopkins has been applied to 
“soft”? polymers (rubbers) over the range 10° to 10° cps (the general 
range of J. D. Ferry’s work at Wisconsin on concentrated polymer 
solutions). 

The relation of these studies to the scientific and technical exploita- 
tion of plastics and rubbers is in knowing what a particular chemical 
composition does to strength, stiffness, ease of molding, impact tough- 
ness, etc. That is, are there qualities of the interaction of saturated 
aliphatic groups that make polyethylene or polyisobutylene have some 
glass-like as well as liquid-like, or rubbery, nature even at room tem- 
perature? If so, conditions causing brittle failures must be watched for. 
How is the storage of molecular strains in injection molded plastics re- 
duced by increasing molding temperature (when the kinetic theory stiff- 
ness per chain actually increases)? These and many similar problems 
may be generalized under the headings below; in each case the chemical 
structure of the macromolecule appears to be reflected in relaxation 
times which combine in different ways to give flow or rigidity, toughness 
or brittleness. 


Extrusion and Molding 


Non-Newtonian flow leading to “frozen-in” stresses, subsequent dis- 
tortion and irregular shapes of plastics’ and rubbers,” implies energy 
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storage in the sheared molecules. The dynamic shear studies will confirm 
this. Also dispersion of carbon black and other pigments is restrained 
by clastic qualities of “liquid” polymers (i.e., instead of ‘‘ mixing”, com- 
pounds just microscopically deform and later re-form.) Likewise, the 
efficiency of compounding’ and extrusion® depend on how quickly the 
molecules relax after straining. 


Impact Strength, Brittleness and Tenacity 


Toughness, mechanical shock resistance, ultimate elongation and 
strength reflect the facility with which the polymer molecules can be dis- 
placed without breaking the piece. Thus, they accommodate to the stress 
by motions presumably similar to those described above. (The situation 
is complicated when crystallites are also displaced.’) In any case, time 
sensitivity in the range 10° sec upward exists.”°'" The discussion by 
Morey” is a valuable survey of these ideas, and explicitly notes the 
significance of multiple relaxation processes on damping of shock waves. 
Evidence of the relation of simple changes in chemical structure to the 
principle relaxation times effective in these physical properties of plastics 
and rubbers is thus another part of the dynamics studies. The “brittle 
point”, or volume-temperature transition of amorphous polymers,” 
apparently reflects directly the correspondence of the time of experiment 
with dominant relaxation time of the polymer.’ ” A few measurements 
on plasticized polymethyl methacrylate (from which, however, no actual 
rigidities were calculated) indeed indicate abrupt stiffening as a function 
of frequency at a given temperature.’° However, the changes measured 
were too small] and indefinite to indicate any particular molecular relax- 
ation. Other work"’ with plasticized polymers is nevertheless concordant 
with the current findings that molecular relaxations and not long range 
order determine embrittlement. The converse of this is, of course, that 
as some ‘‘transition” is approached, hysteresis, heat build up, flex crack- 
ing and fatigue are greatest. 


Creep, Stress Relaxation and Recovery 


Even these “long time”’ qualities of plastics, such as found in cold 
flow, apparently result from integrated displacements of rapidly oscil- 
lating segments of the chains. A most interesting analysis of stress 
relaxation in rubbers employs Kuhn’s suggestion of a particular distribu- 
tion of relaxation times. The present point is that, again, these relaxa- 
tion times reflect processes which should appear directly in reaction of 
the polymer with high frequency shear waves. 
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From these aspects above the current results of dynamics studies will 
be reviewed. 


POLYMER SOLIDS: OVER-ALL MECHANICS 


Solid polymers will denote rubbers and plastics in the state in which 
they are technically used. This is usually their most complex form, with 
inter- and intra-molecular factors undistinguished. Thus, separation and 
identification of the main relaxation processes are difficult or impossible. 
However, it is interesting to consider typical values of modulus and 
viscosity as related to chemical structure, in the range of frequencies 
corresponding to extrusion rates, and stresses in actual use. 

These values of dynamic modulus and viscosity are distinct from the 
usual quantities in the literature: The usual expressions are for longi- 
tudinal (sound) waves, and give dynamic Young’s modulus” 


K* = By — th, 


FE, measures the out of phase part of the force-displacement relation, 
and LE, = w: (“effective viscosity coefficient’’). Now, the general elastic 
constants are \ + 2u, with \ = Lamé’s constant and » = shear modu- 
lus. Here, 


A+ Qu = K + 4p, 
with K = bulk modulus. Alternately, 
3K 3X + Qu 
~X+e" Ne 
However, in general the present results lead to the simpler shear modu- 


lus w. Further the energy losses studied are expressible directly as the 
usual shear viscosity 


iy 


w= 1. 
Previous comprehensive studies of the dynamics of rubbers over sig- 
nificant frequency ranges have yielded loss factors either written as 
E,/E, (see above),” or as a function of the shear viscosity based on 
Stoke’s assumption that the compressional (dilatational) viscosity is 
zero.” But as Nolle” and Ivey, Mrowca and Guth” clearly recognize, 
recent work has strongly manifested the presence of compressional vis- 
cosity in simple liquids” as well as polymeric ones.” * Hence, the pres- 
ent understanding relating molecular structure to viscosity, plasticity 
and visco-elasticity is unsuitable for interpreting mechanical wave mo- 
tion more complex than in shear, unless shear constants are also known. 
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This sums up to mean that the chemical interpretation of basic poly- 
mer mechanics requires shear wave measurements. Nevertheless, fas- 
cinating evidence of the existence of fine-structure relaxations in polymer 
solid has come from longitudinal wave investigations.?'”°*:” °° Also, 
the pioneering shear wave studies of Ferry and collaborators” * on con- 
centrated solutions of polymers have suggested intrinsic relaxations of 
the chain molecules in a highly plasticized ‘‘semi-solid”’ state. 

The more simplified findings cited below will be seen to unify ap- 
proaches in this field. Comment must first be made, however, on formu- 
lation of experimental results in dynamics of polymers. 


Expression of Dynamic Properties 


Alternate and equivalent expressions have been thoroughly surveyed ;” 
all represent combinations of either Maxwell (series) springs and pistons 
(elasticity and viscosity) or Voigt (parallel) springs and pistons. Obvi- 
ously, there is no physical separation of elastic and viscous elements in 
a polymer molecule, so the irrelevance of the detail of the model need 
not be emphasized. However, the models lead to convenient formulation 
of relaxation times which dielectric studies, in particular, have shown 
have clear connections with chemical structure. In this chapter, some- 
times one and sometimes the other model, or combination, will be used, 
with the symbols shown on the next page. 

Other symbols are sometimes used,” but should be easily identified 
in terms of the above. 


Rubbers and Soft Plastics 


In Table I, the shear moduli of rigidity, u, and of viscosity, yp’, are 
shown as calculated for the Kelvin-Voigt model, for polymers having 
the indicated units of structure. The frequencies are from a few hundred 
to a few thousand cycles, hence, in the range of much technical use, 
(flexing of tires ~300 cps) and rates of shear during processing.”’” Data 
are from a general study by I. L. Hopkins” of the Bell Laboratories, 
based on a tuning fork transducer introduced by Rorden and Grieco.™ 
The strains employed were always small, in the range 0.3 to 1.5 per 
cent; » and uy’ were essentially independent of strain, except for some 
loaded rubber stocks. The u values clearly trace the magnitudes to be 
expected in going from the most typical rubber (hevea) to the semi- 
rigid plastics (vinyl chloride-acetate copolymer and plasticized cellulose 
nitrate). As anticipated from steady-stress observations the ‘‘ plastics” 
have » > 10’ dynes/cm’. Increase of » with frequency is also greater as 
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the “plastics” range is approached; a relaxation region is implied. Figs. 
1 to 4 show the dispersion of rigidity with frequency in more detail. 
Especially striking in Figs. 1 and 2 is the small temperature dependence 
(at least between 27° and 66°C) of yu. Because of experimental uncer- 
tainty, » cannot be said to be actually higher at the higher temperatures 
in accord with straight kinetic theory, but at least it is strongly tending 
that way, as also noted for lower frequencies studies on natural rubber.” 
Nothing like this appears for the plastics; in plasticized nitrocellulose 
the 100-cycle rigidity decreases 10-fold from 27° to 66°C. This is, then, 
the second general dynamic quality which reflects the low van der Waals’ 
(dipole, dispersion and induction) forces in hevea rubber and polydi- 
methyl siloxane, as well as their intrachain flexibility. Interchain forces 
in polyisobutylene (Butyl rubber) are low too, but barriers to flexibility 
because of sterically hindered-CH; groups come in. Table I and Fig. 3 





KELVIN-VOIGT 


MAXWELL H 
{ 
o = strain 
S = stress 
t = time 
rt = relaxation time 
r’ = retardation time 
w= G = modulus 
! p’ = n = viscosity 
do 1a 8 do_S_ wo 
dt pdt 4 dt on 7 
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TaBLe I 
Shear Modulus, u, dynes/cm?2 Shear Viscosity 
27°C Poises, np’, 27°C 
Polymer Unit 
100 cycles 5000 cycles wae fore i 
Hevea rubber 3 X 10° | 5.5 X 108 350 40 
CH; 
—CH,—C=CH—CH.— 
Polydimethy] siloxane 0.7 X 105 1 X 108 300 30 
CH; 
eg 
CH; 
Polyisobutylene 5 X 10° | 30 X 10° | 8,000 | 1,500 
i 
—CH:—C— 
CH; 
Polyvinyl chloride (~92%)-acetate | 13 XK 106 | 80 X 10° | 25,000 | 2,000 
(~8%) plasticized by ~36% di- 
octylphthalate 
—CH,-—CH— and 
Cl 
—CH:—CH— 
O 
O—C—CH; 
Cellulose nitrate 60 X 10° | 250 X 108 | 80,000 | 4,500 
CH:ONO2 | 
CH——O 
==0 CH—O— 
CH-———CH 
ONOz ONOs 
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Fig. 1—Viscosity and shear modulus of hevea rubber (cross-linked). 
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Fig. 2—Viscosity and shear modulus of polydimethy] siloxane (cross-linked). 
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Fig. 83—Viscosity and shear modulus of polyisobutylene (cross-linked Butyl). 


emphasize that thinking about the mechanics of a particular chemical 
structure must include the spatial relationships of groups within the 
chains, as well as between them. 

The dynamic viscosities in Table I are also in accord with the se- 
quence of structures. Their frequency dispersion again connotes varying 
relaxation processes. Natural rubber’s low inner friction, for both com- 
pressional” and shear waves’° is famous in its low hysteresis heating. 
(This unique property is geopolitically crucial, because adequate truck 
and bus tires cannot yet be made of any other rubber.) Indeed, it is 
striking that at 100 cycles, a piece of gum rubber has a local viscosity 
of only 350 poises. The silicone rubber gum also has high elastic efficiency, 
and its temperature coefficient of viscosity is very low (see Fig. 2), like 
the thermal coefficients for familiar silicone liquids. It is exciting to 
speculate in Figs. 1 and 2, whether more precise measurements which 
Hopkins is now undertaking will confirm the apparently negative .tem- 
perature coefficients of viscosity at some frequencies. ‘‘ Kinetic theory 
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viscosity” arising from transfer of momentum among thermally agitated 
chain segments, does not seem to have been considered in the theory of 
perfect rubbers. As in gases, it would require an increase of viscosity 
with temperature. 

In polyisobutylene, however, the dynamic viscosity leaps upward in 
both magnitude and temperature dependence. It should be emphasized 
that this is, again, for a cross-linked (Butyl) gum—an infinite network 
like the hevea gum, with presumably infinite macroscopic viscosity. 
The striking thing is that this internal viscosity is not greatly dependent 
on the network, at the degrees of “‘cure” used in rubber technology. For 
instance, recent studies over the frequency range 20-600 cycles, on high 
molecular weight, M, = 1.2 X 10°, linear polyisobutylene,” give, at 
25°C and 100 cycles, u’ = 4800 poises, although the steady flow viscosity 
of this polymer at this temperature is greater than 3 X 10° poises.” 
Then, the infinite network (nsteaay flow — ©) Butyl polymer of Fig. 3 
has at 27°C and 100 cycles »’ = 8000 poises. At 1000 cycles agreement 
appears to be about the same, and is tolerable considering the several 
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Fig. 4—Viscosity and shear modulus of plasticized cellulose nitrate. 


316 THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1952 


per cent of compounding ingredients in the Butyl gum, and possibility 
of a small fraction of low molecular weight uncured polymer in it. 

Also, wide variations in the degree of cure of Butyl gums were studied 
without large changes in yp’. In this regard, the particular sample of 
Fig. 3 had an equilibrium swelling ratio (= volume swollen polymer in 
cyclohexane at 25°C/volume insoluble part of dry vulcanizate) of 4.84. 
This indicates an M,¢ value (average molecular weight between cross- 
links) of <20,000.” Actually many of the dynamic properties can prob- 
ably be found in individual chain units or segments even smaller than 
this. This is a significant point in engineering applications where plastics 
may be cured to reduce creep but where it is desired to retain typical 
‘““chain”’ properties to increase impact toughness. That is, usually some 
optimum condition for this compromise can be found. The later section 
on liquids will suggest that physical properites typically associated with 
chain polymers can indeed reside in even shorter chain sections than the 
M’s observed in usual gum rubbers. 


Filled Polymers 


Marked. effects of carbon black and other pigments are of course 
familiar in both steady and alternating mechanics of rubbers.’ ** °° 
Brief comment on their influence on dynamic shear properties and thus 
relaxation mechanisms involved may be directed toward plastics, also, 
however. Thus, technologically it would be desirable to load thermo- 
plastics with considerable volumes of “‘inert”’ fillers, Just as is done with 
rubber. But, almost invariably strength and toughness decline, instead 
of improving, as in the rubber case. A reason for this appears in inves- 
tigations by Hopkins when carbon black (a standard type of reinforcing 
black) was added to Butyl rubber of the sort described in Table I. It is 
that stiffness seems to rise more rapidly than internal viscosity—e., a 
given strain results in proportionately higher stress than the accompany- 
ing internal viscosity provides means for dissipating the stress (as on 
impact). Hence, the brittleness which fillers normally engender in ther- 
moplastics may represent this change in y vs. pw’ balance. Table II illus- 

















TaBLe II 
Wt. Per Cent Shear Modulus, », dynes/cm? at 27°C Shear Viscosity, u’, Poises, at 27°C 
of Carbon Black 
in Butyl Stock 100 cycles 5000 cycles 100 cycles 5000 cycles 
15.2 8 X 108 60 X 105 11,000 2,000 
28 .6 45 X 106 150 x 105 35 ,000 3,000 

















INTERACTION OF POLYMERS AND MECHANICAL WAVES 317 


trates some values for Butyl compounds. The swelling ratio (SR) for 
the compound containing 28.6 wt. per cent filler has dropped to 3.2, 
implying also considerable reductions in Mg (since theoretically 


a 
(1/SR)Y? = 7 2bM -”’). Thus, the apparent chain segment between 
Cc 


cross-links is shorter than in the unfilled stock (the two were cured to 
give closely similar degrees of primary valence cross-linkage) and corre- 
spondingly the steady-pull modulus is higher. Yet, the internal friction, 
while also higher, seems to reflect new relaxations from interaction with 
the filler, and total shock-absorbing power has declined. 


Microcrystalline Polymers 


The preceding studies at comparatively low frequencies indicated (1) 
magnitudes of shear rigidity and internal viscosity characterizing rub- 
bers and soft plastics. By familiar shifts of temperature or frequency, 
they would also apply to polymers known as hard, amorphous plastics 
at room temperature such as polystyrene and polymethyl methacrylate. 
(2) Dispersion of » and yp’ with frequency affirm that the intrinsic or 
fine structure relaxations have times <10™~ to 10° sec, and so refer to 
chemical units much smaller than the average molecules in the usual 
technical rubbers and plastics. A way to get at what sizes and habits 
these units might have will be by investigation of low molecular weight 
polymer liquids. But, while still in the section on solids, it is recalled 
that microcrystalline polymers such as polyethylene, polyesters (Teryl- 
ene), polyamides (nylons), cellulose esters, polyvinylidene chloride, poly- 
acrylonitrile etc., have mechanical properties dominated by their crystal- 
line-amorphous ratios.”’ ** “°’ ** The amorphous volumes are clearly those 
which donate the flexibility, toughness and shock-resistance of these 
plastics and textile fibers.” “° An interesting point is, how “viscous” are 
the disordered chain segments? In an over-all sense, all kinds of dissipa- 
tion including crystallite friction, analogous to grain friction in metals, 
scattering of longitudinal waves, and stiffening by low temperatures can 
occur in these polyphase systems. Thus, effects of chain orientation as 
well as Jateral order (crystallinity) have been detected in dynamics 
studies.” ** The intrinsically liquid-like or amorphous components of 
this behaviour—and the things which will correlate most simply with 
dipole concentration and other chemical features—are most accessible 
to study at very high frequencies. For, in these polymer solids, unlike 
the essentially continuous and homogeneous amorphous ones first dis- 
cussed, the mechanics reflect small regions having widely divergent 
properties. Thus, methods developed by H. J. McSkimin of Bell Tele- 
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phone Laboratories (described in the last issue), have been used to probe 
for elemental reactions at the upper end of the frequencies presently 
available. Both longitudinal and shear waves were used. In polyethylene, 
a wavelength for the shear waves was 0.0074 cm., at f = 8.55 X 10° 
cycles, and in polyhexamethylene adipamide (the usual 6-6 textile ny- 
lon), the shear wavelength was 0.0125 cm., for f = 8.67 X 10° cycles, 
all at 25°C. 

The important consequence of these experiments so far has been that, 
despite the small strains involved, the viscosity appears to be a “‘poly- 
mer”’ viscosity, rather than an inner friction involving just a few liquid- 
like atoms per unit. Thus, polyethylene of “‘equilibrium”’ crystallinity 
and average molecular weight corresponding to an intrinsic viscosity in 
xylene of [n] = 0.89 (at 85°C), was measured over the range from 0 to 
50°C. The results from both longitudinal and shear wave measurements 














TaBLe III 
Viscosity Poises 
Temp., °C, f, cycles/sec 
yp rv’ + 2p’ xv 
0 8 X 108 15 38 8 
0 25 X 105 5 14 4 
30 8 X 108 15 34 4 
30 25 X 108 5 13 3 








are given in Table III. These viscosities are expressed in this case for a 
Kelvin-Voigt model, of rigidity and viscosity in parallel. The rigidities 
associated with these viscosities are about 3 X 10° dynes/cm’, or not far 
from the value under steady pull of about 1 X 10”. 

- Now this suggests that the rigid plastic polyethylene retains, even 
under mechanical impulse of microsecond duration, a shock-absorbing 
capacity reflected in a shear viscosity of 5-15 poises, and a compres- 
sional viscosity of 3-8 poises. The former, yn’, may roughly correspond 
to the liquid viscosity of a paraffin-like chain of from 50 to 65 c-atoms 
in length. Thus, the dynamics measurements seem to relate to basic 
premises of polymer structure. These are that the amorphous regions 
(whose existence is shown quite independently by x-ray scattering, den- 
sity, heat-capacity, etc.) indeed take up and dissipate sudden stresses 
which the microcrystallites, despite their great strength, would be too 
brittle to sustain. 

These results give hope that further probing of the dynamics of liquid- 
like elements in rigid plastics will eventually lead to precise adjustment 


INTERACTION OF POLYMERS AND MECHANICAL WAVES 319 


of molecular weight, chemical structure (degree of branching in poly- 
ethylene), crystallinity, etc. These quantities, when fitted to a given 
pattern of uw, A, uw’ and ’ at proper frequencies would yield plastics of 
optimum serviceability under the multitude of stresses encountered in 
use. 

A similar liquid-like structure-even where the (crystalline) rigidity is 
much higher” and mobile chain segments smaller—apparently occurs in 
polyamides. Presumably the hydrogen bonding and dipole interactions 
are very imperfect in the disordered regions,” and there the chain inter- 
action is reminiscent of polyethylene. For instance, in polyhexamethyl- 
ene adipamide, measurements in the 8 to 30 megacycle range do indicate 
that the Lamé clastic constant \ is about 5.6 X 10" dynes/em’, but 
only about 3 X 10” for polyethylene. This reflects over-all stiffness 
dominated by crystallites. Nevertheless, the compressional viscosity, \’ 
is 17-6 poises (going from 8 to 30 mc) for the polyamide, but only 5-2 
poises for polyethylene. Of course, since there is dispersion in both cases, 
these relative magnitudes might be quite different at some other fre- 
quency or temperature (all above are at 25°C). Yet it remains that the 
nylon, despite its hardness, also has a liquid-like component more vis- 
cous than that of polyethylene. Similar relations appear in the shear 
viscosities, »’, also determined for these two systems. For the 6-6 poly- 
amide, u’ goes from 19 to 7 poises over the 8 to 30 me interval while 
polyethylene changes from 15 to 5. These quantities indicate again, as 
with the polyethylene, that ‘‘polymer liquids” rather than just a few 
small groups of atoms are the important mechanical elements even at 
frequencies of 10’. Now polystyrene, an amorphous polymer, also has 
rigidities of about 10° dynes/em’ but the »’ and \’ values at room tem- 
perature are far below 5 to 20 poises, and glass-like brittleness (although 
not so bad as silica glass) results. 

So far, then, the two characteristic extremes of polymer mechanics 
have been discussed: (1) the soft rubbers, whose dynamics at low kilo- 
cycle frequencies imply, at ordinary temperatures, predominantly over- 
lapping combinations of relaxation processes whose relaxation elements 
involve many segments per molecular chain; and (2) the hard, micro- 
crystalline plastics, whose behaviour is predominated by relaxation proc- 
esses having times of 10° to 10’ sec because the longer period (slower) 
displacements have been relaxed out at the temperatures of normal use. 
(Likewise, interconvertability by temperature” between these two ex- 
tremes is presumed. Also, a certain correspondence between dielectric 
and dynamic relaxations in these classes is indicated.”*) Next, it is in- 
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teresting to see what are the simplest structures (particularly in terms 
of molecular weight) yielding these effects. In other words, what kind of 
liquid really exhibits “polymer mechanics?” No detailed answer to this 
can be given below, but results on some polymer liquids of low average 
molecular weight will indicate that the mechanisms in rubbers and plas- 
tics are probably more general than previously supposed. 


POLYMER LIQUID MECHANICS 


By techniques described in detail elsewhere,””’ a series of polyisobutyl- 
ene liquids have been investigated. These polymers were made by ionic 
catalyzed mass polymerization at reduced, but not very low, tempera- 
tures. While no great care to purify the monomer was used, such poly- 
merizations require fair purity to go at all. Seemingly, the resulting liq- 
uids do represent a polymer homologous series, although head-to-tail 
sequence of the monomer units, some single ethyl rather than paired 
methyl! side groups, etc., may differ slightly from the higher molecular 
weight forms in Butyl rubber and polyisobutylene gum. Whatever are 
these details, it appears that the polymers represent a linear hydro- 
carbon chain, with essentially two methyl groups on every other chain 
atom: 


a eee 
—CH:2-C -; CH,-C—,— 
CH; | ele 


By contrast, polyethylene, with the nominal chain —CH 2:CH:2:-CHe:- 
CH;— and to a lesser extent polystyrene, . 


Cy: = ; 


have chains in which rotation about the bonds is less sterically hindered. 
The final section, on isolated polymer chains (in dilute solution), will 
consider this aspect further. However, some results will be reported be- 
low on a low molecular weight poly-a-methyl styrene, which may be 
considered structurally a cross between the rubber, polyisobutylene, and 
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the plastic, polystyrene, 
CH; 


tae 22 


Other studies in progress on liquid polybutadiene, polyisoprene, poly- 
propylene, and polypropylene sebacate from which further information 
about intra-chain stiffness may be derived, will also be noted. 
Properties of the polyisobutylenes studied are summarized in Table 
IV, some additional molecular weights in this range appear as extra 
points in some of the high-frequency graphs. The molecular weights M, 


, Sata : : 42, 42 , ‘ 2 : 
are “intrinsic viscosity” averages’ “” and, with reasonable estimations 


1 ee : : = 
of the il, ratio, check with cryoscopic number average, M, , values on 
n 
such materials, which are in turn listed in the table as expressed by melt 
viscosity relations of Fox and Flory.” These molecular weights repre- 

















TaBLe IV 

iso- 25°C Voigt 'Maxwell 
mie DP, My Mn viscosity Maxwell Bp B! Freq. Cycles 

Polymer Poises mn Poises | Poises 
A 10 565 318 0.37 | 3 x 108 0.6 0.6 | 14 & 108 
A” 30 1660 697 39.6 6.2 X 10° | 16.5 | 18.8 2X 104 
1.7 X 10° 7.9 | 10.0 | 14 & 108 
B 45 2520 | 1070 216 3 x 10° | 15.2 | 24.2 | 14 & 108 
C 56 3300 | 1720 737 3.6 X 109 | 20.2 | 47.9 ; 14 & 108 
D 74 4170 | 25380 | 1840 4.5 X 10° | 23.4 | 78.9 | 14 & 108 
E 147 8240 | 4850 | 4600 5.3 X 109% | 27.2 | 92.38 | 14 * 108 


























sent reasonable averages rather than absolute values for these hetero- 
geneous polymers. The DP values are just the number of isobutylene 
units per average chain. The y values are the steady flow viscosities at 
low rates of shear—usually determined by a falling ball. 


Rigidity and Viscosity Magnitudes 


The properties of these liquids ranging from polymer A having only 
forty times the viscosity of water to E, which begins to approach fluid- 
ities of technical polymer melts (polyamides, for instance), were explored 
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in the kilocycle range with shear waves generated by torsional crystals 
and in the megacycle region by shear waves with the reflectance method 
and by longitudinal (ultrasonic) waves from a pulse propagation tech- 
nique.” The results have been expressed in two ways. First, in earlier 
reports,’ ” a trend corresponding with experiment was given by two 
Maxwell elements arranged in parallel. This result is too simple com- 
pared to the distributions of relaxation times previously proposed for 
high molecular weight polymers to reproduce detailed observation. 
Nevertheless, perhaps because of the smaller molecules involved, there 
seems to be clear indication that two principal relaxations predominate 
the mechanical reactions of these liquids over the range of frequencies 
of present interest, 10° to 10’ eps. For example, for polymer D, these are: 








First Relaxation Second Relaxation 
fe ~ 4 X 103 cycles fe ~~ 5 X 108 cycles 
p~ 4 X 107 dynes/cm? uw ~ 6 X 10° dynes/cm? 





(In accounting for the second main relaxation, a hysteresis component 
had to be introduced whose significance has been suggested.”) 

Second, specific values of shear rigidity » (Maxwell) and yp (Voigt), 
shear viscosity up’ (Maxwell) and y’ (Voigt) as well as the constants for 
related compressional wave systems, \ + 2y (elastic) and \’ + 2y’ 
(viscous) have been calculated for particular frequencies. Unlike in the 
first way of expression, these latter quantities are all highly frequency 
dependent. However, they describe conditions at various frequencies of 
interest, and are thus often worthwhile. 

Both ways of looking at the data lead, as implied by the figures above, 
to the proposal that typical polymer stiffness (shear rigidity of ~10' 
dynes/cm’) is present at M, ~ 1600, with DP, ~ 30, or an average 
chain length of about 60 carbon atoms. This appears when the straining 
is done in 10 ° to 10“ sec. In the 10 ° to 10° sec range, rigidity occurs 
for even an average chain length of 20 atoms as shown in Table IV. 


STRUCTURAL FACTOR IN LIQUID MECHANICS 


The main relaxations in the kilocycle range in polyisobutylene liquids 
seem to lead to quasi-configurational elasticity. This is where the kinetic 
theory tendency for a most probable separation of chain ends is retarded 
by viscous interaction of segments between and within the chains. Hence, 
the middle dashed curves of Fig. 5, showing shear elasticity for some of 
the polymers of Table IV, decrease exponentially with increasing tem- 
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perature. While pure kinetic theory elasticity would give a modulus 
increasing linearly with increasing temperature, these systems, like all 
practical rubbers and plastics, actually grow softer with rising tempera- 
ture when deformed dynamically. It is striking, nevertheless, that a 
modulus of ~10' dynes/em’ seems characteristic of the visco-elastic 
energy storing of these simple polymer structures. As noted below this 
is 10° less than the crystal-like, close-packed, stiffness found for these . 
same molecular frequencies above their second principal relaxation time. 
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Fig. 5—Shear elasticities and viscosities of polyisobutylene liquids assuming a 
Maxwell model with relaxation frequencies 103-104. (Lower dashed curves for fre- 
quency-dependent models.) 
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Thus, it seems that the former, 10’, modulus is typical of the struc- 
tural arrangements in polymers said to be above their second order 
transition temperatures’ '* while the second, 10°, modulus reflects in- 
teractions below the freezing-in. 

These conclusions obtain regardless of the particular expression of the 
data. But, for comparison, curves are shown on Fig. 5 for a polyiso- 
butylene A” in which the dynamic modulus u at 20 ke is computed for 
both Maxwell and Voigt elements. The two points denoting the steady 
flow viscosity of polymer A” rank it with respect to the others in the 
Series. 

Apparently even very fluid polymer melts, chain molecule plasticizers, 
and small segments of long molecules must be expected to show appre- 
ciable rigidity when stressed rapidly. Referring to the introduction it is 
reasonable that rough extrusions, frozen-in molding stresses and the 
like are so easily produced. The lines of Fig. 5 are not, of course, implied 
to be linear over any considerable temperature range. In the region rep- 
resented, the temperature coefficient for viscous flow is about 16 kcal 
for the B, C and D liquids (about 12 for A). This agrees roughly with 
the steady flow values found for very high molecular weight polyiso- 
butylene.’ °* The temperature coefficient for the rigidity is less, as 
would be expected, since the whole center of gravity of the chain need 
not be displaced, but only local segments. 

This quasi-configurational elasticity is increased by molecular weight 
(although kinetic theory elasticity of chain segments in a network is de- 
creased by increasing segment length). The log u vs density at 25°C 
plotted in Fig. 6 indicates that the chief influence is the number of 
chains per cc, since the points for all the molecular weights now lie on 
a single line. It should be repeated that the elasticity modulus plotted, 
u, is again for a roughly frequency-independent or “absolute” model.” ” 
The same is true for the three solid lines on Fig. 6, showing y» in the 
second, or 10’ cycle, relaxation range. Here effects of detailed liquid 
structure come out; the three average molecular weights no longer lie 
so nearly on a single line. This elasticity is presumably from the crystal- 
like interaction of nearest-neighbor segments. If temperature is adjusted 
so that densities are the same, it is seen that the lower average molecular 
weight liquid has the higher elasticity modulus. This difference is not 
large, and should not be interpreted as showing an equal segment inter- 
action, for a polymer of lower specific volume (B compared to D). 
Rather, it emphasizes in this relaxation range, approaching the “glass” 
behaviour, that the relaxation rate is vastly more temperature dependent 
than the specific volume change alone, and structural variations in the 
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Fig. 6—Principal shear elasticities of polyisobutylenes as related to their densi- 
ties. 


Ob 
@ 
i) 
2 
@ 
bh 
° 
o 
o 


packing of segments in the liquid (coordination number, etc.) become 
important. 


“ Soft” and “Hard” Inquid States; Second Order Transitions 


A recent noteworthy study of volume-temperature and viscosity- 
temperature changes in polystyrene (with a note on polyisobutylene) 
brings out many points in common with ideas of polymer liquid struc- 
ture indicated by the dynamics work.”’” Particularly, the fact that 
according to steady state measurements, the ‘‘local configurational ar- 
rangement of the polymer segments’ below 7’, remains fixed accords 
with the postulations from dynamics work. That is, above the second 
main relaxation, it seems to be just the interactions in these fixed ar- 
rangements which cause the glassy (or ‘“‘crystalline”) dynamic modulus 
of 10° to 10° dynes/cm’. Further, the point that 7’, is not an isoviscous 
state for polymers” agrees with the dynamics result that macroscopic 
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viscosity of the polymer has relatively little to do with the actual values 
of dynamic viscosities. These would be at frequencies where the response 
of the polymer liquid to the mechanical field is determined only by 
motions within the local fixed arrangements mentioned above. 

Fig. 7 illustrates this, where on one scale the macroscopic viscosity is 
plotted according to the familiar log-log relation with molecular weight. 
Two extremes of average molecular weight, M,, and M, are used for the 
liquids, to show that the molecular weight distribution does not alter 
the general conclusions. (M, is an upper limit weight average figure.) 
On the other scale, the dynamic viscosity yu’, in this case for a single 
element frequency-dependent Voigt model, shows low values and marked 
curvature. These betoken the relaxation in which molecular weight, 
through its effect on free volume and other structural factors, is signif- 
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Fig. 7—Comparison of steady flow and dynamic viscosities (at 25°C and 8 mc) 
of polyisobutylene liquids of different molecular weights. 
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icant for displacements superficially quite different from those in macro- 
scopic viscosity. 

The compressional viscosity, \’ is also plotted, for the same model, in 
Fig. 7. It is, within experimental error, zero for polymer A’, as deter- 
mined by shear and compressional wave studies at 8 mc frequency.” 
This is a rare case, then, where the attenuation of sound waves through 
a liquid has been quantitatively accounted for by the shear viscosity. 
But, as soon as the average molecular weight rises to 1000 or so, 0’ 
comes in clearly, and the new mechanism for dissipating compressional 
or dilatational stresses is developed. As this presumably represents di- 
rectly free volume or coordination number changes in liquid struc- 
ture,” its detailed study near 7’, ,” and in connection with brittle 
points of rubbers, may eventually be especially fruitful. 

Another depiction of influence of average molecular weight in these 
liquids on dynamic viscosities occurs in Fig. 8. Here, the \, curve is 
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Fig. 8—Dynamic viscosities of polyisobutylene liquids as function of average 
molecular weight (25°C). 
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for, again, a crude model attempting to show compressional viscosity 
over the whole frequency range, while the other viscosities are Voigt 
expressions at 8 (or 14) mc. Extremes of molecular weight averages are 
shown. 

Comparison of the ‘‘soft”’ or quasi-configurational rigidities, expressed, 
like the » of Fig. 5 as relatively frequency independent yp. , with the 
“hard” or glassy rigidities is given in Fig. 9. The \ and yp values are for 
the Voigt model at 8 me. The graph does not show the bend-over of the 
“soft”, ue, curve with molecular weight, but that happens more grad- 
ually. The “hard” rigidities \ and » quite readily show this inflection. 
As before, the relaxing segments must be <100 chain atoms, according 
to the behaviour of the molecules at room temperature. 

Concerning influence of molecular weight on engineering ‘‘brittle 
points” of such importance in rubber technology, the present studies 
agree with earlier proposals. Thus, although the 7’, or v-7' second order 
transition point always decreases with decreasing molecular weight,” 
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Fig. 9—Dynamic rigidities of polyisobutylene liquids as function of average 
molecular weight (25°C). 
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the brittle point tends to increase as the molecules get smaller. This was 
supposed to be because the ultimate elongation (doubtless due to visco- 
elastic or quasi-configurational elasticity and not kinetic theory elas- 
ticity as sometimes said) declined with chain length, so that the speci- 
men broke at lower and lower strains,” even though it was not really in 
the glassy state. Now, the results above” ” demonstrate that the shear 
modulus at a given temperature does fall off for average molecular 
weights of polyisobutylene below ~5000. Hence, the mechanical em- 
brittlement of the low molecular weight samples is not because they are 
stiffer, but because they are weaker. 


POLY-a-METHYL STYRENE: A ‘“‘PLASTIC”? LIQUID 


Most of the liquid studies have been on polyisobutylene polymers, 
made “hard” or “soft”? by temperature or frequency, but under use 





TABLE V 
zp’, 14 me, Poises u, 14 me dynes/cm? 
Bee n Poises 
- Maxwell - Voigt Maxwell Voigt 
I a2 =| Md 23.6 |5.1xX 10°14 x 10 
II 4340 502.7 14.3 7.6 X 109 | 7.4 * 109 





conditions, considered rubbery. If one of the methyls in polyisobutylene 
is replaced by a phenyl, poly-a-methyl styrene, a hard plastic is pro- 
duced. Low molecular weight polymers of this composition are, however, 
liquids at room temperature. Hence, it is interesting to compare their 
reaction to mechanical waves with that of polyisobutylene liquids of 
similar macroscopic viscosity. Table V lists a few properties at 25°C. 
Polymer I has roughly the steady flow 7 of polyisobutylene B. Also, 
the Voigt yu’ at 14 mc is similar: 15.2 poises compared to the 23.6 poises 
of the poly-a-methyl styrene. However, the Voigt u is already 50 times 
higher for the phenyl substituted chain. (This shows the shift of the 
second relaxation range, of course, where nearest neighbor interactions 
rule.) Even more striking, the temperature coefficient of a ue , calculated 
as the second principal shear viscosity in a frequency independent 
model,” as before, is about 24 keal, compared to about 12 for yu, for 
polyisobutylene in its similar model. Thus, although there was no ap- 
parent difference between the mechanical properties of these two poly- 
mer liquids in their room temperature state, their dynamics diverged 
remarkably. This was when they were studied with shear waves whose 
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frequencies approached the glass-into-rubber relaxation times. Clearly, 
again, individual interaction of chain-like chemical units and not any 
micellar or other special aggregation of them, predominates polymer 
mechanics. 

It still remains, however, to separate interactions of the basic units 
within and between chains. Most likely, the model plastic vs rubber 
liquids just discussed differ in the high frequency region substantially 
only because of znter-chain forces between phenyl vs methyl groups. 
However, especially in the high-frequency region, questions of intra- 
chain structure, such as the steric hindrance of adjacent pairs of methyl 
groups in polyisobutylene, restricted rotations about bonds, etc., come 
in. Obviously where configurational or quasi-configurational displace- 
ments are important, as in all cases of elongation >20 per cent (this is 
certainly an upper limit), flexibility of single chains needs to be under- 
stood. This is built deeply into chemical structure; plasticizers pre- 
sumably may change over-all configuration as well as modify interaction, 
but they are impotent to vary flexibility. Accordingly, problems of rub- 
ber, usable in the Arctic, and of wire and cable insulation bendable at 
low temperatures always come back to whether the polymer chain bonds 
have free rotation. Some examples of the combinations of effects within 
and between chains can indeed be shown in several other polymer liquids 
which are rubber models. 

This influence of small changes in chemical structure is compactly 
illustrated by comparing a few other hydrocarbon polymer liquids with 
polyisobutylene. Also, rather dilute dipolar groups have been introduced 
in the linear polyester liquid polypropylene sebacate, whose structure is 
otherwise like that of hydrocarbons.” In Table VI, liquids of the given 
structure with some (unknown) distribution of molecular weights, were 
studied with shear waves at 77 and 142 ke at a temperature where each 
had the same steady flow viscosity. The figure chosen was 700 poises, 
and the temperature range required to adjust to it in the series was 
10.9° to 85°C, meaning that the liquids had comparable consistencies at 
ordinary temperatures. 

Despite these similarities under steady stress, the retardation times, 
r’, vary three-fold, with the highly substituted hydrocarbon chains, 
polyisobutylene and polypropylene, the highest. Despite the intermo- 
lecular action of the dipoles in polypropylene sebacate, the low polymer 
has a short retardation time, although its “brittle point”? with decreas- 
ing temperature is far above that of polybutadiene or even polyiso- 
butylene. Presumably the flexibility around C—O—C bonds rather 
compensates for increased dipole interaction. Where both low polarity 


TaBLeE VI—Dynamic Properties of Polymer Liquids of Varying Structure. 





Den- 


Polymer °C ay Pics 
Polyisobuty]- 
ene C 25.410.8857| 700 


Polypropylene |85.0| .8248) 700 
Polybutadiene |10.9| .8767| 700 


Polypropylene 
Sebacate 41.7/1.0421| 700 









































Voigt Maxwell 
77 ke 142 ke 4 77 ke 142 ke n 
r= G sec T= G sec 
G, G, G, 

Poines dynes/cm? Poises | G» dynes/em? 142 ke Poines dynes/cm? dynes/cm? 142 ke 
165.1/8.5 X 107)184.2| 6.4 X 107 | 21 x 10-7/197 |2.2 « 108 3K 108 | 6x 10-7 
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and chain flexibility obtain, as in polybutadiene and the silicones, dy- 
namic properties apparently accord with brittle points in implying small 
temperature coefficients of relaxation times. In fact, the temperature 
coefficient for dynamic viscosity of polybutadiene is only about 1.5 keal, 
whereas a comparable figure for polyisobutylene and polypropylene is 
12 keal. 

The frequency range in which the structural comparisons above were 
made, resides, as discussed earlier, in the zone of configurational visco- 
elasticity. That is, over-all shape changes, rather than just nearest 
neighbor interactions, are predominant even at these comparatively 
short average chain lengths. Now, other recent studies of polyisobutyl- 
ene liquids, at 5 to 100 eps frequency, exhibit no rigidity at 25°C and 
above, although they become non-Newtonian rapidly as temperature is 


TaBLEeE VII. Shear Dynamics of Polyisobutylene A’’. 
(1, = 1660; 742° = 39.6 Poises) 




















Voigt Maxwel 
T, °C Freq., eps 

uw dynes/cm? n Poises pw dynes/cm? n poises 
25 266 3.8 X 103 39.2 1.2 X 108 39.3 
25 1601 4.8 X 104 38.4 3.1 X 108 39.0 
27 25300 5.4 X 105 19.9 1.9 xX 10° 20.5 
27 41390 1.5 X 108 19.9 1.9 X 107 21.5 
27 53060 1.5 X 108 18.1 2.6 X 107 19.3 

reduced.” The questions are, where does the configurational elasticity 


drop out, as frequency is reduced at 25°C; and does it seem reasonable 
that this dispersion correlates with a shift in frequencies at lower tem- 
peratures. Partial answers are given by very recent studies of I. L. 
Hopkins of Bell Telephone Laboratories. He has equipped the tuning 
fork vibrator described earlier® with two parallel vanes filled in between 
with a film of polymer liquid. Pure shear properties can be derived from 
the response of this system. Table VII lists a few typical figures obtained 
on polyisobutylene polymer A”. These indeed show that the kilocycle 
relaxation zone (some new data by McSkimin’s torsional pulse method 
are given for it) extends smoothly down to where dynamic and steady 
stress viscosities are equal. Seemingly there are no new “extra long 
time” relaxation mechanisms; probably the slow relaxation times some- 
times indicated for high molecular weight rubbers are just displacements 
of this configurational relaxation to long times because of high molecu- 
lar weight and internal viscosity. 

By contrast to the conclusions associated with the data of Table VII, 
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some observations at low frequencies on isoviscous properties of polyiso- 
butylenes A” and C indicate nearly identical retardation times. Thus A” 
at 25°C and C at 61°C have yn, = 39.6 poises. The 7 values at 266 and 
1600 cps are also about 39 poises, » at 266 cps is 3800 dynes per cm” 
for both liquids, and at 1600 eps is from 3.5 X 10° to 4.8 X 10° dynes 
per cm’, 

In the final section, mechanical waves have been used to explore 
dilute polymer solutions, to see how isolated molecules behave, free of 
interaction with each other. 


DILUTE POLYMER SOLUTIONS 
Physical Principles in Measurements 


Precise information on dynamics of solutions approaching infinite 
dilution (and thus complete separation of the polymer chains) is desired 
here. Again these must be shear dynamics; bulk rigidity of ordinary 
liquids is so high that a few polymer molecules added cause little effect. 
Dilution is emphasized because even at 1 per cent by volume, high 
polymer molecule coils frequently interact, especially in “‘ good” solvents. 
Thus, several workers have detected shear rigidity in polymer solutions, 
in one case for polymethyl methacrylate of average molecular weight 
320,000, at 1 per cent concentration in o-dichlorobenzene.”® Very low 
frequencies used (~10 cycles) there and in an earlier study”’ suggest, 
however, that even here, appreciable entangling of the molecules created 
a temporary network such as studied by Ferry.””’ * Such was certainly 
present in the 5 to 18 per cent solutions of cellulose acetate in dioxane 
measured in one of the earliest observations of shear rigidity in polymer 
solutions.” | 

Accordingly, since strictly linear, and hence non-interacting, mechanics 
are sought for the macromolecules in dilute solution, careful evaluation 
of experiments is essential. Since already it appears that important 
over-all (quasi-configurational) relaxations occur for, say, polyisobutyl- 
ene in the kilocycle range, and it is suspected that not all of the inter- 
actions involved are between chains, the torsional crystal techniques 
are attractive. The absolute viscosity of these solutions is very low, so 
the ammonium dihydrogen phosphate crystal whose piezoelectric qual- 
ities are appropriate for polymer liquids in the circuits previously 
noted”: ” is advantageously replaced by quartz. 

Detailed electromechanical behaviour of such crystals in the pure 
liquids cyclohexane and benzene is of first concern. The electric field 
applied to electrodes on the suspended crystal produces mechanical 
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torsion generating pure shear waves. These waves may be modified by 
the environment around the crystal (vacuum, gas, liquid, solid) and 
react back. Thereby a mechanical resistanee, Ry, and a mechanical 
reactance, X, are imposed on the electrical properties of the crystal 
element in the circuit. This connection comes out as: 


AR» = KiR x 
Af = -K.Xu, 


where Al, is the increase in measured electrical resistance of the crystal 
element in the medium compared to in vacuum (or practically in dry 
air or nitrogen). The decrease in resonant frequency of the crystal ele- 
ment under these conditions is Af. Thus Ky and K, are electromechanical 
constants, which fundamentally may be calculated from the dimensions 
and piezoelectric constants of the crystals.’ Now, in simple, Newtonian 
liquids, 


Ru = Xu = Vaxfnp 


Thus, by carefully measuring ARz (or Af) on a liquid of accurately 
known density p and viscosity 7, at a given frequency f, and a given 
temperature, the constants K, and K. may be evaluated without as- 
sumptions and approximations of deriving them. Their constancy will 
then reflect the electromechanical stability of the system. Their be- 
haviour under various conditions will be illustrated below. One further 
point is that when a liquid or solution does exhibit shear rigidity, or, in 
other words, if the single large molecules in a dilute solution are able 
to store energy, then Ry» > Xs. Hence, in this case, the observed 
quantities ARz , and especially Af require particular precision. 

In this regard, typical magnitudes of change of fp between dry air 
and pure cyclohexane, at various temperatures, appear in Fig. 10. 
Questions often arise as to the arbitrariness of suspension of the radi- 
ating crystal, by the fine supporting and lead wires. The effects with 
the plain wires, in the solid curves of Fig. 10, are somewhat, but not 
radically, changed when a metal bead is put on, heavily loading vibra- 
tions in the wires, as shown by the dashed curves. In Fig. 11, a some- 
what larger influence of the loaded support wires is shown for the Rx 
values, but both curves, by their smoothness and shape over a tem- 
perature range where the thermal expansion and other elastic constants 
of the metal support wires are quite different from those of the quartz 
crystal, affirm reliability of mounting and electromechanical coupling. 

Fig. 10 shows, even for an 80-ke crystal, that Af for an organic liquid 
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Fig. 10—Temperature variation of resonant frequency before and after add- 
ing weights to mounting wires. 


°o 
Ww 


(or dilute polymer solution) is bothersomely small. An excellent oscillator 
at 20 ke can hardly be expected to drift less than +2 cycles, but at 20 
ke the Af like that between the sets of curves on Fig. 10 might be only 
10 cycles, so 20 to 35 per cent error could come in. Hence, a different 
scheme for measurement of fg than that in earlier systems” was 
evolved. The tenth harmonic of the (say 80 kc) resonant frequency was 
beat against the 79th harmonic of a controlled standard 10-ke frequency. 
An interpolation oscillator accurately readable to 1 cycle then supplies 


the many hundred (roughly 1000) difference between these two high 
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Fig. 11—Temperature variation of resistance at resonance before and after 
adding weights to mounting wires. 
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Fig. 12—Temperature variation of crystal constants K, and Ke at 80 ke be- 
fore and after adding weights to mounting wires. 


harmonics. In this way, and in about 30 sec a balance can be conven- 
iently achieved and the required ten-fold gain in accuracy attained. 

By these means, and with best literature values of viscosity and 
density (which were checked in the laboratory at several temperatures) 
for purified solvents, curves for K, and Ke were obtained as exhibited 
in Fig. 12 for 80 ke. Behaviour of K, at different frequencies over a tem- 
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Fig. 13—Temperature variation of crystal constant K, over a frequency range 
with benzene and cyclohexane as standard fluids. 
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perature range is shown in Fig. 13, and of K2, in Fig. 14. Fig. 14 brings 
out the significant point that in the present arrangement, where the 
oscillating crystal is immersed in the liquid studied, the dielectric prop- 
erties of the liquid are important. Apparently the dielectric losses even 
of these purified hydrocarbons are different enough so that K, at 80 ke 
is quite different for benzene and cyclohexane. (Dielectric studies have 
previously indicated difficulty in preparing benzene having theoretically 
expected loss.) It is also possible that slight differences in wetting the 
crystal cause Ky to vary with the liquid used. 

The K, and K» values determined for all the various conditions above 
were then used under these conditions for measurements on the polymer 
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Fig. 14—Temperature variation of crystal constant K»2 over a frequency range 
with benzene and cyclohexane as standard fluids. 
solutions in the kilocycle range. In the megacycle range, the balanced 
shear wave reflectance technique” gave satisfactory results over certain 
concentration zones which could fairly well be extrapolated to high 
dilutions. Thus, over the whole spectrum, there seems to be no doubt 
about the reality of the effects described below. That is, their magnitude 
far exceeds experimental uncertainty, as demonstrated in this section. 


POLYISOBUTYLENE SOLUTIONS; DYNAMICS OF SEPARATE CHAINS 


Solutions of polyisobutylene of M, = 1.2 X 10° from about 0.1 to 
1.0 wt. per cent concentration in cyclohexane yield Ry and X y curves 
as shown in Fig. 15. The points coincide for the pure solvent, as they 
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Fig. 15—Electromechanical interaction of solutions of polyisobutylene (M, = 
1.18 X 10°) in cyclohexane with crystal vibrating torsionally at 20 ke. 


should for a liquid having only viscosity. But, apparently as soon as any 
polymer chains are added, the curves diverge. A stiffness coming from 
separate chain molecules is being displayed.” Qualitatively, theoretical 
expectations of Kuhn”'” and others seem justified, at least that there 
is a relaxation mechanism for isolated chains. 

The usual question of how best to express the dynamical results arises. 
The procedure of earlier sections for polymer solids and liquids will be 
followed. In general, a frequency dependent modified Maxwell element 
as sketched on Fig. 16 will be used. However, a frequency-independent 
analysis has also been carried out for one sample system, and, from this, 
basic mechanical constants of single “‘average’” molecules are obtained, 
if it is reasonable to relate the mechanical models for the liquid con- 
tinuum to the discrete chains dissolved in it. 

Fig. 17 shows typical results from the simple scheme of Fig. 16, where 
the pure solvent viscosity, 74, has been considered to be in parallel 
with a Maxwell element. The total shear rigidity of the solution (at a 
given concentration) is represented by yz . The viscosity of the polymer 
molecule coils in solution with the solvent streaming through them is 
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' Fig. 16—Relations for calculation of shear stiffness and viscosity of dilute 
polymer solutions. 
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taken to be 7; . Thus, the steady flow viscosity, 7» = na + ns. Also, 


+ 
fe Nr or 2 = Nap 
NA NA 


under steady flow, or, alternatively, approximately a ‘dynamic intrinsic 


viscosity” 
[= | 
Na  C_jc-+0 


can be written for any given frequency. 

The curves in Fig. 17 are frequency dependent, however, although it 
turns out that 72 is only slightly so. Nevertheless, the considerable rise 
of n above the pure solvent viscosity, as the concentration is increased, 
indicates other mechanisms are being lumped into y4. As usual, some 
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Fig. 17—Rigidity and viscosities of polyisobutylene (M, = 1.18 X 10) in cy- 
clohexane, at 25°C and 20 ke. 
extensive distribution of relaxation times is probably responsible. How- 
ever, from the chemical point of view, it is best to see if some principal 
mechanisms related to known structures can be identified. If so, they 
could be associated with new ideas about the details of polymer intrinsic 
viscosities, as well as the form of isolated molecules.” ”* *” 

First, the frequency dependence of the yz of the model of Fig. 16 is 
as shown on Fig. 18. Striking regions of dispersion appear, although 
more points are needed to define the 10° cycle zone. Actually, many sets 
of data have been obtained in the 10° cycle zone. Recently, an immersed 
quartz tuning fork has given the approximate value shown for 2300 
cycles. The experiments of Fig. 18 were on a polyisobutylene having 
M, = 3.9 X 10°, dissolved in cyclohexane. Values of 74 and 1g were, 
of course, also obtained. The results were then analyzed for a system of 
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Fig. 18—Frequency dependence of shear stiffness of 1 per cent solutions of 
polyisobutylene in cyclohexane at 25°C. 
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Fig. 19—Schematic diagram of possible sources of rigidity of single chain 
molecules in solution. 
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three Maxwell elements in parallel with again, as in Fig. 16, a solvent 
viscosity, this time called m (truly absolute solvent viscosity) in parallel 
with them. The ys curve of Fig. 18, running through the observed 
points could then be calculated, by some special trial methods, with 
which Messrs. H. T. O’Neil and O. J. Zobel of Bell Telephone Labora- 
tories kindly helped. 

This analysis, identifying three principal relaxation regions for the 
motions of polyisobutylene chains in cyclohexane, gave for a 1 per cent 
solution (taken as linear part of concentration curve and hence equiva- 
lent to high dilution). 


Principal rigidities uw. = 890 dynes/em” 
3,190 dynes/cm” 
84,000 dynes/em? 


B3 


I 


ba 
Principal viscosities 7, = 0.0082 poise (pure cyclohexane) 
n. = 0.255 poise 
ns = 0.006 poise 
ns = 0.004 poise 
Principal relaxation frequencies fo = 550 cycles 
fs = 8.45 X 10° cycles 
fs = 3.52 X 10° cycles 


Tentatively, these mechanisms may be schematically described as on 
Fig. 19. Here, the polymer coil, subjected to shear waves in dilute solu- 
tion, exhibits rigidities uw. , ws and py, all shown on different scales. ye is 
the configurational elasticity because of actual changes in root mean 
square separation of chain ends, as from FR to R’. It is retarded by vis- 
cous drag through the solvent, 72, which is presumably the main source 
of characteristically high y, of chain polymer solutions. The relaxation 
frequency for this mechanism is low—a few hundred cycles. It may 
come in significantly in work on more concentrated solutions at low 
frequencies,’ *’ *: °’ where chain entanglement is nevertheless the dom- 
inant factor. 

#3 is When segments of the same chain in the molecular coil temporarily 
entangle with each other. Striking evidence has recently been given by 
Fox and Flory” that because of mutual interference, the theoretical 
random flight configuration of a chain gives very much too small a 
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molecular coil volume, V,. This suggests that thermal agitation tending 
toward a smaller V. , and excluded volume or repulsions forcing a large 
one, will cause collisons or entanglements which might last long enough 
to give a van der Waals cross-bond as denoted by crosses on the ps 
sketch. (The actual forces in these would somewhat resemble those be- 
tween different molecules in the concentrated solutions of Ferry.” *’) 
This mechanism has the reasonable (based on Ferry’s and others’ work) 
relaxation frequency of 8.45 X 10°. A small viscosity, 73 , may comprise 
friction of slippage at the entanglement points, with both the polymer 
and associated solvent molecules. 

In Fig. 19, us is a relatively high stiffness presumed to be some average 
hindrance to rotation of one segment with respect to another. In the 
sketch, close-packed spheres. representing methyl groups in polyiso- 
butylene are portrayed. Their force fields overlap more in some places 
than others, in the meandering of the chain to form the molecular coil 
(of course, some tail-to-tail structures may be important here; they have 
all been shown head-to-tail in the sketch). Thus, this total internal steric 
restraint on chain flexibility, with a relaxation frequency of 3.5 X 10°, 
contributes greatly to the large dispersion of rigidity in the megacycle 
range noted in Fig. 18. The related viscosity, 4 , is again low. 

There is no doubt a considerable distribution of relaxation character- 
istics associated with each and all of these mechanisms. 


Physical Properties Per Molecule 


Since the viscosities and rigidities in the dilute solutions indeed seem 
to be additive with the number of molecules present, values of these 
properties, for the hypothetical mechanisms, can be expressed per aver- 
age chain. Of course, the measured quantities are expressed as constants 
per cc of solution, but it may be useful to think of in terms of one aver- 
age chain in each cc. Then, the shear deformation of this chain could be 
denoted by a force constant. The associated viscosities remain, however, 
dependent on solvent surroundings. Thus, for the polyisobutylene of 
M, = 3.9 X 10°, in cyclohexane solution, at 25°C the molecular quan- 
tities are: 


[fe] = 17 X 107" dyne cm [nm] = 1.6 X 107° poise 
[f] = 6 X 10°” dyne cm [ns] = 3.9 X 10~* poise 
[fi] = 16 X 10°" dyne cm [ns] = 2.4 X 10° poise 


In the section on polymer liquids, the high-frequency modulus pu was 
attributed to a nearest-neighbor glass or crystal-like interaction (since 
the actual values were indeed typical of the hardest organic solids). 
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However, in polyisobutylene (and to some degree in poly-a-methyl 
styrene), it is especially difficult to distinguish inter-chain from intra- 
chain crowding of methyl groups. Thus, while average center-to-center 
separation of methyls is ~4 A in adjacent chains,”° it is <2.5 A within 
chains, in polyisobutylene. This crowding is apparently strong; the ob- 
served AH pn is only 12.8 kcal per mole instead of the 19.2 expected.” ”° 
The energy of steric hindrance thus amounts to almost half of the actual 
heat of polymerization. It is reasonable that a large part of the hardness 
of a mass of polyisobutylene chains, such as in the liquids, should there- 
fore reflect the same mechanism as that for ys (Fig. 19) in the dilute 
solutions. A rough check on this can be made. A polyisobutylene having 
considerably lower molecular weight than 3.9 X 10° and thus inter- 
mediate between the “liquid” and “‘solid” ranges, had a Maxwell shear 
modulus in the megacycle region (14 me) of » = 5.3 X 10°, at 25°C. 
The number of molecules/cc, with individual [f,] given above, necessary 
to give the observed density of this polymer was multiplied by [f:], giving 
wu = 2.8 X 10° dynes/cm’. Accordingly, about half of the observed high 
frequency rigidity of polyisobutylene, at 25°C, may be calculated from 
a ‘‘molecular constant’? embodying intra-chain stiffness. 

Much more refined and detailed treatments are required to generalize 
these ‘‘molecular constants”? which are after all, as shown below, de- 
pendent on using a thermodynamically “inert” solvent. However, much 
as structurally significant dipole moments can be derived from measure- 
ments in dilute solutions, it seems hopeful that macromolecular me- 
chanics can be so elucidated. Also additional structures, such as poly- 
propylene and polydimethy] siloxane compared to polyisobutylene, are 
currently being studied. 


Temperature Variation 


Some further behaviour at different temperatures and solubilities of 
separate chains in dilute solution may now be considered against this 
background of possible mechanisms. Practically, these studies will bear 
on processing and properties, lacquers, paints, and casting solutions of 
polymers, as well as on the other qualities outlined in the introduction. 
Results may be conveniently discussed in terms of the modified Max- 
well single element, with factors 4, ns, and uz (Fig. 16). Mostly, the 
kilocycle range, reflecting molecular coil changes, will be of interest. For 
comparison, it may be noted that at 20 ke, the polyisobutylene whose 
us = 1061 dynes/cm’ in 1 per cent solution in cyclohexane receives 889 
dynes/em’ of this from pz , the retarded configurational mechanism; 169 
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Fig. 20—Temperature variation of rigidity of 1 per cent solution of polyiso- 
butylene (M, = 1.18 X 105) in cyclohexane. 


dynes/cm’ from the entangled segments stiffness, 3 ; and only 3 dynes/ 
em’ from the intra-chain stiffness, us. Thus, chain configuration me- 
chanics, including the associated viscosities, can be well enough thought 
of in the following paragraphs, in terms of us, ne and na. 

The exponential decrease of » with temperature familiar for polymer 
solids and liquids is much suppressed in the yes vs JT curve of Fig. 20. 
While the ys, internal rotation, mechanism for single polyisobutylene 
molecules probably has a considerable activation energy, that for the 
we, configurational, rigidity should be very small. Then, without re- 
tardation, the intrinsic chain modulus would rise with rising tempera-. 
ture. These influences seem to combine to give the modest decline of ys 
appearing in Fig. 20. If these rigidities are plotted against 1/7, the 
temperature coefficient is 2.3 kcal. This is much less than the familiar 
values for the stiffening of rubbery solids, and emphasizes that inter- 
chain action reigns then. 


Solvent Variation 


Effects of solvents of different (mostly positive) heats of mixing on 
state of polyisobutylene molecules in solution have been nicely estab- 
lished by Fox and Flory.” Especially, this work has clarified principal 
factors in the intrinsic viscosity expression 


[n] = —e oe = KoM "08g. 
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Here, V. = effective volume per molecule (and hence as determined by 
chain configuration), M = molecular weight, a represents change in 
linear extent of molecule because of mutual interference of segments” 
and ¢ expresses the hydrodynamics interaction of solvent and molecular 
coil (including varying degrees of “straining through” the coil).°" 
Interpretation of the mechanical properties of chains in dilute solution, 
with reference to the rough concepts of Fig. 16, arouses particular in- 
terest in the factor a*. For a high molecular weight polyisobutylene, 
intrinsic viscosity theory” indicated that a* the ratio for volume of 
actual coil divided by volume for ideal random flight coil was 3.81 in 
cyclohexane but only 1.42 in benzene, both at 30°C. This striking alter- 
ation in equilibrium chain configuration, a variable which is not readily 
introduced into polymer liquids or solids, appears in the inherent vis- 
cosity vs c curves in cyclohexane, Fig. 21, and benzene, Fig. 22. The 
large difference in [7] at 25°C, 6.00 in cyclohexane vs ~1.5 in benzene, 
indeed emphasizes the different solvent powers.” Likewise, the large in- 
crease of [yn] with temperature in Fig. 22 accents the poor solvent qual- 
ities of benzene.*° Too, empirically, polymer molecules which are either 
tight coils or are actually chemically cross-linked to form microgel mole- 
cules characteristically show positive slopes of inherent viscosity vs ¢ 
plots.°” Accordingly, all this evidence for large changes in the conforma- 
tion of chain molecules in ‘‘good” vs ‘‘poor’’ solvents should show up 
in dynamics of dilute solutions. Also, technically, quite different physical 
properties are found for polymer-plasticizer compounds where compat- 
ibility is high (good solvent) than where it is low (poor solvent). Here, 





CONCENTRATION IN GRAMS PER !00 ML. OF SOLUTION 


Fig. 21—Inherent viscosity of polyisobutylene (M, = 3.87 X 10%) in cyclo- 
hexane at 25°C. . 
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Fig. 22—Inherent viscosity of polyisobutylene (Jf, = 3.87 X 10°) in benzene, 
at various temperatures. 
more flexible compositions are often produced with low compatibility 
plasticizers—indeed, sometimes with those on the verge of phase sepa- 
ration than with those with highly favorable heats of solution.” This 
would mean that the bad solvents would compress the chains so that 
they would be more easily strained than if they were in a “free chain” 
or even extended configuration. If single chain, visco-elastic stiffnesses 
are acting this way, the dynamic ys would then actually decline as heat 
of mixing become more positive. 
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Fig. 23—Rigidity of polyisobutylenes in cyclohexane and benzene at 25°C and 
C. 
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This seems to take place, as indicated by the lower compared to the 
upper curve on Fig, 23. Here, the ws of the usual modified Maxwell 
model, at 20 ke, is plotted against c for the polyisobutylene of M, = 
3.9 X 10°. Also, the middle curve shows ys for a polymer of about a 
third of this molecular weight; while there is a small reduction in ys 
with M, in this range, it is much less than the reduction caused by 
tightening up the polymer coil. 

The ys values per average molecule, [fs], fall from 18 X 107° dyne 
cm in cyclohexane to 7 X 10° in benzene. (Of course, [fs] for the inter- 
mediate molecular weight polymer in cyclohexane is only 5 X 10° be- 
cause so many more molecules are present in solution.) 

The temperature dependence of yz also becomes nearly zcro at least 
over the narrow range from 25 to 50°C, in benzene compared to cyclo- 
hexane. This seems to accord with the indications previously, from Fig. 
20, for a lower molecular weight polymer, that different mechanisms are 
competing. These may be the configurational, with » « J, and relaxa- 
tion, with » varying in some complicated way with 7’. Thus [7] increases 
markedly with 7’, and presumably denotes an expanding molecular coil 
tending toward the ‘“‘normal”’ configuration in cyclohexane. At the same 
time, the relaxation processes with rising temperature tend to cause the 
decrease in yz typical of the upper, solid, curves on Fig. 24. In engineer- 
ing use, often times poorly compatible plasticizers give compounds 
- which stiffen more gradually with temperature than do “solvent’’ plas- 
ticized ones. 

For similar reasons, the dynamic molecular coil viscosity, 7: , ought 
to vary less with temperature in thermodynamically poor than in good 
solvents. This is indeed seen in Fig. 25. On the other hand, y4 for the 
modified Maxwell element has been described as the solvent viscosity 
with segment hindrance and restricted rotation terms from the polymer 
molecules lumped in with it. These latter terms are presumably little 
affected by over-all configuration (u2 term; the »; mechanism will be 
somewhat affected, but not the us, on Fig. 19). Thus, 74 should have 
comparable temperature dependence in both good and bad solvents, as 
seems to be indicated by Fig. 26. 


Microgel Molecule Solutions 


The statistical coil of linear polymer molecules may be replaced by a 
chemically fixed, cross-linked network in microgel molecules.” These 
may be made completely rigid, like Einstein spheres, or highly swell- 
able. The latter are hybrids between rigid spheres and coiled chains. In 
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Fig. 24—Temperature variation of rigidity of 1 per cent solution of polyiso- 
butylene (M, = 3.87 X 108) in cyclohexane and benzene. 
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Fig. 25—Temperature variation of yz for 1 per cent solution of polyisobutylene 
(M, = 3.87 X 10°) in cyclohexane and benzene. 
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synthetic rubber, they confer unique flow properties, causing the excel- 
lent processibility of GR-S 60. However, dynamic tenacity, such as in 
flex crack growth, is degraded by their presence. Now presumably the 
excellent extrusion qualities of synthetic rubber composed of from 60 to 
80 per cent microgel molecules are because of their individual shear 
stiffness. Thus, if a wire coating, for instance, is extruded at high rates 
of shear, chain molecules are deformed, and store energy just as dis- 
cussed in the earlier sections on liquids. After emerging from the extru- 
sion die, they relax, and cause the gross retraction, shrinkage and rough- 
ness shown in the wire insulation of the upper photograph of Fig. 27. A 
polymer with about 70 per cent microgel molecules gives the smooth 
covering shown in the lower specimen of Fig. 27. Here, the shearing 
stresses of extrusion seem insufficient to distort the tiny networks of 
the microgel molecule; in any case, the covering does not roughen or 
relax. Similar effects have been found for microgel plastics. Neverthe- 
less, unlike gross or macro gelation, the whole melt can flow. 

On this basis, dilute solutions of microgel molecules ought to indicate 
high shear rigidity per molecule. The mechanism y; of Fig. 19, in which 
now the junction points are not temporary, but are primary valence 
cross-links, should be predominant. Fig. 28 shows, for a polybutadiene 
microgel in cyclohexane,” that ys has indeed risen, compared to equal 
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_ Fig. 26—Temperature variation of n4 for 1 per cent solution of polyisobutylene 
(Ma = 3.87 X 10°) in cyclohexane and benzene. 
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Fig. 27—Effect of microgel molecules in synthetic rubber on smoothness of 
extruded wire insulation. Rough covering is from high-speed extrusion of GR-S 
without microgel. 


weights of chain molecules. Further, accompanying the extremely high 
average molecular weight of the microgel (18.6 X 10°), the [fs] per 
average molecule is 42 X 10” dyne cm or about twenty-five times that 
of the polyisobutylene with M, = 3.9 X 10°. Also, the temperature 
coefficient for uz of polybutadiene microgel is low. 

Of course, polybutadiene; as chains or as microgel molecule segments, 
has many double bonds. These will surely influence the ps, or internal 
rotation mechanism. Further work remains to show just what is their 
effect in the microgel case. But, it is interesting to compare us values 
for Hevea rubber chains with those for, say, polyisobutylene, which has 
only single bonds in the chain. 
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Fig. 28—Rigidity of 0.5 per cent solution in cyclohexane of polybutadiene 
microgel (Mw = 18.6 X 10°) at 20ke. 
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Hevea Rubber Solutions 


The comparison of equal weight concentrations of natural rubber in 
cyclohexane with polyisobutylene in cyclohexane is surprising: 


Hevea rubber M, = .23 X 10° ys = 1350 dynes/em’, 1 per cent 
solution (corr.). 


Polyisobutylene M, = 1.2 X 10° us = 1000 dynes/cm’, 1 per cent 
solution (corr.). 


Both results are at 20 ke. The higher value for natural rubber may be 
because of the double bonds causing stiffening of the chain. On the other 
hand, maybe easy rotation around single bonds raises the u3 part. Cer- 
tainly the viscous retardation within natural rubber chains is very low, 
as noted in the section on solids. However, its interaction with, or con- 
figuration, in cyclohexane may be peculiar. The [fz] per average mole- 
cule is, however, low, being 15 X 107“* dyne cm at 25°C. 


Polystyrene Solutions 


Much work, on light scattering and other properties, has indicated 
appreciable intra-chain stiffness for polystyrene,” but still much freedom 
compared to polyisobutylene.”* However, this work, as well as AH pin of 
17 kcal compared to ~19 kcal calculated for no steric hindrance, sug- 
gests comparatively small restraints on ideal flexibility. This needs to be 
checked by a frequency analysis of dilute solution mechanics, but poly- 
styrene seems to be a reasonable example of ‘plastic’ behaviour at room 
temperature because of interaction between the chains. (It is recalled 
that, earlier, a-methyl styrene polymer was cited as plastic model show- 
ing both intra- and inter-chain stiffness. Unlike in polystyrene, the 
intra-chain factor shows up in a AH,,, of 9-10 keal, a third less than 
that calculated if there were no steric hindrance.) Thus, no evidence of 
unusual stiffness appears in Fig. 29, when, indeed, the us values are 
considerably lower, for equal weight concentrations, than those for nat- 
ural rubber. The highly milled rubber studied had M, very nearly that 
of M, = 0.234 X 10° of the polystyrene, so the [fs] per average poly- 
styrene chain, 4.5 X 10 “ dyne em is less than a third that of the rub- 
ber. No wonder that at high temperatures, where the phenyl group 
interaction between chains is much reduced, polystyrene makes a good 
rubber. Also, in Fig. 29 are shown data for a polymer of M, = 1.2 X 10°, 
made in emulsion and having [y] = 4.850 in benzene at 25°C. 
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The polystyrene solutions discussed above were in benzene, a good 
solvent. Here, the situation is converse to that for polyisobutylene; for 
polystyrene, cyclohexane is a poor solvent and benzene, good. Hence, if 
the previous interpretation of reduced single chain quas?-configurational 
(ue) stiffness is general for solvents of more endothermic mixing, the 
“plastic” molecule polystyrene should show it in cyclohexane. This is 
indeed evident in Fig. 30, showing one of the same polystyrenes of Fig. 
29, measured at 20 ke (normalized to 1 per cent concentration). Also, on 
Fig. 30 are shown the inherent viscosity (practically, the intrinsic vis- 
cosity, in this case) and the absolute viscosity of the 1 per cent solution 


A Mp = 234,000 
Oo My =1,204,000 










is 4 | 
A pe. —_| 
x | 
0) 10 20 30 40 50 60 70 80 90 {0 
FREQUENCY IN KILOCYCLES PER SECOND 


Fig. 29—Change of shear stiffness, uz , with frequency, for 1 per cent solutions 
of polystyrene in benzene. 
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under steady flow, 7; . These are all plotted against temperature down 
to phase separation, at about 26° to 27°C. 

The marked positive slope of the fy,/c curve denotes the large con- 
traction in molecular coil volume preceding phase separation or insolu- 
bility. The absolute viscosity, y;, however, rises with declining tem- 
perature because it is dominated by solvent viscosity, but when the 
polymer phase comes out, 7; abruptly falls off. 

The us values are consistent with this steady flow behaviour, except 
that the rise of us at the turbidity point seems to be because a layer of 
swollen polymer-rich phase forms on the torsional crystal surface. This 
condition is seen in Fig. 30 to coincide nicely with the abrupt changes 
in steady flow viscosity. 

The slight maximum in the ys curve at about 35°C may not be real. 
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It does come near the point of minimum interaction for the whole sys- 
tem. In any case, as discussed before, the average temperature coefficient 
of wz in the poor solvent is very low compared to the good solvent. The 
values of uz are roughly 2 to 4 those in benzene. 


GENERAL THEORY OF SINGLE CHAIN MECHANICS; KUHN AND KIRKWOOD 


As noted before, much of the present understanding of stress-strain 
properties of polymer chains, in dilute solutions, liquids or solids, has 
come from W. Kuhn’s long interest in this subject. Many supplementary 
contributions have been stimulated by Kuhn’s work, and new points of 
view have been introduced by others. For instance, recently new and 
different proposals have been made about the flow birefringence and 
non-Newtonian viscosity of solutions of deformable spheres.’ These 
ideas could be tested on suitable microgel solutions. 

Recently, moreover, an especially significant general theory of visco- 
elastic behaviour of polymer in solution has been constructed by Kurk- 
wood.” It explicitly considers the hydrodynamic conditions leading to 
the rigidity now observed for high-frequency shear waves. It formulates 
definitely the configurational changes of isolated chains in solution when 
strained in shear. As this theory is advanced to forms where simpler 
calculations can be made, it may answer many of the questions raised 
by the new experiments on single chain properties. 
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Fig. 30—Temperature dependence of absolute viscosity, 7g , inherent viscosity, 
énnr/c, and shear stiffness, uz of, 1 per cent solution of polystyrene in cyclohexane 
down through turbidity point. 
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CONCLUSIONS 


To leave some impression of the elemental chemical structures which 
move around when wood, rubber, plastics, textiles and finishes are used 
mechanically—that has been the aim of this study. Polymer viscosities 
have been found in a variety of “solids’’; rigidities have been demon- 
strated for very fluid “liquids” and solutions. Studies of these solid and 
liquid extremes have given some chemical reality to the classical spring 
and dashpot models. 

Existence of compressional viscosity has been shown for polymer 
liquids and solids. It may comprise a new quality for investigation of 
polymer structure. At present, too little is known of its origin to inter- 
pret further the effects of intense ultrasonic irradiation of polymer solu- 
tions. Experiments of Schmid and co-workers” early indicated degrada- 
tion of molecular weight of polystyrene, so irradiated, but whether this 
is chemical, from local heating in the solvent, or actual physical coupling 
with the wave field, is still unsettled. However, these workers also con- 
sidered a compressional stiffness of the polymer molecules in the solu- 
tions,” and showed that if there was coupling, it was not inertial (by 
dissolving polystyrene in solvents of exactly the same density, no reduc- 
tion in effect was observed). A point of general interest arises here; 
impact fractures of plastics presumably actually fracture some primary 
valence bonds. This is certainly true for many thermoset materials, and 
probably for chain compounds. Hence, if the detailed mechanism of how 
compressional waves move and perhaps rupture polymer segments were 
known, information on the baffling problems of ultimate strength would 
be gained. The observations above on dependence of \ and \’ on molec- 
ular weight and structure provide only the barest start on this but a 
new goal is in view. Too, basic questions of how rapidly molecules being 
formed in a polymerization equilibrate in temperature with their sur- 
roundings are elucidated by compressional wave propagation constants. 
For instance, absolute rate measurements on velocity of chain growth 
cannot be said to be isothermal if they seem to be faster than the ther- 
mal relaxation times which the ultrasonic measurements indicate can be | 
~10~° to 10 ” sec. 

Likewise, more thorough understanding of velocity and dispersion of 
compressional waves in polymer solutions would clear up anomalies in 
velocity measurements for a wide variety of polymers,” some of which 
have been tentatively attributed to chain branching. 
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The Reliability of Telephone Traffic Load 
Measurements by Switch Counts 


BY W. S. HAYWARD, JR. 
(Manuscript received October 15, 1951) 


The switch count method of telephone traffic measurement is subject to sam- 
pling errors. The nature of these errors is discussed and formulas are de- 
rived which describe the extent of the errors under normally encountered 
traffic conditions. 


INTRODUCTION 


Of prime importance to the telephone traffic engineer is the deter- 
mination of the busy season busy hour load carried by groups of trunks 
or other circuits of a telephone switching system. Three direct methods 
of measuring such loads are found in the field today. These are: 


a. Peg Count and Holding Time Method 


The number of calls carried by the circuit group during the observa- 
tion period is counted. This number multiplied by the average holding 
time per call (in hundreds of seconds) and divided by the length of the 
observation period (in hours) gives an estimate of the group load in 
units of hundred-call-seconds per hour (CCS). The major drawback to 
this peg count method is that it requires a separate determination of the 
average holding time per call for the group under observation. R. I. 
Wilkinson! has analyzed the sources of errors of holding time measure- 
ments. In addition, correlation between load and holding time introduces 
an error which has not been studied. 


b. Switch Count Method 


At fixed intervals the circuit group is scanned and the number of busy 
circuits is counted. The total number of busy conditions counted divided 
by the number of scans is, then, an estimate of the load on the group in 
units of average simultaneous calls or erlangs*. This estimate is generally 
converted to CCS (1 erlang = 36 CCS) by traffic engineers since the 


* The name ‘“‘erlang’’ for average simultaneous call was adopted at a plenary 
meeting of the CCIF at Montreux in October, 1946. 
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load entries of most traffic tables are in terms of CCS. For theoretical 
studies the erlang is a more convenient unit and will be retained here. 


c. Continuous Method 


The busy condition of each circuit is represented by a fixed increment 
of electrical current through an ampere-hour meter. The instantaneous 
current is then analogous to the calls simultaneously present so that the 
meter, which integrates the current, may be calibrated to indicate hun- 
dred-call-seconds or erlang-hours directly. Although this method is po- 
tentially the most accurate, practical difficulties have limited its use. 


In addition to these direct methods, there are several methods of in- 
direct load measurement which, relying more heavily on traffic theory, 
make use of partial load indications, such as duration of group busy or 
the number of calls finding the group busy. Such measurements are less 
reliable than the direct measurements particularly when applied to un- 
derloaded groups. 

This paper is concerned with the reliability of switch count load meas- 
urements since this method appears to have prospects of considerably 
wider adoption in the future. Main emphasis will be placed, both quali- 
tatively and by the application of error formulas, on the relative effects 
of various measurement and traffic parameters on the accuracy of switch 
count measurements. Where long derivations of formulas are required 
they are deferred to the Appendix. 


SOURCES OF ERROR 


As has been described, switch count measurements yield the average 
number of calls found present when a group of circuits is scanned at 
fixed intervals during an observation period. Usually only that period 
of the day during which the load is greatest is of interest to the traffic 
engineer. Because the load during such periods also fluctuates from day 
to day, measurements of the loads for several days must be averaged to 
provide a useful load estimate. 

There are two main sources of error, therefore, in switch count esti- 
mates of telephone traffic loads: 

1. Each individual count of busy circuits is separated from the next 
by a time interval during which changes in load are not detected. Con- 
sequently, the load indicated by measurement may differ appreciably 
from the actual load carried. This difference can be decreased by de- 
creasing the interval between scans. 
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2. Even if the load carried during a measurement period were known 
very accurately, it is still only a sample of the many loads that might be 
offered by the same source of traffic under statistically identical con- 
ditions. Therefore, the average of several load readings may be expected 
to be somewhat in error as an estimate of the true average of the traffic 
source. The latter will be referred to as the source load to distinguish it 
from the carried load. 

Mechanical and human errors are likely to be present as well but, since 
' they are not inherent in the switch count method, they will be neglected 
here. 


SWITCH COUNT ERROR 


As shown in the Appendix, for periods of observation which are rela- 
tively long with respect to average holding time made on traffic with 
certain assumed characteristics, the average error of switch counts in 
estimating traffic load carried in the same period is zero. The coefficient 
of variation of the error, which is the standard deviation of the error 
expressed in per cent of the traffic load carried, is given by: 


Vz = 100 V/ | retan(5) — 2| —_ (1*) 


PX babe. : r 
etnh() ae ea ie hyperbolic cotangent of 5 


re = T/t > 20 





where r = ratio of scan interval to holding time 


oH 
I 


average holding time 
a’ = carried load in erlangs 
ec = number of switch counts 
T = length of observation period 
N = number of observation periods 


and where the following assumptions are made: 

a. Calls originate individually and collectively at random.t 

b. Holding times are exponentially distributed. 

c. Congestion loss from the group is negligible. 

* T have recently learned that these carried load formulas have been published 
by Conny Palm in Tekniska Medelanden fran Kungl. Telegrafstyrelsen, 1941. nr. 7-9. 


Tt See T. C. Fry, Probability and Its Engineering Uses, D. van Nostrand Co. 
Inc., New York, p. 216, for a definition of this condition. 
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As shown in the Appendix, this formula simplifies, when r < 2, to 


_ 100 i 
Yes c/T / aT 2) 


where c/T = rate of scan in cycles per time unit. From equation (2) it 
is apparent that if the scan interval is of the order of a holding time, the 
error of an estimate of traffic carried is inversely proportional to the rate 
of scan and inversely proportional to the square root of average load, 
holding time and hours of observation. For example, take the case where ~ 
switch counts are made during the busy hour, five minutes apart on a 
trunk group carrying calls with an average holding time of 3 minutes 
and an average load of 5 erlangs (180 CCS). What is the error in the 
estimated load carried if the readings for ten days are averaged? (As- 
sume conditions (a), (b) and (c) are met.) We have 





N = 10 observation periods 


T = 1 hour 
é = 1/20 hour 
a’ = 5 erlangs 


° 
I 


12 scans per observation period ~ 
re = T/t = 20 average holding times per observation period 


From equation (2) since T/t = 20 andr = T/c = 1.7 


100 / 1 
ae 12 VY 6-5-10-1-1/20 21870 


If, as proposed in the Appendix, it is assumed that the error has a 
normal distribution, there is 90 per cent assurance that observed values 
will fall within 1.64V, , or in the example within 3.52 per cent, of the true 
average*. Note that this error limit would be halved if the rate of scan 
were doubled or if four times as many hours of observation were taken. 

The coefficient of variation of the switch count error for constant values 
of 7'/t as a function of r is plotted on Fig. 1 for one observation period 
of a one erlang load. For loads other than one erlang the coefficient of 
variation is found by dividing by +/a’N. Thus in the example we have, 
using the dotted curve, 


15 


* This assumes that a sufficient number of observations are taken so that a 
priori information may be neglected in making an estimate of the universe. 
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or more accurately using the solid curve, 
ae 23% 
== 5/510 ~ 287 


The error of using equation (2) is seen to be negligible for most purposes 
even when 7'/t is less than 20. The probability of an observation oc- 
curing within a given number of standard deviations is widely published 
for the normal curve. A few values are given below: 








Z Pz Probability of exceeding +: zo or + 2V 
0.6745 0.50 
1.44 { 0.85 
1.64 0.90 
2.00 0.9545 
3.00 0.9973 





Tig. 2 is a plot for 40 observations of measured load vs carried load. 
Hach observation was made for a half hour period on a panel line finder 
group with switch counts made at the start and middle of the period. 
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EXACT FORMULA 





Tr, RATIO OF SCAN INTERVAL TO AVERAGE HOLDING TIME 


; [r-ctnh (S) -2| (rc-Kr+2b) +(Kr2-2kr4 2b) 


ro 
0.1 : 
0.8 1.0 3 4 5 6 8 10 20 30 40 60 80100 200 
VefaN 5 COEFFICIENT OF VARIATION OF ERROR IN ESTIMATING A ONE 
bs “ERLANG CARRIED LOAD IN ONE OBSERVATION PERIOD 


. Fig. 1—Accuracy of switch count estimate of load actually carried. 
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@ =1 OBSERVATION PERIOD 
T =1/2 HOUR 

t = 3 MINUTES 

i =15 MINUTES 

C = 2 SCANS 


Oy=0.615Va FROM FIG. 1 


NK | a 
S 


nN) 
a 


tw 
nN 





nN 
(o} 








© 


vA 
7 


ig 

LA 
7a 
Lae 
Pes 
Age 
eae 


7 
7 
7 


SENS 


Q 


CARRIED LOAD IN ERLANGS ESTIMATED BY 2 SWITCH COUNTS 
R 








10) 2 4 6 8 10 12 14 16 18 20 22 24 26 28 
CARRIED LOAD IN ERLANGS ESTIMATED BY 60 SWITCH COUNTS 


Fig. 2—Accuracy of switch count estimate of true average load. 


This is compared with the average of switch counts made every 30 sec- 
onds which has a relatively negligible error. The average holding time 
per call for the group was 176 seconds. The accuracy of only two 
counts is surprisingly good and the observations are seen to lie satisfac- 
torily between the 2c limits. 


ERROR OF TRAFFIC IN A GIVEN PERIOD AS AN ESTIMATE OF THE SOURCE 
LOAD 


The average traffic carried in two different periods but generated by 
the same traffic source is subject to statistical variation. As a result, any 
measurement of load, even if measurement errors are eliminated, is only 
a sample of the wide range of traffic loads that might have been gen- 
erated by the same source of traffic under identical circumstances. 
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J. Riordan has shown? that the standard deviation of the average 
traffic load for any one period is given by 


oy = 4/7 (; at on eri) (3) 


where a = the average source load 








t = average holding time per call 


wh 


length of observation period 


(Assumptions are as before with an additional one that all periods are 
in statistical equilibrium) 
When t/T' <1 this reduces to the form also given by F. W. Rabe’ 


pets oe (4) 


or expressed in a per cent of the average 


2t 
y= — 5) 
Vy = 100 ae (5) 


When N periods of length 7 are observed the coefficient of variation is 


reduced further to: 
7 Ot 
Vy = 100 Va (6) 


In the example of the previous section, 
N = 10 
T=1 
t = 1/20 
a=5d 


_ [EAPO 


COMBINATION OF ERRORS 





Evidently if switch count readings are used to estimate the average 
which may be expected in other periods, the two errors described above 
should both be taken into account. The errors are probably correlated 
but this correlation is weak and at present no method of allowing for it 
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is evident. Such a refinement would probably change the equation for 
standard deviation only slightly from that derived for the independent 
case; therefore independence will be assumed. The standard deviation 
of the sum of two independent variables is the square root of the sum of 
the squares of the component standard deviations: 


Os = Vo + o, (7) 

r | fa’ 2ia 
= a ee pee dhe ee 8 
/ [reinh(5) 2| nT + aa (8) 


Assuming - is approximately unity, that is, that carried load is approxi- 





mately equal to source load, 


t r 
Veni ant = 
00 a retnh( 5) (9) 


In the example given, 
Vs; = 4,96% 


There is, then, 90 per cent assurance that the source average is within 
1.64 X 4.96 = 8.1 per cent of the observed average. Note that doubling 
the switch count rate (which halves the switch count error) reduces the 
total error only to 7.6 per cent (about 6.7 per cent improvement), while 
doubling the number of hours of observations reduces the error to 5.9 
per cent (about 30 per cent improvement). Plots of the coefficient of 
variation of a one hour observation of a one erlang load versus scan rate 
for various average holding times are given in Fig. 3 for a wide range of 
holding times. The coefficient of variation of error in estimating other 
loads may be found from Fig. 3 by dividing the unit load coefficient 
by VaNT. In the example, the unit load coefficient is found, by entering 
Fig. 3 with ¢ = 3 minutes and rate of scan = c/T' = 12/1 scan cycles per 
hour, to be 35.0 per cent. Dividing by +/5-1-10 gives a coefficient of 
variation of 4.96 per cent as before. It is evident from Fig. 3 that in- 
creasing scan rates is not a universal way to improve the accuracy of 
source load estimates. 





CHOICE OF SCAN RATES 


What then governs the choice of scan rate? Evidently increasing the 
rate increases the accuracy of carried load estimates to any point de- 
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sired. This is far from true if source load is being estimated. If the cost 
of making a scan is constant, increasing the number of observation 
periods and decreasing the scan rate will improve accuracy of source 
Joad estimates without changing measurement costs. The number of 
hours available for measuring, of course, limits this procedure, while 
the increase in accuracy becomes negligible as r becomes large. On the 
other hand, if the cost of each observation is only slightly affected by the 
cost of making additional scans, a high scan rate might be justified. 

In applying the above relationships to traffic measurements, the 
usual question raised by the traffic engineer will be either how many 
hours of data need he take to be reasonably sure of his estimate or, 
conversely, how sure is he of an estimate based on available data. 
Assuming as before that the error distribution is normal, the per cent 
plus or minus error limits within which a proportion, P, , of the estimates 
will fall is given by zV; ; the value of z corresponding to any selected 
P, may be found from tables of the normal probability distribution. 
“Reasonably sure” is often taken to mean that there is 90 per cent 
assurance that the error does not exceed 5 per cent. When P, is 0.90, 
zis 1.64, so that under this condition 1.64V, = 0.05, or V; = 0.0305. 
Given scan rate and holding time, V, is proportional to 1/+/aNT accord- 
ing to equation (9) or Figure 3. When J; is held constant, aNT is con- 
stant so that the plot of log NT against log a is linear, as shown in 
Figs. 4 and 5. The number of hours needed to meet any chosen ero 
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Fig. 3—Efficiency of switch counts for usage measurement. 
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Fig. 4—Hours of measurement required for 90 per cent assurance that error in 
estimating source load does not exceed plus or minus 5 per cent when measuring 
interval is 120 seconds. 


requirements may then be read directly from such graphs. In the second 
type of question, z, NT’, scan rate and holding time are fixed so that 
zV; is proportional to 1/+/a. Plotting log zV; against log ~/a again 
gives a linear plot as shown on Fig. 6. 

In the numerical example above, the limits of error corresponding to 
90 per cent assurance may be read from Fig. 6 which is plotted for the 
appropriate assurance, average holding time and scan interval. Reading 
the error limits at the point where the 10 hours measured line crosses 
180 CCS (5 erlangs) gives + 8.1 per cent as before. Fig. 5 may be 
entered to find the total number of hours required to reduce this error to 
5 per cent. Reading at the point where the 180 second holding time line 
crosses 180 CCS gives 26 hours. 


QUALITATIVE EXTENSION OF THEORETICAL APPROACH 


The original traffic assumptions made in deriving the theoretical re- 
sults above are: 
a. Calls originate collectively and individually at random. 
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Fig. 5—Hours of measurement required for 90 per cent assurance that error in 
estimating source load does not exceed plus or minus 5 per cent when measuring 
interval is 300 seconds. 
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Fig. 6—Limits of error reached with 90 per cent assurance in estimating source 
oad. 
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b. Holding times are exponentially distributed. 

c. Congestion loss from the group is negligible. 

d. Observation periods are in statistical equilibrium. 

How do departures from these assumptions affect the reliability of 
usage measurements? 


a. Holding Time Distribution 


Experience in application of delay and loss formulas has shown that 
theories based on exponential holding times are often applicable to other 
holding time distribution cases-which have a wide range. However, for a 
constant holding time distribution special theories often are called for. 
The average and standard deviation of switch count estimates of carried 
load when holding time is constant, are given in part 2 of the Appendix. 
It is shown there that for estimates of carried load, 


z= 0 
Z 
> a i oe 
r>1 Vz = 100 W/ wit (r — 1) (10) 
r<1 minimum V,=0 (r = 1, },3,h, ete) 
‘ r ; 11 
maximum V, = 100 5 ont (r = 2, 2, 2, ete.) (11) 


Since constant holding times found in practice are often very short, the 
case of r > 1 is the most likely to be met. For all values of r greater than 
one, the error given by formula (1) for exponential holding times is some- 
what greater than the error given by formula (10) for constant holding 
times, so use of formula (1) for the constant holding time case is con- 
servative. For values of r less than 1, the error is an oscillating function 
of r. The coefficient of variation varies from zero to 23 per cent above 
that for exponential holding times. Where 7 may not be accurately known 
the formula for exponential holding times again seems appropriate. 

In making estimates of the source load when the holding time is con- 
stant, if r > 1, each scan is uncorrelated with any other, since no call 
can be counted twice, and may be considered a random sample of traffic. 
There area total of Nc scans which have anaverage scan of a andstandard 
deviation »/a. The average error in estimating a is, therefore: 


5=0 
with coefficient of variation 


= / ane = / ar 
V. = 100 yg ~ 100 net (12) 
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Equation (12) may also be derived with the procedure used for equation 
1 
2 
= 1 equation (12) is approached by equation (9). For smaller values 
of r (but with r still greater than 1), Vs for constant holding times is 
- less than V, for exponential holding times. When r = 1, there is no 
carried load error. For values of r less than 1, the coefficient of varia- 
tion of error in estimating source load average will vary from 


t t r? 
/ an = mat +5) 


depending on the exact value of r. It is interesting to note that V; for r 


(9) using o; = o; + o,. For values of r large enough to make etn 


= 0.5 is the same as for r = 1.125. 


b. Loss 


The effect of loss in the group depends upon the disposition of the lost 
calls. In general, accuracy in measuring carried load increases with in- 
creased loss because under these circumstances fewer load changes occur 
between scans. This is evident in the extreme case of a group which is 
100 per cent loaded; a single switch count gives a correct reading for any 
length period. Obviously load readings at 100 per cent occupancy are 
not very useful in estimating offered loads since the amount of lost load 
cannot even be guessed at. However, in the cases of lost calls held 
(Poisson) or cleared (Erlang B), the offered load may be estimated from 
the carried load (less and less accurately as occupancy increases) and in 
the case of lost calls delayed the offered and carried loads are likely to be 
the same even at high occupancies. With high loss, therefore, estimates 
of source load are subject to errors not considered in deriving equation 
(99); however, switch count error in estimating carried load will be 
materially less than predicted by equation (1). 


c. Random Call Origination 


On trunk groups which are alternate routes, calls may no longer be 
considered as originating at random. The resultant grouping of call orig- 
inations wil] tend to decrease the accuracy of switch count measurements 
in estimating carried load; however, there is a corresponding decrease in 
accuracy in estimating the source load from the carried load so that ac- 
curacy in estimating carried load may be less worthwhile. 
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d. Statistical Equilibrium 


Statistical equilibrium may be thought of as the absence of trends in 
subscriber calling rates or holding times with the passage of time. The 
effect of trends on switch count accuracy in measuring carried load is 
very small except where the changes in traffic level are frequent and 
abrupt with respect to the scan frequency. Such traffic behavior is rare. — 

Trends within the busy hour complicate the problem of estimating the 
average source load. However, it can be shown that if the trends aresmall 
(in the order of 10 per cent to 20 per cent) little error is introduced by 
assuming that no trend exists. Large trends (in the order of 100 per 
cent), however, may indicate that the traffic source is so unstable that 
more hours of traffic data should be taken in order to insure that the 
sample is representative. 

Trends from day to day do not affect the source load estimates in the 
same way as within hour trends. The source loads are seldom exactly the 
same on any two days although in most offices a load pattern is repeated 
from week to week. The traffic engineer may be interested in the average 
source load of either a typical week day in the busy season or, some- 
times, of the average of the two highest days in the week. As long as the 
source load of each particular day remains close to the average for that 
day of the week, the general average for several different days of the 
week, will be known with about the same accuracy as if they had all 
come from a common source. If, however, there is no stable pattern in 
the source load, a third error in estimating the average is generated. 
There is some difficulty in determining whether or not variations in load, 
as indicated by measurements, are due to sampling variations or to an 
unstable source. Quality control methods might be used to detect in- 
stability but gathering and processing sufficient data for such an analysis 
might prove uneconomical. In general, if a traffic engineer feels that 
his source load is unstable he will need more hours of data than indicated 
by formula (9) to meet a given criterion of reliability. 


CONCLUSIONS 


A theoretical approach to the problem of the accuracy of switch count 
measurements in estimating carried load and average source load has 
been explored. It is believed that the assumptions made are satisfied 
sufficiently often in practice to enable fairly wide application of the re- 
sults of this exploration to traffic measurements. However, it should be 
kept in mind that where the assumptions are clearly not valid, special 
allowances will need to be made. In any case, the confidence placed in 
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usage measurements by a traffic engineer is a function of his experience 
and judgment. It is hoped that the results of this study will add to the 
knowledge essential to sound traffic engineering. 


APPENDIX 


DERIVATION OF SWITCH COUNT ERROR IN ESTIMATING CARRIED LOAD— 
WITH EXPONENTIAL HOLDING TIMES 


This derivation is based on a similar derivation by R. I. Wilkinson!. 
However, since load rather than holding time is of interest here, the em- 
phasis has been somewhat shifted. 

Assume that switch count measurements are being taken on traffic 
with: 

a. Calls originated individually and collectively at random 

b. Exponentially distributed holding times 

c. Negligible loss 


Let 7 = interval between scans 


t = average holding time 


, 


a’ = traffic carried, in erlangs 


T = length of observation period 


sel ie number of holding times in a scan interval 
T : 
c = = = number seans in observation period 
; 
T Oat : ; 
ae ae number of holding times in observation period 


N = number of observation periods. 


Consider that the observation period begins with the first sean and 
ends 7 time units after the last scan. It is desired to find the error in esti- 
mating the true load carried by averaging the number of circuits found 
busy on each scan. Following Wilkinson’s method we will first estimate 
the error of the switch count method in measuring the contribution of a 
single call to total usage and then modify it to take account of n calls. 
Calls of two types must be considered, those originating outside the in- 
terval and extending into it, Type I, and those originating within the 
interval, Type II. Both types may be subdivided depending on whether 
or not they extend beyond the end of the observation period. These are 
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indicated in Fig. 7. Only that part of a call which falls within the obser- 
vation period contributes to the usage of that period. First the error 
made by switch counts in measuring this contribution will be derived. 


Type I 


Consider a call which is already in progress at the start of the obser- 
vation period. Its duration beyond that point, according to theory, will 
be exponentially distributed about an average of 1. 

If this duration, t, is between 0 and 7, the call will be counted once (a 
measured contribution of 7 erlang hours) and a positive error of s = 7 — t 
will be made. The same error will be made if t = 22 — x so that the call 
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Fig. 7—Graphical indication of the two types of calls with their two sub- 
divisions. 


is counted twice and so forth. Summing all the ways of making an error 
x, we have: 


P(x) dz = fi — a) + f(2i— 2) + +++ f(a — x) (1) 
where f(t — 2x) is the probability of t = 7 — a and 
f®) = ro ae 
Calls lasting beyond cz neither start nor end in the observation period 
so that their contribution is measured without error. For these: 
P(0) = Pt > ct) =e” (2) 
Therefore: 


1 4-2 
Pepolt) da = 5 i dx 








€ i dz (3) 
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Letting 
e—1 , 
= a : as = = _— '. = ree = 5 
yr; € b=1-0'; K 2° ee 
Pzso(v) dx = e’e "K dy (4) 


The moment generating function My,(e) of y is: 


M(x) = P(y) dye” + ¢™ 


te « eS (5) 
e" —e 
eee ae RE oe 
Neglecting terms of order higher than a?, 
2 
Mia) = 1+ a(rK — 0’) + 5 (Kr? + 2b’ — 2rK) (6) 


Type II 
Calls of Type II may have either positive or negative errors given by: 


t+ 2 


Preo(x) dx = ; [fo(—a) + file — x) + fo(2i — 2) 





+o + fea((o — i — 2] 7) 
+ gol—2%) + ult — x) + go(2i — x) +--+ + gulle — It — 2] 


Pig) ges" hese) tO as eGo) 


1 





where f,(nt — x) = probability that a Type II call has length nz — x and 
ends before the end of the observation period. 


ni-2z 


c¢c— 1 
c 


dx 





ste t 
t 


(,(nt — %) = probability that a Type II call starts nt — 2x before 
the end of the observation period and ends after 
the end of the observation period. 


| 


1 ni-—z 


pe t dx 
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Equation (7) becomes: 

















1 ni—z = = 
Peodt) x = [PEED TS Ste Eee fae 
1 ni—2z (8) 
Parole) de = | ge ol ae 
Letting 
; =u) K= > e “" as before 
and noting that 
i —nr K — cb 
2, = a l —e” 
raternel(e 8) (25822) + S]a 
r c l—et TC 6) 





Pz>o(x) dx = € "e" (: — ¥) (x - as = *) | dy 
r el =e 


The moment generating function of this pair of equations is the sum 
of their separate m.g.f.’s: 


0 r 
Mula) = | Py<olye’ dy + Pysoly) dy e*” dy 


_g¢- Ker, o- Ky = Ker, (10) 
ro+tK +c Dares Wie ee le =e 
+ a(k + re — e Ke“) 
7 re(1 + a)” 
Neglecting terms of order higher than a”, 
Mr(a) 
(11) 


= fret a = 04) +2 [ (romn(2) -2) Gere +20] 


Now the number of Type I calls present in an observation is a vari- 
able—with average ‘a’’ and a Poisson distribution. Similarly the num- 
ber of Type II calls is a variable, independent of the number of Type I 


: fi ; ae ae 
calls, with an average of “‘a oe or “arc” and a Poisson distribution. Ac- 
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cording to the laws governing the compounding of variables the moment 
generating function of the sum of n variables y, when v is also variable 
with generating function G(d), is G(/(e)) where M(a) is the moment 
generating function of y. 
The generating function of a Poisson variable with average “a” is 
6 7*' so that 
G(My(a)) = @ heer) 


G(Mn(e)) = g eter Mi (e) 


These independent variables may be added by multiplying their mo- 
ment generating functions to give the m.g.f. of the total measurement 
error of the carried load 

M(a) = ge atare)-taMy(a)-tareMz7() (13) 


(12) 


From (6), (11) and (18) the following parameters are found: 
g=0 (14) 


K , 20’ K K , 20’ 
oh = aro[(rten( 5) - 2) (1-7 +5) + (2-20 + 5) 


If, now, rc is sufficiently large 


oy = arc  retan(5) — 2| (15) 


It is more convenient to deal with the standard deviation expressed 
as per cent of the carried erlang load, the coefficient of variation. This 
is done by multiplying both sides of equation (15) by # to convert the 
time dimension from holding times to hours, dividing by 7’ to convert 
from erlang-hours to erlangs, dividing by a’ to convert to proportion of 
carried load, and multiplying by 100 to convert to per cent. Assuming 


a, ‘ 
qs approximately a: 


Vz = 100 





When JN observations are made this reduces further to 


“ am Ae 
Vz = 100 NT | retnn(5) | 


3 5 
Now etnh(a) = <+ 5-4 3 — (x? < 4) 
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and forr < 2 
Therefore 
retnh (5) = 2 i 
Ve = 100 wane = / on (r <2) (16) 


The error in carried load may be considered as the sum of a large num- 
ber of independent errors. Its distribution may, therefore, be expected 
to approach the normal distribution. Comparison of the third and fourth 
moments of the normal distribution with those of the error distribution 
(which may be obtained from equation (13)) show good agreement for 
values of a’ greater than 1. 


DERIVATION OF SWITCH COUNT ERROR IN ESTIMATING CARRIED LOAD WITH 
CONSTANT HOLDING TIMES 


Wilkinson has shown! that, for constant holding time, switch count 
error in measuring the holding time of one call has an average 


and standard deviation 





where JT >t 
T>1 
2, = negative error 
X2 = positive error 


Divide the problem into two parts: 


1. Forr > 1 
t= —t 
t=i-—t 





n= Vii -P =ivr—1 (17) 
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2. Forr <1 
Min. oz = 0 for r = 1, 3, 4, 4, ete. 
t t 1 -f 
. =| -—-——  —_—_—_-—— ll Cee il -——- Paes 18 
Max. oz (5 55 bs (18) 
for r = 4, 2, 2, etc 


Expressing this error in terms of carried load and proceeding as in 
Part I of the Appendix 


t 
oe ee | Ve = 100 4/ oir — 0 (19) 
2.r< 1 Min. V, = 0 
_ r t 
ell (/ ees 
c/T 4a'NTt 


Equation (20) compares favorably with the exponential holding time 
coefficient of variation of error of 


100 / i 
cJ/T VY ba'NTI 
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Network Representation of Transcendental 
Impedance Functions 


BY M. K. ZINN 
(Manuscript received November 5, 1951) 


The purpose of the paper ts to show that the admittance or impedance of 
certain continuous structures, such as, for example, a finite length of trans- 
mission line of any sort, or resonant cavity, can be represented exactly at 
all frequencies by a network comprising lumps of constant resistance R, 
inductance L, conductance G and capacitance C. The network will contain 
an infinite number of branches, in general, although a finite number may 
be used if it is destred to represent only certain modes. 

The procedure ts based upon a proposition known to students of function 
theory as “Muttag-Leffler’s theorem,” which amounts, roughly, to an ex- 
tension of rational functions to apply to transcendental functions of the 
type encountered in the theory of continuous structures. 

Several illustrative examples of the network synthesis are given. 


GENERAL 


Students of network theory are familiar with the fact that the im- 
pedance at a pair of terminals in a linear network comprising a finite 
number of resistors, inductors and capacitors, connected in any manner, 
is a rational function of the frequency having, in general, the fractional 
form of one polynomial divided by another. They are also familiar with 
the partial fraction rule whereby the function can be broken up into a 
series of elementary fractions, each of which exhibits one of the poles of 
the original function. This form is sometimes useful in the problem of 
network synthesis, where the impedance function is given and the ob- 
ject is to find a network having this impedance. 

The purpose of the present paper is to show how a similar procedure 
can be carried out for certain transcendental impedance functions per- 
taining to structures having distributed constants, such as, for example, 
a resonant cavity or a piece of transmission line. The method employs a 
well-known proposition of function theory, which is usually referred to 
as Mittag-Leffler’s theorem. This theorem provides a tool for breaking 
up a transcendental meromorphic function into an infinite series of 
simple fractions in much the same way as the partial fraction rule is used 
to break up a rational meromorphic function. The series representation 
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provides a means of determining a network of resistors, inductors and 
capacitors that will have an impedance equal to the specified transcen- 
dental impedance function. This process will be referred to as obtaining 
a “network representation” of the function. If the given function is the 
impedance of some continuous (i.e., non-lumped) electric structure, the 
result will be an equivalent network for the structure. For other pur- 
poses, such as, possibly, analogue methods of computing, the given func- 
tion may not arise from any electrical structure. In either case, the net- 
work representations to be derived are possible only if the function 
satisfies certain restrictions, which are stated in the section immediately 
following. 

The discussion is confined to transcendental impedance functions be- 
cause of the technological interest in the electromagnetic structures with 
which they are associated and because they have not received as much 
attention as rational functions in the literature dealing with network 
synthesis. The problem with which this paper is concerned can then be 
stated as follows: given, a transcendental impedance function satisfying 
certain conditions: to determine a network comprising elements of 
constant resistance, inductance and capacitance whose driving-point im- 
pedance function, at a pair of terminals, will equal the given function at 
all frequencies, real and complex (except at the poles). 

For illustration of the procedure, three examples are given. The first 
is the impedance of a short-circuited or open-circuited transmission line 
in which the distributed primary constants, R, L, G and C are assumed 
to be invariable with frequency. The second and third examples are the 
impedances of resonant cavities driven in two different modes. In these 
examples the variation of resistance with frequency, due to ‘“‘skin-effect,”’ 
is taken into account. 


IMPEDANCE FUNCTIONS 


The functions under discussion will be referred to as “impedance 
functions” with the understanding that the term is meant to include 
“admittance functions” as wel]. By reason of the duality principle that 
runs through all electric circuit theory, any general proposition devel- 
oped for one must apply to the other. The functional designation, F'(p), 
will be used to denote either an impedance or an admittance function. 
When a distinction is necessary, the impedance will be designated by 
Z(p) and the admittance by Y(p). The independent complex variable p 
is the generalized radian frequency. (For sustained sinusoidal currents 
and voltages, p = tw = 2zif where f is the real frequency.) 

For the applications contemplated, F(p) is a transcendental mero- 
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morphic function, which term implies that the function is given by the 
ratio of two entire functions, one or both of which is transcendental, and 
that the singularities of the function are ordinary poles, except for the 
point at infinity, which is an essentially singular point. In order to realize 
the particular network developments to be given, it will be supposed 
that the function satisfies the further restrictions given below: 

(1) All the poles lie in the left half of the p-plane with none on the 
imaginary axis. 

(2) F(p) = F(p). (The superbar denotes the complex conjugate of 
the unbarred symbol.) 

(3) Real part [F(7w)] > 0 for all real values of w. 

These three conditions are necessary to insure that the function is the 
impedance of a possible linear, passive electric circuit structure. Inter- 
preted physically in terms of this possible equivalent structure, the first 
condition specifies that the structure shall be stable; that is, every natu- 
ral mode of oscillation dies away exponentially. The second condition 
specifies that the natural oscillations are real functions of time. The 
third condition specifies that if a sinusoidal current flows at the driving- 
point terminals of the equivalent structure, the average real power de- 
livered to it will be positive. Since these three conditions, or their equiva- 
lents, are frequently mentioned in discussions of network theory, it is 
assumed that they are understood without more detailed explanation. 

In addition to the above restrictions on the form of the impedance 
function, the following two conditions, while not necessary, will be im- 
posed to limit the scope of the discussion: 

(4) All the poles of: F(p) are simple. 

(5) F(p) = 0(1), exactly, as |p| — © everywhere except at the poles. 

Condition (4), while limiting the scope of the exposition required, does 
not restrict the application of the results in any important way, because 
most impedance functions for which a network representation may be 
required have only simple poles. 

Condition (5) implies that as p increases along any straight line drawn 
through the origin and not passing through any pole of F(p), the modu- 
lus of F(p) either approaches a limit or oscillates between finite limits. 
The physical implication of this condition is that the response of the 
network as a function of time to a suddenly applied cause begins with a 
discontinuity of the same degree as that of the cause. For example, the 
current response of the network to an applied step of voltage begins with 
a finite discontinuity. This behavior is a characteristic of continuous 
(non-lumped) electromagnetic structures, which furnish the principal 
application of the network developments to be described. 
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MITTAG-LEFFLER’S THEOREM? 
Let the poles of the given function F(p) be pi, pe, pz +++ , Where 
0<|m| <|pe| <|psl--- 


and let the residues at the poles be A1, A2, A3 - +: , respectively. Suppose 
that it is possible to draw a sequence of closed contours, C,, such that 
C’, encloses pi, pz, °** Dn but no other poles and such that the minimum 
distance of C, from the origin tends to infinity with n. Suppose also that 
F(p) satisfies conditions (2), (4) and (5) above. Then Mittag-Leffler’s 
theorem gives the following series development for F'(p): 


eee A A 
F(p) = F(0) + Limit ny An 1 
(p) = FO) + Limi p24 (. oa As) (1) 


The notation here used employs the convention that 





pn = Pn and An = An, 


since, by virtue of condition (2), the poles occur in conjugate complex 
pairs. The value, n = 0, then allows for a pole on the negative real 
axis. 

Given any suitable function, the procedure is to determine its value 
for p = 0 and the location of its poles. The residues are next determined 
by 


A, = Limit (p — p,)F(p). 
P—Dn 
Then the Mittag-Leffler expansion can be written down at once. 


NETWORK REPRESENTATION 


In the:series (1) the terms occur in pairs with conjugate complex poles 
and residues. The object is to obtain a network representation of each 
such pair of terms. If F(p) is taken as an admittance, the branches rep- 
resenting the pairs of terms will all be connected in parallel; if F(p) is 
taken as an impedance, they will all be connected in series. 

Methods for obtaining a network representation for a rational func- 
tion, such as the one comprising a pair of terms in the series (1), are well 
known. It is only necessary to describe certain procedures of particular 
application to the present problem. Brune’ has stated that the necessary 
and sufficient condition for a network representation of a rational func- 
tion of p to be realizable is that it be a ‘“‘positive real function,” that is, 
a function that is real for real values of p and whose real part is positive, 
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or zero, when the real part of : is positive, or zero. In view of conditions 
(1) and (2) above, only one test” need be applied to each pair of terms 
of the series (1): the sum of a pair of terms will be a positive real func- 
tion if, and only if, the real part of their sum is greater than, or equal to, 
zero for all purely imaginary values of p. 

The general term pair for which a network representation is sought is 
P— Pr P— Pn Pn Pn 
Evidently two cases can be distinguished at the outset, depending upon 
whether P,,(0) is positive or negative. If P,,(0) is positive, the network 
branch, in order to be realizable, should be designed to represent P,(p). 
The left-over negative term, —P,(0), then can be absorbed in the posi- 
tive first term, /’'(0), of the series (1); more will be said of this later. If, 
on the other hand, P,,(0) is negative, the network branch should repre- 
sent the whole term, P,(p) — P,(0). This procedure insures that the 
real part of the branch impedance will be positive, or zero, at zero and 
infinite frequencies. To guarantee that the resistance is positive at all 
other frequencies requires further tests now to be specified. 

Let the real and imaginary coefficients of the poles and residues of 
the n™ term be 








F.(p) = = P,(p) — P.(0) (2) 


Pn = An ae Bn Pn = —Qn — 1B n 
An = Gn + ibn, An = Qn — in 


(With this notation, a, and 8, are always positive; a, and b, can be either 
positive or negative.) Then (dropping the subscripts) 


2(aa — bB) + 2ap 


PO) = at BE tap +p 
. 7 _ 2(aa — bB)(a” + 6°) + 2w'(aa + bf) 
BON EF BEE Dhe — B) Fal ia 


Olan = 3a°bB — 8aa08" + bB°)w" — 2(aa — bB)w* 
(a? =P B*) Ca? = 87)? 2 (2 Ba? > a] 
The necessary and sufficient conditions” for the real part of a rational 
function of p to be positive, or zero, for purely imaginary values of p are 


that the function be positive for p — +7 and have no imaginary roots 
of odd multiplicity. When this test is applied to the functions P(p) and 


R[P(iw) — P(0)] = 
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P(p) — P(0), as given by (8), the following conditions are obtained: 
P(p) will be a positive real function if, and only 1, 


aa — bB > 0; Le. P(0) > 0 (4) 
and 

aa + b8 > 0 
P(p) — P(O) will be a positive real function tf, and only if, 

aa — bB <0; i.e. P(O) <0 (5) 
and 


ao® — 30°bB — 3aap’ + bp* < 0. 


If all terms of the series satisfy one or the other of these conditions, 
network branches can be devised to represent all the terms and all the 
R, L, G, C elements of the branches will be positive. 

In case all the terms are of the type where P,,(0) is positive, so that 
the network branches are made to represent P,,(p), the left-over constant 
terms can be collected and added to the first term, F(0), of the series. 
This collection of terms then must be represented by a final branch of 
pure resistance, or conductance, of value, 


FO) — >0 P,(0) 
n=0 

If the sum of the variable terms approaches zero for p — -i~, the 
final constant term supplies the high frequency resistance of the func- 
tion F(p) and since this must be positive, if condition (3) is satisfied, the 
final resistive element will be positive. If the series converges non-uni- 
formly, the sum of the variable terms can have a value other than zero 
as p — -kio~ in spite of the fact that every term approaches zero indi- 
vidually. In that case (see example 1) all or part of the high frequency 
resistance may be supplied by the sum of the variable terms. 

In case all the terms are of the type where P,,(0) is negative, so that 
the network branches are made to represent the sum, P,(p) — P,(0), 
of the variable and constant terms and the series is uniformly conver- 
gent, all the high frequency resistance is provided by the branches rep- 
resenting these terms. The first term, F(0) then supplies the de re- 
sistance, which is positive by condition (8). Non-uniform convergence 
can modify this division of high- and low-frequency resistance, however. 

Cases can arise in which the series contains terms of both types. In 
such a case the de resistance, or high frequency resistance, or both, of 
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the given function might be less than the sum of the variable terms for 
these frequencies, with the result that the final resistance branch would 
be negative for either the series or parallel type of network development. 

To make the procedure as concrete as possible, particular forms of 
networks are described in the section following with explicit formulas for 
computing their elements. 


NETWORK FORMULAS 


Simple forms of network branches are shown in Figs. 1 and 2. Those 
of Fig. 1, referred to as branches of “‘the first kind” are suitable for con- 
nection in parallel where the given function F(p) is an admittance, Y(p), 
while networks of “the second kind,” shown in Fig. 2, are suitable for 
connection in series to represent an impedance, F’(p) = Z(p). The net- 
works of Figs. la and 2a apply where the value P,,(0) of the general term 
is positive, while Figs. 1b and 2b apply where P,,(0) is negative. Figs. 3 
and 4 illustrate, respectively, networks of the types of Figs. la and 2a 














Cn Ln 
n n n Cn 
1 
or Gn 
(a) (b) 
Fig. 1—General branches of the first kind. 
Fig. 1a Fig. 1b 
(use where F'(p) = Y(p) (use where F(p) = Y(p) 
and Y,(0) > 0) and Y,(0) < 0) 
wo _ Balan + Bn) Can + Da) 
Dg = 5a, Ly = WE 
ee ei Mee M? 
Let (* v 1) LC, ~ Ba, + 84) 
Gy 1 _ Ann — OnBn (6) 
C, = a, (QnQn DrBn) Galan = VF 
Rn _ 1 N 
G, = Y(0) — >» Y,(0) G, = Y(0) 
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connected to form the completed network with the final non-reactive 
branch, G, or R, , in place. 

Formulas for the network elements are obtained by equating the poles 
and residues of the network impedance function to the given poles and 
residues of the general term of the series. Since both poles and residues 
occur in conjugate complex pairs, and since equality of real and imagi- 
nary parts is involved, there are four equations, which are necessary and 
sufficient to determine the four constants, R, L, G, C, of the network. 
The formulas that are obtained by solving these equations are given 
beneath Figs. 1 and 2. 

The values given for G, and R, in each case assume that all the terms 
of the series are of the type specified for that case. 








Cn Ln 
Gn 
000 
Re be Rn Cn 
(a) (b) 
Fig. 2—General branches of the second kind. 
Fig. 2a Fig. 2b 
(use where F'(p) = Z(p) (use where F(p) = Z(p) 
and Z,(0) > 0) oa Z;; a : _ 
aes ae Palen of Seen te) 
1 be 1) Te ies MM 
Lea LC, ~ Ba, +B) 
Rr = 1 = = __ AnOn = baBn 
E. _ a, (Anon bnBn) BaCx _ oat | ne (7) 
Gn 1 N 
Ga EN ae Ee 
R, = Z(0) — >> Z,(0) R, = Z(0) 


n=0 


where M = a,(8;, — of) + 2etnBnda 
N = —ane®, + 307,biBn + 3dnanB, — dnB* 
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In the case of the parallel-type networks (Figs. la and 1b), pr = 
—an + iB, is a pole of the ddmittance, Y(p), and A, = dy + 7b, is the 
corresponding residue. In the case of the series-type network, the same 
symbols represent a pole and residue of the impedance, Z(p). 

The networks specified by Figs. 2a and 2b are duals of the networks 
of Figs. 1a and 1b, respectively, and are obtained from the latter merely 
by replacing LZ, by C,, Ra by Gn, and vice versa. 

The formulas are intended to apply to complex poles. They can be 
applied to real poles by taking b, and 8, equal to zero and doubling the 
residue, @, , but this procedure is unnecessary, because the network rep- 





Fig. 3—Network of the first kind (branches 1a). 





Fig. 4—Network of the second kind (branches 2a). 


resentation of a real pole can be found readily enough by inspection of 
the impedance terms involved. (See Example 1.) 

The above discussion is intended to sketch a general picture of the 
procedure. Individual cases may involve considerable detail that can be 
understood more readily by reference to the next section. 


APPLICATIONS 


Example 1a: A transmission line with its far terminals short-circuited 
affords a simple illustration of the equivalent network theory. Let it be 
assumed that the parameters, R, L, G and C of the line are constants. In 
the more advanced examples to follow, the variation of these parame- 
ters with frequency for a particular kind of line will be taken into con- 
sideration. 
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The impedance of the short-circuited line (Fig. 5) is 
Z = Z, tanhT (1-0) 


where Z, is the characteristic impedance and I is the total propagation 
constant of the line. We have 


1/2 
r= ([(R+ pL\(G+ pcyy? (1-2) 


Rk, L, G and C being given for the total length of line. 
To obtain a development in terms of network branches of the kind - 
shown in Fig. 1, we consider the admittance function, 


Y = Yo cothT (1-3) 


where Y = 1/Z and Y» = 1/Z,. Our first task is to find the poles of this 
function and the residues. Since the complex frequency variable p occurs 


R,L,G,C 
Q 
=f 
0) 


Fig. 5—Short-circuited transmission line. 


under square roots in both Zp and I, it might be suspected, offhand, that 
the singularities of the function are branch points rather than poles. 
Such is not the case, however. There are no branch points and all the 


poles are simple. 
The singularities of Y are to be found among the zeros of tanh I, 


which occur at 
I = trn, n=0,+1,4+2,+3,-:: (1-4) 
To determine them, we solve 
I’ = (R+ pL)(G@ + pC) = —2'n’ (1-5) 
and find these roots: 
Dn = —On + Bn, Dp-n = Pn = —On — Bn 


where 


G R 
= 50 oy, 


2°22 2771/2 
Bn = ae = (S = ail (n>0) (1-6) 
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For n = 0, the above would give 


ee em 
Po i Ps C 


But if we let T — 0, so that tanh T — I, we find that only the point, 
—R/L, is a singularity of Y; the other point, —G/C, is a regular point. 
Therefore Y has only one real singularity. 

To find the nature of the singularities of Y, we next calculate 


= Dn : 
Lat at ni Tone As oe 


and find that at each p, the limit exists and has the value 


; R 
Wan —-S 
1 1 ( a r) ; 
a —— (p,.)  L BrL at — 


where I’ (pn) = 5 I'(p), evaluated at p = p,. The fact that this limit 


exists shows that all the singularities are simple poles. The values of An 
are then the residues at these poles. 

When we now apply formulas (6) to determine the elements in the 
general branch of the equivalent network of Fig. la, we obtain, for 
n> 1, 

2.2 
L, = so lt = as G, = G Rn = ie (1-9) 
2 LiGa LC C, C os L 

The network then comprises an infinite number of such branches in 
parallel. Each branch has the same elements #,, and L, , equal, respec- 
tively, to half the total resistance and inductance of the transmission 
line, but the elements G, and C,, decrease from one branch to the next in 
inverse proportion to the squares of the integers. 

The Q of the n™ branch, which can be regarded as the Q of the asso- 
ciated resonance of the short-circuited line, is 





(1-10) 


il Go  /re @ 
oa pe - B= -& (1-11) 


where 
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Thus, for small dissipation, the resonances would become sharper in di- 
rect proportion to the frequency (if the parameters R, L, G, C, were 
invariable with frequency, as assumed). 

The above described branches of the equivalent network account only 
for the complex poles (n > 1) of the admittance function. Two more 
branches remain to be calculated. One is for the real pole (n = 0), which 


occurs at pp = —R/L, with residue, Ay = : The required branch for 
this pole is 

Ao _ t 

Pp — Po R + oL 

The other is the final conductance branch, which is calculated as follows: 


= Y(0) + 4 “" = * = 4/S coth VER - } 


n=—-% 





(1-12) 


(1-13) 
= 1 
PaO De aa GR 


so that, for this example, the conductance branch vanishes. The network 
is drawn in Fig. 6. . 
A series type of network, as shown in Fig. 7, can be determined by 





nIO 
nla 





2R 2L 2R 2 
mea)’ we($/° iad 9 aia 


Fig. 7—Network of the second kind equivalent to the short-circuited line of 
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similar means. Since, however, it is a dual of the parallel network of Fig. 
9 for the open-circuited line, next to be discussed, it can be drawn im- 
mediately, without further calculation, once the latter has been found. 

Example 1b: We now calculate a network for the same line with its 
far terminals open (Fig. 8). To obtain a network of the first kind, with 
branches in parallel, we deal with the admittance function, 


= Y)tanhT (1-14) 


The singularities of Y are found among the zeros of coth I, which occur 
at 


TI = in(n + 4), n=0,+1,+2,4+3,--- (1-15) 


The points p = —R/L and —G/C are both regular points this time. 
(—G/C is a zero of Y.) The singularities are simple poles, as before, with 
residues, 

1 


Tale) oe 


as before. 

The network branches for the complex poles are therefore obtained 
merely by putting n + 4 in place of the 7 in all formulas of the short- 
circuit network. There is no branch corresponding to the branch R + pL 
of the other network and the conductance branch is again found to be 
zero. The complete parallel network is drawn in Fig. 9 and the series net- 
work, in Fig. 10. 

It will be observed that the series network of Fig. 10 is the dual of the 


Riise 





—Z=y 


eo) 





Fig. 8—Open-circuited transmission line. 





Fig. 9—Network of the first kind equivalent to the open-circuited line of Fig. 8. 
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parallel network of Fig. 6 and the series network of Fig. 7 is the dual of 
the parallel network of Fig. 9. These dual relationships are of course a 
result of the fact that the impedance of an open-circuited line is the dual 
of the impedance of the same line when short-circuited. 

Example 2: Short-circuited Concentric Line (or Toroidal Cavity with 
E| Radial). The preceding example considered a fictitious transmission 
line of invariable parameters, R, L, G, C, having a perfect short circuit 
at one end. The present example has to do essentially with the same 
problem but considers it from a more practical point of view. The vari- 
ation of & and L with frequency is taken into account and the impedance 
of the ‘‘short-circuit”’ is no longer neglected. - 

Let the line be the piece of coaxial cable plugged at both ends with 
conducting material as illustrated in Fig. 11. Considered from an alter- 
native point of view, our line is now a toroidal cavity oscillating in the 


Cc Cc 
2 2 





ee Fig. 10—Network of the second kind equivalent to the open-circuited line of 
ig. 8. 


mode where the electric force / is directed radially and the magnetic 
force H lies in planes at right angles to the axis. If we assume the cavity 
to be excited, or “driven,” from one end,* the impedance that is effective 
in defining the selective characteristic of the cavity with respect to fre- 
quency is the total impedance at that end, that is, the sum of the im- 
pedance Z;, viewed into the cavity, and the impedance, Z:2, of the ad- 
jacent end-plug. Therefore, we have to deal with the impedance, 


Z=F,+ Ze. (2-1) 


By “impedance” is here meant the same thing that one considers in look- 
ing at the problem from the point of view of transmission line theory, 
namely, the complex ratio, for exponential oscillations, of the voltage 
between the inside and outside cylindrical surfaces to the total current 

* For determining the ‘‘natural frequencies”’ of oscillation of the cavity, it is 
immaterial at what point along it the impedance is taken; the total impedance 
at every point has the same roots. The impedance is, nevertheless, not the same 


at all points so that the behavior of the cavity, when driven, will depend to some 
extent on the driving point. 
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flowing axially in the inner conductor at the same point. The zeros of Z 
define the natural frequencies of oscillation of the cavity and their’ asso- 
ciated damping constants, or Q’s. Our task is to develop an equivalent 
network for this Z. 

We have 


_ _7 (itoe™ , 1+p 
Zim ty + y= ty (tte 4 ite (2-2) 





2hc 3, 3, 
7 2n2 (- 59% * Th; ) 





AL 1 1 

= tSag- *onae) 
baal oe ae 3_) 
Won C 2d, 2hdn 











Ho b ma rs 8 
L=S2loc= ve = 3 (109) 
27 a VUE 
— 27€0_ G,= Won HY 
= 5 = 
LOGS 


-9 ; 
(€.9.) Mo= 47 (10"7), €9= 35, FOR AIR IN M.K.S. UNITS 
M=Uo, J=5.8(107), FOR COPPER IN M.K.S. UNITS 
(N= 1 FOR FUNDAMENTAL MODE) 
Fig. 11—Toroidal cavity, # radial. 
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where 


1/2 
Ly = ( : r) , y= (R+ ply (G+ pcy”? (23) 


G+ pC 
ZL, — Zo 
= 2-4 
BS ear 7. (2-4) 
i Mig? : 
Zo = 5 log = (2-5) 


— 9 In(ca) 7 Ko(ob) Pio b 
TPE Gee a) 2b GG) ee ea. 





Z 
G+ pC = aoe (go + peo) (2-7) 
log : 


h, a, b = cavity length, inner radius, outer radius, as shown in Fig. 11, 


all measured in meters 
ae (2)" 
g (2-8) 


o = (pug) ” 


H, g are permeability, conductivity of the conducting material of the walls 
(for copper: » = 47(10-’), g = 5.8(10’) in M.K.S. units). 

Ho, Jo, €o9 are permeability, conductivity, dielectric constant of the dielec- 
tric material occupying the cavity (for air: wo = 47(10~), go = 0, 
€) = (10° °)/36z in M.K.S. units), p = generalized frequency vari- 
able. 

T(z), [1(z) are Bessel functions of the first kind for imaginary argument 
and of order 0, 1. 

K,(z), Ki(z) are Bessel functions of the second kind for imaginary argu- 
ment and of order 0, 1. 

Except for ignored small deviations of the field around the corners of 
the cavity, the above formulas are exact. To arrive at results that are 
sufficiently compact to be useful, we make these approximations, at the 
start: 

1/2 
Zy = Ky = H = 7 log (2-9) 
where 
1/2 
no= || = 1207 ohms (2-10) 


€9 
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From this, 





1 — No 
= 2-11 
Oe eee (2-11) 


Having in mind microwave applications, where the moduli of the argu- 
ments of the Bessel functions are >3000, we take 

Iy(z) = Ko(z) _ 

I,(2) Kyx(2) 








so that 
n 1 1 PKHo b 
R+opb = —(-+ 5 -- Gq SG (2-12) 


2r \a T 
Also, we have in mind only air dielectric and assume any loss therein to 
be negligible; that is, we assume G = 0. 
All further approximations that are made are either 


1 


a 1/2 
pS Le oor (1+ 2A)" =1+A 


where, for an air-space enclosed by copper walls, and for frequencies on 
the order of 30,000 megacycles, A is on the order of 10 *. For cavities 
made of other materials, the results obtained may not be sufficiently ac- 
curate and the problem would have to be reviewed from the start. In 
particular, the results do not hold for a cavity having walls of magnetic 
material, because we assume here that the permeability of the metal walls 
is the same as that of air;ie., u = uo. 

To obtain an equivalent network of the first kind, we deal with the ad- 
mittance, which is, from (2-2), 


_1_ iy (— (1 — pe) 
Y oa % = Ho 2(1 = pre27") (2-13) 





where Hy = 1/Ko. 

The poles of Y are then the zeros of 1 — pe”, which are obtained by 
successive approximations. We first make a close estimate of the zeros by 
assuming that the impedance of the short-circuiting plugs is zero; that 


—2yh 


is, we assume, Z_. = 0, whence p = —1. To obtain this estimate, we have 
to solve 
h 9 1/2 
yh = h(i 42) = rin (n = +1, 42,43---) (2-14) 


where 
_ 2ab log (b/a) 
a+b 
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and v = 3(10°) meters per second. The approximate solution is 


1 
n —=Po0n 1 
Pl Po ( a ig x) 





where 
imnv 
Pon = _ and Con = (ponug)” 
Next we improve our estimate of the zeros by the well-known method in- 


volving the derivative of the function, 1 — pe °”, with respect to p, 


evaluated at pi,. This now takes account of the actual impedance of the 
end-plugs. The values of the zeros, so obtained, are 

1 a i 2 ie wa |S) Don = Pn = — Qn — 182 
where 





] 1 
Qn = Won & si =) 


1 1 


where 6,* is the real part of oo, . That is, 
bn = (wonug/2)? 


(2-15) 





where 


Tv 


On = >. 


h 


As an incidental matter of interest, the above gives the Q of the cavity 
at any resonance, namely 








1 
Qn = aes (2-16) 
2an (eet 2d 
h 
For example, the dimensions, a = .5 cm., b = 1.0 em., h = .5 em. pro- 


vide a cavity that resonates at about 30,000 megacycles. Then the Q’s at 
the first three resonances would be as follows: 


W0n 
1 30,000 x 10° 4250 
2 60,000 x 10° 6010 
3 90,000 x 10° 7360 


* For any frequency, 6 = (wug/2)/? issometimes referred to as the ‘“‘skin depth’’ 
because it is the depth of metal at which the current density falls to 1/e times its 
value at the surface of the metal. 
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The importance of including the effect of the end-plugs in determining 
Q is shown by the fact that, if they were assumed to have zero impedance, 
Q at the first resonance would be 12,120 instead of 4250. 

To determine the residues at the poles, we write 


Yeah (1 — p)(1 — pe") _ F(p) 





2(L = pret) Gp) oe 
and then the residue at a simple pole p, is 
Ae AMD) (2-18) 


G'(pn) 
This limit is found to exist, showing that the poles are, in fact, simple. 
The value found for the residue, A, , is 


Aptech Dig Ag = Ay = ae = th 








Qn, = Hoon 1— aS os a) 

"en Qdin 2hép (2-19) 
= Hywon ( 1 1 

ee mn \2d6, T an) 


When formulas (6) are applied to determine the elements of the tuned 
branches of the equivalent network of the first kind, the results are, for 


the n™ branch, 
~ Korn 1 ] 
ere (1 a * a) 


3 
| 








c 7 a =) 
ae Won 1 + =— Le. 
Lin a hin (2.20) 
Gn _ 0m 
Ci 2nd», 
Tun. 1 3 
T om (2 + ge) 


In terms of the R, Land C of the picce of coaxial line, the elements of the 
** branch are as follows: 


Ln = Mat oe +o) 


2dé, | 2h6, 
hk 3d 
By = (a + 


(2-21) 


| 


2he 3 3 
Cr ea - sx + ax) 


_ wonC = 3 3 
= ns, (1 ods, 7 ) 
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where 
_ Ho, 6 — 1 (oon)? (1 ‘) 
ba glog?, Ras (SH) (G4 
2rre nv ae 
C = n> a On — one 
lo b a h’ ( 2 
5G 


The network is shown in Fig. 11. 

It will be found that a “leakage” element, G, , appears in the equiva- 
lent network, although the air dielectric in the cavity was assumed to 
have no leakage (G = 0). This element arises from the end-plugs and is 
necessary to account for the dissipation in them. 

To obtain a network exactly equivalent to the cavity at all frequen- 
cies, we should add a branch corresponding to n = 0, as was done in 
example 1. This branch would make the equivalence hold down to and 
including zero frequency. But, inasmuch as the approximations that have 
been made hold only for the high frequencies, where the resonances oc- 
cur, it would be inconsistent to add this branch. What has been arrived 
at, then, is a partial network representation that gives a close approxi- 
mation to the impedance of the cavity at high frequencies, only. 

Example 3: Toroidal Cavity with E Axial. For further illustration, we 
consider another mode of oscillation of the short-circuited concentric 
transmission line investigated in the previous example. This time it is 
assumed that the radial electric force vanishes while the axial electric 
force between the end-plugs exists. The magnetic force is directed in 
circles concentric with the cylindrical central conductor, as before. This 
situation is illustrated in Fig. 12, which is the same as Fig. 11, except for 
the new disposition of the /-vector. 

For the new mode of oscillation, where the wave is a cylindrical one 
propagated back and forth between the inner and outer conducting 
cylinders, the oscillatory space is naturally thought of as a “toroidal 
cavity,’ while, in the previous example, where the wave was propagated 
axially back and forth between the terminal discs, the space was called 
a “concentric line.” Actually, the cavity itself has the same geometric 
form in the two cases. A practical distinction may exist, however, in that 
the axial mode of oscillation could be more easily excited in a cavity 
whose axial length is large compared to its radius, while the cylindrical 
mode would arise more easily in a flat “pillbox”? cavity whose radius is 
large compared to its axial dimension. 

The approach to the problem will be that of transmission line theory, 
as before. This time, the “line”? comprises two circular discs between 
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which the cylindrical wave is propagated. The series impedance and 
shunt admittance of such a line are functions of the radius and so will be 
designated Z(r) and Y(r), respectively. Their values are given below: 


Zy) = A ne (3-1) 


V(r) = Meare (3-2) 


These formulas take into account the losses in the flat walls but assume 
the conductance of the air between them to be negligible. Losses in the 
inner and outer “short-circuiting” cylinders will be taken into account by 
the boundary conditions. 







Lg HY) 











2(b-a)C ('- 34 3 ) 
7 ene 2hén 2 (b-a) dn, 





(b-a)L =a 
(1 + Shae ‘cs i) 








(b-a)R h 
> (1+ Spal) 





Won C (1+ Shd; Bb=aldp) 














=o 2g Won b-a 
. Moh 1 
V= = 3 (108) 
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277 ao 6, = \ Pond 
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SEE FIGURE 11 FOR fo, &,/59 
(N=1 FOR FUNDAMENTAL MODE) 


Fig. 12—Toroidal cavity, # axial. 
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If V is the voltage between the flat faces of the cavity at a radius r and 
I the total current in the lower face at this radius, we have - 


dV 





— = —I[Z(r) 
dr 
dl oa 
By differentiating, 
dV dZ dI__ (1 dZ\aVv 
ie ae oa oe) et PAY oa 
But 
ZY = (2n + iaph) 52 = 7’ 
which is a squared propagation constant, independent of r, and 
1 G8 S20 
Z dr r 
Therefore, 
OV 2 av Gb. 
a ea (8-5) 


is the differential equation for the voltage. The usual solution of this 
equation is a linear combination of Jo(yr) and Ko(yr) but since, in this 
case, the arguments will be almost purely imaginary, it is more con- 
venient to employ the pair of functions, Jo(—7zyr) and No(—7tyr). 

The solution for the voltage between the upper and lower surfaces at 
radius 7 is 


V(r) = AJo(—tr) + BNo(—-11) (3-6) 
and, from this, the total radial current in the lower surface, at that 
radius, is 
an dV(r) 
Z dr 





10) == = —1Y(r)[AJi(—tyr) + BNi(—tyr)] (3-7) 


where 
Yo(r) = 1/Zo(r) = (Y()/Z@)I” 
The impedance at the inner radius a, looking outward, is then 


V(a) AJo(—tya) + BNo(—tya) 


ta(a) = Fy AJi(—tya) + BNi(—iya) 


= 1Z(a) (3-8) 
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The total impedance at a (inward + outward) for which we require an 
equivalent network is 
Z=Z\(a)+ Z 
where Z, is the impedance of the central plug to axial current, viz., 


== nh Io(oa) 
* ~ 2na Iy(ca) 





(3-9) 


To evaluate the constants A and B, the following boundary conditions 
are imposed at radii a and b: 


ata: V = V(a), a given voltage 
at b: V = I(b)Z, 


where Z, is the impedance of the other “short circuit,” comprising the 
outer cylindrical wall. It is given by 
h Ko(ob) 
7,8 3-1 
° 2b Ki(ob) en) 
Except for ignored small deviations of the field around the corners 
of the cavity, the above expressions are exact. The process of finding 
the singularities of Z by successive approximations results in expressions 
that are too long to write down here. To obtain results sufficiently com- 
pact for engineering use, we resort to the following asymptotic approxi- 
mations for the Bessel functions: 


9 \12 
Jo(z) ~ (2) cos (¢ — 7/4) 





Jz) ~ (2) cos (2 — 3/4) 


1/2 


No(z) ~ (2) sin (2 


a/4) (8-11) 


Ni(z) ~ es sin (2 — 37/4) 





Iof2) Kol@) 4 
hh(z) ‘ Ki (2) 
Also, with an error on the order of 107“, 
ho 


Zo(r) OA = Ko(r) = 1/Ho(r) 
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These substitutions result in the following asymptotic formula for the 
total impedance Z at radius a 


2n . me ee 
— cos kx + 7(1+ 4) sin kr 
Z = Kj,(a) DO ee NE (3-12) 


cos ka + 7 sin ka 
No 


where k = ° —landz = —-vya. 


To find an equivalent network of the first kind to represent Z, we deal 
with the admittance, Y = 1/Z. It is instructive and saves much work 
to put Y in the form of exponential functions, with the substitution 


7 — No 
n + no 


which is the reflection coefficient at both inside and outside cylindrical 
surfaces of the cavity. By this means we obtain 





p= 


Y = Hla) 5 (3-13) 


This is now identical in form to the formula (2-13) of example 2, where 
the #-vector was radially, instead of axially, directed. In fact, since 


tke = y(b — a) 


. 1/2 
v ho 


comparison with the similar formulas of example 2 shows that all the 
results of that example can be made to apply to the present one merely 
by changing the dimensional parameters as follows: 


and 


Example 2 Example 3 
(E radial) (EZ axial) 
h goes into b—a 
_ 2ab log (b/a) ' 
d = eh goes into h 
The first result of interest is the value of Q, which is 
1 
Qn = hin ia Ge (3-14) 





— a 


l+¢5 
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where, as before, 6, is the ‘‘skin depth” equal to the real part of o,. That 


is, 
bn as Wong 
/%3 


To gain an idea of numerical magnitudes, consider the same cavity used 
in example 2. The dimensions are, as before, h = .5cm.,b — a = .5 cm. 
For the square cross-section chosen, the first resonance again occurs at 
30,000 megacycles, very nearly, and we can make the following direct 
comparison of the Q’s for the two modes of oscillation: 




















Qn 
n won/2ar 
Ex. 2 (E radial) Ex. 3 (E axial) 
1 380,000 * 108 4250 4370 
2 60,000 x 108 6010 6180 
3 90,000 x 108 7360 7560 











Due to the asymptotic approximations used, the results for example 3 
are not as accurate as those for example 2; the two sets of results show 
only that the Q of the cavity is substantially the same for the two differ- 
ent modes of oscillation. 

The poles of Y are given by 


Pn = An -++ 1B x p11 = Dn = ~"Aa 1B n 
1 1 
= om | ie + ae | (3-15) 


1 1 
om [1+ ae - Gwe 


and the residues are 


I 


Bn 


An = Qn + ibn, An = An = Gn — Abn 


eo Ho(a) won E Se ds tes f | 
7 an 2hé, 2(b — abn 


_ HolGons| 1 1 | 
= en Dh 2 we. 


(3-16) 


Applying formulas (6) gives the following values for the n™ branch of 
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the network of the first kind: 





| rn | 1 
In = Kale) [14+ et oe | 
12. ext ne” 
EO oe E a he =| 
(3-17) 
Gn mn 
Ci 2(b — ajédn 
Ry _ 1 3 
i E aa 2(b — =| 


in all of which wo, = anv/(b — a) and v = 1/(woeo)” = 3(10°) meters 
per second. 

The results can be put in the same form as those obtained for the other 
cavity mode, dealt with in example 2, by employing the “primary con- 
stants” of the cylindrical transmission line, viz.: 


wo) = A [em ]” 140) = 








7a 29 27a 
GG) 26 Cy = si 


In terms of these constants, the elements of the n** branch of the equiva- 
lent network of the first kind are 


_ (6 — a)L(a) if 1 
a 2 (1 + 2hs, a 2(b — an) 


_ (b — a)R(a) 3h 
ee (1+ 5g5) 


_ 24b— aC) (, 3 3 ) 
arr rn? (1 2hbn T 2(b — a)éb,, 


—_ wonC (a) = 3 3 ) 
Bae rns, (1 2Nbn £ 2(b — a)in 





(13-8) 








The network is shown in Fig. 12. 

As in the preceding example, a leakage element arises, in spite of the 
fact that we assumed initially that go of the air in the cavity is zero. 
This element accounts for the losses in the inner and outer cylindrical 
walls. 
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A number of people with whom the above material has been discussed 
have given helpful comments and criticisms. I wish to acknowledge my 
debt in this respect to H. Nyquist, 8. A. Schelkunoff, R. M. Foster, 8. O. 
Rice, J. Riordan and W. H. Wise. 
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October, 1951, Part II, entitled The DT-2 Microwave Radio System 
reprinted as Monograph 1921). 

To meet the stringent requirements of the 4,000-mile transcontinental micro- 
wave relay system, a number of new developments had to be included in the 
design of the repeater stations. The circuits of these uuattended stations, and 
how they are maintained, are the subject of this article. 
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p-m or f-m mobile telephone transmitters, such as those used in 30 to 44 -me 
highway and 152 to 175-mc urban service. Measures r-f power output, audio 
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D. A. Auspere’. Proc. I.R.E., 39, pp. 1393-1400, Nov., 1951. (Mono- 
graph 1911). 

The impedance of a two-terminal network is defined completely by the inser- 
tion loss and phase shift it produces when inserted between known sending and 
receiving impedances. Recent advances in precise wide-band phase and trans- 
mission measuring circuits have permitted practical use of this principle. Reac- 
tive and resistive impedance components are read directly from a simple graphi- 
cal chart in which frequency is not a parameter. The basic principle described 
promises attractive possibilities in many cases of impedance measurements 
where present methods are inadequate. 


Electron-Vibration Interactions and Superconductivity. J. BARDEEN’. 
Revs. Modern Phys., 23, pp. 261-270, July, 1951. (Monograph 1912). 


The Copper Oxide Rectifier. W. TH. Brarrain’. Revs. Modern Phys., 


23, pp. 203-212, July, 1951. 
It is shown that the conductivity in the ohmic part of the cuprous oxide layer 
can be explained with the usual band picture of semiconductors only by assum- 
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ing the presence of some donor-type impurities in addition to the usual acceptor 
type. The energy difference between the acceptors and the filled band is 0.3 
electron volt, and the total number of impurity atoms is about 10“ to 10'* per 
em, the number of donors being less than but of the same order as the number of 
acceptors. One finds that the density of ion charge in the rectifying layer is of 
the same order of magnitude as the difference between the donors and acceptors 
found from the conductivity. The field at the copper-cuprous oxide interface is 
about 2 x 10 volts/cm; the height of the potential at the surface as compared 
with the oxide interior is about 0.5 volt; and the thickness of the space charge 
layer about 5.0 x 10-5 cm. The diffusion equation for flow of current through this 
space charge region can be integrated to give the current in terms of the field 
at the interface and the applied potential across the space charge layer. Two 
currents are involved, one from the semi-conductor to the metal (J,) and one 
from the metal to the semiconductor (Jm) which is similar to a thermionic emis- 
sion current into the semiconductor. The net current is, of course, J = In — Is. 
One can get this “emission” current (J) by dividing the true current by the 
factor 1 - exp (—eV./kT), where V, is the applied potential. This emission current 
depends on the absolute temperature and on the field at the copper-cuprous 
oxide interface. At high fields the logarithm of the current is proportional to the 
square root of the field, and at low fields the current decreases more rapidly 
indicating a patchy surface having small areas of low potential maximum from 
which all the emission comes when the field is large. 


Effect of Packaging on Corrosion of Zinc Plated Equipment. K. G. 
Compron’, 8. M. Arnot’, and A. Menpizza’. Corrosion, 7, pp. 365- 
372, Nov., 1951. 


Physics as a Science and an Art. K. K. Darrow’. Phys. Today, 4, pp. 


6-11, Nov., 1951. (Monograph 1914). 

The last of six invited papers presented on October 25th during the symposium 
on ‘physics today’ which keynoted the 20th Anniversary Meeting of the 
American Institute of Physics in Chicago. 


Ionization by Electron Impact in CO, N2, NO, and Oz. H. D. Hacstrum’. 
Revs. Modern Phys., 23, pp. 185-208, July, 1951. (Monograph 1916). 


Tonization by electron impact in diatomic gases has been studied in this work 
with a mass spectrometer designed to measure m/e, appearance potential, and 
initial kinetic energy for each ion observed. Results have been obtained for the 
gases CO, N2 NO, and O: with some confirmatory work in H;. Discussion is 
included of the nature and identification of dissociative ionization processes and 
of the retarding potential and appearance potential measurements. Values of 
important quantities such as the dissociation energies of CO, Nz, and NO; the 
sublimation energy of C; the electron affinity of O; and the excitation energy of 
O~ are determined again by electron impact in this work. 
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Equivalent Temperature of an Electron Beam. M. E. Hings’., Letter 
to the Editor., J. Appl. Phys., 22, pp. 13885-1386, Nov., 1951. 


Bell System Cable Sheath Problems and Designs. F. W. Horn! and 
R. B. Ramsey’., Elec. Eng., 10, pp. 1070-1075, Dec., 1951. (Monograph 
1917). 


Engineering Planning. H. S. Osporne’. J. Eng. Educ., 42, pp. 121- 
125, Nov., 1951. 


Acceptance Inspection of Purchased Material. J. &. Paumer® and E. G. 
D. Parerson’. Ind. Quality Control, 8, pp. 23-27, Nov., 1951. 


A Note on the Partial Differential Equations Describing Steady Current 
Flow in Intrinsic Semiconductors. R. C. Prim’. Letter to the Editor. 
J. Appl. Phys., 22, pp. 13888-1389, Nov., 1951. 


General Theory of Symmetric Biconical Antennas. 8. A. SCHELKUNOFF'. 
J. Appl. Phys., 22, pp. 1830-1332, Nov., 1951. (Monograph 1922). 


This paper presents the input admittance of a symmetric biconical antenna 
of an arbitrary angle as the limit of a certain sequence of functions. The first 
term of this sequence approaches the exact expressions for the input admittance 
as the cone angle approaches either zero or 90°. For this reason our conjecture is 
that this term represents a good first approximation for all angles. 


Artificial Dielectrics for Microwaves. W. M. Suarpuess’. Proc. I. R. E., 
39, pp. 1889-1393, Nov., 1951. (Monograph 1923). 

This paper presents a procedure for measuring the dielectric properties of 
metal-loaded artificial dielectrics in the microwave region by the use of the 
short-circuited line method. Formulas, based on transmission-line theory, are 
included and serve as guides in predicting the approximate dielectric properties 
of certain loading configurations. 
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